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ABSTRACT 


Recent  surveys  of  actual  subgrade  moisture  conditions 
beneath  impervious  highway  and  airport  pavements  show  that 
many  subgrades  do  not  attain  complete  saturation  in  the 
course  of  time,  but  tend  towards  a  particular  terminal  value. 

If  this  terminal  moisture  content  of  a  subgrade  under  a 
given  set  of  conditions  could  be  predicted,  the  design  of 
the  pavement  could  be  based  on  the  strength  of  the  soil  at 
this  moisture  content.  For  making  such  a  prediction,  an 
understanding  of  the  various  mechanisms  by  which  water  may 
move  through  soils  is  necessary. 

The  case  of  flow  through  fully  saturated  soils  under 
pressure  potentials  is  relatively  well-understood,  and  the 
solutions  of  problems  in  this  category  are  based  on  the 
empirical  Darcy's  law.  On  the  other  hand,  the  case  of  flow 
through  partially  saturated  soils  presents  many  difficulties, 
and  various  theories  have  been  proposed  by  different 
investigators  to  account  for  the  phenomenon;  these  theories, 
all  of  them  based  on  some  simplifying  assumptions,  are 
reviewed  in  some  detail.  Of  particular  interest  to  the 
problem  on  hand  is  the  fact  that  the  theory  based  on  the 
concept  of  suction  pressure  of  soil  water  has  been  successfully 
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applied  to  the  case  of  water  accumulation  in  subgrades  in 
England. 

In  general,  electro-osmosis  may  not  contribute  to  the 
accumulation  of  water  beneath  pavements?  however,  a  brief 
review  of  the  basic  principles  involved  is  included  in  this 
report,  since  this  is  one  important  mechanism  of  moisture 
migration,  and  since  some  investigators  believe  that 
thermo-osmosis  is  only  a  special  form  of  electro-osmosis. 

In  the  field  of  thermo -osmosis,  there  exists  a  great 
deal  of  controversy  as  to  the  nature  of  the  mechanism 
involved.  Our  knowledge  of  the  overall  subject  is  still 
immature  and  so  a  complete  and  accurate  picture  of  the 
phenomenon  cannot  as  yet  be  given.  The  available  literature 
on  the  subject  is  reviewed  at  some  considerable  detail,  and 
the  present  state  of  our  knowledge  is  summarized. 

The  experiments  conducted  in  connection  with  this 
thesis  were  aimed  at  verifying  the  existence  of  a  return 
flow  in  the  liquid  phase  of  the  water  accumulated  at  the  cold 
end  of  specimens  due  to  thermo-osmosis ;  such  a  return  flow 
has  been  assumed  by  Hutcheon  in  his  hypothesis  for  the 
phenomenon . 

Silica  sand,  all  passing  a  60-mesh  sieve,  was  mixed  with 
different  amounts  of  distilled  water  to  give  six  different 
initial  moisture  contents,  and  was  compacted  into  specimens 
inside  Lucite  moulds  16  centimetres  long  and  4U4 
centimetres  diameter,  employing  a  suitable  technique  to  give 
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as  uniform  and  consistent  a  density  as  possible.  In  all, 
thirty  such  specimens  -  five  each  for  each  moisture  content  - 
were  prepared  and  subjected  to  a  temperature  gradient  of 
4.5  degrees  Fahrenheit  per  centimeter  and  twenty  four  more- 
four  each  for  each  moisture  content  -  were  subjected  to  a 
gradient  of  2.25  degrees  Fahrenheit  per  centimeter.  The  first 
specimen  of  each  group  was  sampled  after  an  elapsed  time  of 
77  hours,  and  the  subsequent  specimens  were  sampled  at 
intervals  of  24  hours  thereafter.  The  specimens  were 
sectioned,  and  moisture  contents  and  dry  densities  were 
determined  for  every  2  centimeter  section  of  the  specimens. 

The  experimental  data  do  reveal  clearly  that  the  water 
accumulated  at  the  cold  end  tends  to  redistribute  itself  in 
the  liquid  phase  as  far  as  conducting  water  films  are  present. 

In  specimens  having  a  fairly  high  initial  degree  of 
saturation,  this  redistribution  extends  as  far  as  the  hot 
end  and  thus  a  state  of  dynamic  equilibrium  may  not  be  reached 
at  all.  The  data  also  reveal  that  there  is  an  optimum  initial 
degree  of  saturation  for  maximum  thermo -osmotic  transfer,  but 
no  exact  quantitative  measure  of  this  value  could  be  obtained 
since  the  system  had  not  attained  a  state  of  dynamic  equilibrium. 
For  the  soil  tested,  this  value  was  found  to  range  from  15  to 
20  per  cent  saturation  for  the  higher  temperature  gradient, 
and  from  13.5  to  17.2  per  cent  saturation  for  the  smaller 
temperature  gradient. 

Recommendations  for  further  investigation  of  the  topic 


are  included. 
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CHAPTER  I 


INTRODUCTION  OF  THE  PROBLEM 


1.1.  Statement  of  the  problem: 

"The  accumulation  of  water  in  subgrades  and  bases 
(beneath  impervious  highway  and  airport  pavements)  ...  is 
a  fact  often  observed  and  always  dreaded  in  airport  and 
highway  engineering.  This  water  accumulation  is  not 
restricted  to  regions  of  humid  and  subhumid  climates,  nor 
to  the  presence  of  high  ground  water  tables?  rather,  some 
of  its  most  impressive  examples  have  been  observed  in 
deserts  and  at  considerable  distances  from  the  ground 
water  table"  (Winterkorn  and  Eyring,  1945) .  Such 
accumulation  of  moisture  beneath  pavements  may  cause  a 
large  reduction  in  the  bearing  capacity  of  the  subgrade, 
thus  leading  to  the  failure  of  the  pavement. 

In  the  past,  highway  and  airport  engineers  have 
leaned  towards  the  assumption  that  all  subgrades  will 
eventually  become  saturated  due  to  capillarity  and  that 
the  pavement  should  be  designed  to  withstand  this  extremely 
unfavourable  subgrade  condition  (Spangler,  1945)  .  However, 
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recent  subgrade  moisture  surveys  of  existing  highway  and 
airport  pavements  in  U.S.A.  have  shown  that  many  subgrades 
do  not  approach  this  limiting  condition  of  complete 
saturation;  rather ,  they  seem  to  approach  a  fairly  constant 
terminal  value  which  is  appreciably  less  than  saturation 
(Kresten  1944,  1945) . 

Naturally  the  question  then  arises:  what  will  be 
the  final  equilibrium  moisture  distribution  beneath  a  given 
pavement  under  specified  boundary  and  initial  conditions, 
and  how  can  this  be  determined?  Broadly,  there  are  three 
ways  of  seeking  an  answer  to  this  important  question.  One 
may  try  to  view  the  problem  from  the  standpoints  of  accepted 
classical  laws  of  hydrodynamics  or  of  thermodynamics  and 
arrive  at  a  final  equilibrium  condition  compatible  with  the 
nature  of  the  problem.  This  may  be  called  the  classical 
approach.  All  known  forces  that  may  influence  the  final 
state  are  considered  simultaneously,  and  the  condition  of 
dynamic  equilibrium  between  these  forces  governs  the  final 
state  of  the  soil  beneath  the  pavement. 

Alternatively,  one  may  consider  each  force  on  its 
own  merits  and  obtain  the  final  equilibrium  condition  by 
superimposing  the  results  of  all  these  studies.  This  is 
the  phenomenological  approach  to  the  problem.  Although 
less  general  than  the  classical  approach,  this  method 
permits  one  to  obtain  a  better  understanding  of  the  problem. 
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In  the  third  method,  one  may  list  the  various  sources 
from  which  water  may  enter  the  subgrade  and  try  to  evaluate 
the  magnitudes  of  the  quantities  involved.  Moisture  may 
find  access  into  the  subgrade  through  the  following 
sources  :  i)  rain  water  may  percolate  downward  through  cracks 
in  the  pavement;  ii)  water  from  the  shoulders  may  be  drawn 
into  the  subgrade  laterally;  iii)  water  may  be  sucked 
upwards  from  the  water  table  by  capillarity;  iv)  water  may 
diffuse  upwards  in  the  vapour  phase;  v)  there  may  be 
gravitational  seepage  of  water  from  the  adjoining  high  grounds; 
and,  vi)  the  water  table  itself  may  rise  and  reach  the  subgrade 
under  favourable  conditions  (Low  and  Lovell,  1959) .  A 
complete  and  precise  evaluation  of  the  contribution  of  these 
various  sources  appears  to  be  a  formidable  task.  To  the 
writer's  knowledge,  this  approach  has  not  so  far  been 
employed  in  practice. 

1.1.2.  Application  in  other  fields: 

It  is  perhaps  fortunate  that  the  highway  and  airport 
engineers  are  not  the  only  ones  faced  with  this  problem;  the 
accumulation  of  water  in  subgrades  is  only  a  part  of  the 
overall  problem  of  flow  of  water  through  soils,  or,  in  more 
general  terms,  of  flow  of  fluids  through  porous  media.  These 
problems  have  very  important  applications  in  many  different 
fields  of  engineering  and  science.  It  will  not  be  possible 
to  list  all  the  various  fields  of  application;  only  a  few  of 
the  more  important  cases  are  mentioned  below. 
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In  civil  engineering  practice,  seepage  of  water 
through  and  beneath  earth  dams  and  similar  water  retaining 
structures,  construction  of  filter  beds  for  municipal  water 
supply  systems,  the  general  subject  of  ground  water 
hydrology,  the  rates  of  settlement  of  structures  founded  on 
compressible  soils,  the  control  of  destructive  frost  heaving 
of  pavements,  and  many  other  similar  cases  may  be  cited, 
wherein  the  basic  problem  involves  the  flow  of  water  through 
soils . 

In  petroleum  engineering  practice,  the  yield  of  oil 
wells  obviously  depends  upon  the  ease  with  which  natural  gas 
and  oil  can  flow  through  rocks  and  soils.  In  ceramic 
engineering,  electrophoresis  (reverse  process  of  electro¬ 
osmosis;  see  Chapter  III)  is  employed  for  obtaining  pure 
kaolinite  from  dilute  suspensions,  and  this  involves 
dewatering  of  the  kaolin  clay  by  means  of  applied  electric 
potentials  (Speil  and  Thompson  1942) .  In  foundry  engineering 
practice,  flow  of  gases  through  moulding  sands  is  of  great 
importance . 

The  importance  of  the  phenomenon  in  the  realms  of 
agriculture  and  soil  conservation  is  too  well  known  to  warrant 
special  mention  here. 

1.2.  Various  efforts  to  solve  the  problem. 

1.2.1.  Classical  hydrodynamics  applied  to  soil-water  systems: 

Since  the  system  involves  the  dynamic  interaction 
between  a  fluid  and  a  (porous)  solid,  it  must  obey  the  classical 
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laws  of  hydrodynamics  -  the  law  of  conservation  of  matter  and 
the  law  of  viscous  flow.  The  law  of  viscous  flow  has  been 
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formulated  in  hydrodynamical  terms  as  a  set  of  three  differential 
equations  known  as  the  Navier-Stokes 1  equations  (see  Appendix 
A . ,  Vallentine  1959,  Muskat  1937) .  These  equations,  in 
conjunction  with  the  law  of  conservation  of  matter  and  an 
equation  describing  the  state  of  the  fluid  (liquid  or  gas)  and 
the  dynamic  nature  of  the  flow  (adiabatic  or  isothermal) , 
theoretically  suffice  to  solve  for  the  five  unknown  variables 
in  any  given  problem,  namely,  the  fluid  pressure,  the  fluid 
density,  and  the  three  components  of  the  fluid  velocity  in 
the  three  co-ordinate  directions. 

In  practice,  however,  the  Navier-Stokes1  equations  are 
so  complicated  in  their  formulation,  that  a  complete  solution 
is  not  possible  except  in  special  cases  with  very  simple 
boundary  conditions  (Vallentine,  1959) .  When  trying  to  apply 
these  laws  to  the  soil-water  systems,  one  is  faced  with  the 
following  complicating  factors:  (1)  in  most  cases,  the 
phenomenon  involves  the  interaction  of  three  different  phases 
of  matter  -  solid,  liquid  and  gas;  (2)  there  are  a  large 
number  of  distinct  force  potentials  that  may  cause  the  flow, 
and  the  extent  to  which  each  of  these  mechanisms  contribute 
to  the  mass  flow  is  unknown;  (3)  the  individual  pore  channels 
are  so  tortuous  and  their  shapes  and  cross-sectional  areas 
are  so  variable  in  space,  that  an  adequate  geometric 
description  of  the  boundaries  of  the  flow  region  is  impossible; 
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and  (4)  chemical  and  electrochemical  reaction  between  the 
percolating  water  and  the  colloidal  fraction  of  the  soil, 
growth  of  micro-organisms  in  the  pores  of  the  soil,  and 
similar  factors  exert  some  influence  on  the  characteristics 
of  the  flow  region. 

With  so  many  complications  involved,  it  is  not 
surprising  that  Darcy  resorted  to  the  experimental  approach, 
when  he  was  interested  in  the  flow  of  water  through  filter 
sands,  and  observed  the  exceedingly  simple  relationship:  the 
filter  velocity  is  proportional  to  the  first  power  of  the 
hydraulic  gradient  causing  flow.  After  much  experimental 
verifications  by  numerous  workers  in  the  various  fields, 
Darcy's  empirical  law  has  now  become  to  be  accepted  as  valid 
over  a  wide  range  of  conditions  of  flow,  and  not  merely  for 
the  case  of  flow  of  water  through  filter  sands. 

The  inability  to  apply  the  classical  laws  of 
hydrodynamics  to  the  problem  has  prompted  a  few  investigators 
to  accept  Darcy's  law  as  the  macroscopic  equation  of  motion-in 
place  of  the  classical  Navier-Stokes 1  equations-and  combine 
this  equation  with  the  physical  equations  of  continuity  and  of 
state,  in  order  to  obtain  a  final  mathematical  equation 
describing  the  flow  (see  Appendix  A,  also  Muskat  1937, 
Scheidegger  1957) .  At  present,  this  semiempirical  approach 
is  widely  used  in  all  problems  involving  flow  of  liquids 
through  saturated  porous  media. 
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1.2.2.  Classical  thermodynamics  applied  to  soil-water  systems: 

The  other  classical  approach  to  solve  the  problem 
stems  from  the  principles  of  thermodynamics,  which  is  concerned 
with  the  various  forms  and  transformations  of  energy;  while 
classical  hydrodynamics  is  concerned  with  the  potential 
gradients  and  the  principles  of  conservation  of  matter,  this 
new  approach  is  concerned  with  energy  gradients  and  the 
principle  of  conservation  of  energy. 

It  is  generally  conceded  that  Buckingham  (1907) ,  with 
his  introduction  of  the  term  ^capillary  potential5,  was  the 
first  to  recognize  the  possibility  of  solving  the  problem 
from  considerations  of  the  energy  state  of  the  soil  water. 
Subsequently,  this  energy  concept  of  soil  water  has  been 
utilized  by  numerous  workers  in  soil  physics  and  soil  mechanics 
(Gardner  1919,  Schofield  1935  ,  Russel  and  Spangler  1941, 

Croney  1952,  and  many  others) .  In  particular,  the  energy 
concept  has  been  successfully  applied  to  the  problem  of 
predicting  the  equilibrium  moisture  distribution  beneath 
impervious  pavements  in  England  (Coleman  1959) „  Further 
research  aimed  at  evolving  a  generalized  method  of  analysis 
that  will  be  applicable  to  cases  of  polyphase  migration  of 
water-liquid  flow,  vapour  flow  and  film  flow  -  is  reportedly 
under  way  in  the  Road  Research  Laboratory,  England  (Croney  1958) 

While  conceding  the  fact  that  the  thermodynamical 
principles  are  applicable  to  the  soil-water  system,  it  is  well 
to  point  out  certain  factors  that  must  be  clarified  before  a 
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unified  theory  could  be  evolved.  In  the  first  place ,  one 
must  realize  that  the  interaction  between  soil  and  water  is 
not  a  thermodynamically  reversible  process.  Such  a 
realization  is  not  lacking,  but  the  same  authors,  who 
explicitly  state  that  there  are  hysteresis  effects,  assume 
a  single-valued  relationship  between  the  moisture  content 
of  the  soil  and  the  pressure  deficiency  (suction)  in  the 
soil-water  (see  Section  2.7.) .  In  the  second  place,  the 
thermodynamic  nature  of  the  adsorbed  films  of  water  needs 
further  clarification;  it  is  generally  stated  that  the 
properties  of  the  water  in  these  films  are  different  from 
those  of  the  bulk  water  outside  the  zone  of  influence  of 
the  adsorptive  force  field.  It  is  not  definitely  known  -  as 
far  as  the  writer  is  aware  -  whether  the  energy  state  of  these 
films  is  similar  to  that  of  free  water,  except  for  the 
modification  necessitated  by  the  adsorptive  force  field. 

Finally,  the  energy  relationships  of  three-phase  interfaces 
are  still  obscure.  The  writer  feels  that  these  and  other 
related  puzzles  must  be  satisfactorily  explained  and  properly 
accounted  for,  before  we  can  claim  to  solve  the  problem. 


"Within  the  adsorbed  layer  the  physical  properties  of  the 
water  are  very  different  from  those  of  free  or  normal 
water  at  the  same  temperature.  Near  the  surface  of  the 
particle  the  water  has  the  properties  of  a  solid.  At  the 
middle  of  the  layer  it  resembles  a  very  viscous  liquid. 

As  the  outer  surface  is  approached,  the  properties  of  the 
water  becomes  normal."  (Terzaghi  and  Peck,  pp.  11,  1958) . 
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1.2.3.  Phenomenological  approach  to  soil-water  systems: 

The  first  of  these  phenomenological  studies  is  the 
experiment  of  Darcy  in  1856 ,  on  the  flow  of  water  through 
filter  sands  under  hydrostatic  pressure  potentials.  This 
was  followed  by  Helmholtz's  mathematical  analysis  of  the 
electro-kinetic  phenomenon  at  solid-liquid  interfaces,  that 
led  to  the  analysis  of  flow  of  water  through  fine  grained 
soils  under  electric  potentials.  Later  (1907),  Lippman 
found  that  liquids  could  be  made  to  move  through  semipermeable 
membranes  under  the  influence  of  applied  thermal  potential 
forces,  and  he  coined  the  term  1thermo-osmosist  to  describe 
the  phenomenon  (Hutchison  et  al .  1948) .  Although  numerous 
experiments  have  since  been  conducted  -  either  as  a  fundamental 
research  on  thermo-osmosis  itself,  or  as  an  inevitable  part 
of  some  other  research  problem  -  the  phenomenon  is  still 
defying  an  adequate  theoretical  treatment.  The  numerous 
hypotheses  that  have  been  put  forward  by  the  various  workers 
are  briefly  discussed  elsewhere  in  this  dissertation. 

Besides  these  major  divisions  of  the  overall  phenomenon, 
there  may  be  many  other  specific  cases  of  flow.  For  example, 
the  transpiration  of  water  by  plants  involves  differential 
osmotic  pressure  exerted  by  the  roots;  whether  a  similar 
effect  of  considerable  significance  is  present  in  the  case  of 
flow  of  water  through  soils  (from  an  engineering  standpoint) 
has  not  been  systematically  investigated.  Such  physico-chemical 
properties  as  the  differential  energy  of  hydration  of  the  ions. 
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the  differential  heat  of  wetting  and  of  adsorption,  and  many 
others  may  have  significant  influence  on  the  movement  of 
water  through  soils.  As  our  knowledge  of  the  physicochemistry 
of  the  soil-water  systems  is  still  incomplete,  these  factors 
have  not  yet  been  adequately  investigated. 

103.  Concluding  remarks: 

A  phenomenological  approach  has  been  employed  in  this 
dissertation  also,  since  the  classical  approaches  are  far  from 
being  complete.  The  flow  of  water  through  soils  under  pressure 
potentials  is  described  in  the  next  chapter;  the  subjects  of 
electro-osmosis  and  of  thermo -osmosis  are  discussed  in  the 
succeeding  two  chapters;  the  basic  principles  involved  in  the 
classical  methods  of  hydrodynamics  and  of  thermodynamics  are 
briefly  outlined  in  Appendices  A  and  C  to  this  dissertation. 

As  the  experimental  work  performed  in  connection  with  this 
thesis  is  in  the  field  of  thermo-osmosis,  a  rather  detailed 
review  of  the  existing  literature  on  this  subject  is  included. 
The  writer  makes  no  similar  claim  in  the  case  of  other 


mechanisms  of  water  movement. 
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CHAPTER  II 


FLOW  OF  WATER  UNDER  PRESSURE  POTENTIALS 


2.1.  Introduction: 

In  this  chapter,  the  flow  of  water  through  soils 
that  is  initiated  by  fluid  pressure  differentials  will  be 
considered.  The  subject  matter  divides  itself  into  two 
distinct  parts  -  one  dealing  with  flow  through  fully 
saturated  soils  the  other  dealing  with  flow  through  partially 
saturated  soils.  Many  of  the  theoretical  analyses  described 
in  the  following  sections  are  based  on  the  well-known  Darcy's 
law  and  so,  a  brief  review  of  the  law  itself  is  given  first. 

2.1.1.  Darcy's  law: 

Based  on  the  results  of  his  experiments  on  the  flow 
of  water  through  filter  sands,  Darcy  gave  the  empirical  equation 

q  =  v . A  =  k.A.h/L  (1) 

where,  Q  =  the  discharge  through  the  filter  bed  in  unit  time, 
v  =  the  filter  velocity  or  the  discharge  velocity, 

A  =  the  area  of  cross-section  of  the  filter  bed  in  a 
direction  normal  to  the  overall  direction  of  flow, 

As  the  subject  of  the  dissertation  is  primarily  concerned 
with  the  latter,  only  a  brief  outline  of  the  former  case 
is  given. 
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k  =  the  coefficient  of  permeability  or  the  permeability 1 
of  the  bed, 

h  =  the  hydraulic  head  lost  in  flow  through  the  filter,  and 
L  =  the  length  of  the  filter  bed. 

Although  Darcy's  experiments  were  concerned  with  uniform 
sands,  the  equation  itself  has  been  verified  by  numerous 
investigators  for  widely  different  soils  with  varying  degrees 
of  uniformity.  As  a  result  of  the  work  of  these  investigators 
it  is  now  known  that  the  empirical  Darcy's  law  is  valid  over 
a  considerable  range  of  practical  interest,  subject  to  the 
limitations  discussed  in  a  subsequent  section  of  this 
dissertation  (cf.  Section  2.3.). 

2.1.2.  Units  of  permeability: 

From  a  dimensional  analysis  of  equation  (1)  ,  it  is 
seen  that  the  coefficient  of  permeability  has  the  units  of 
velocity.  This  is  the  way  in  which  the  coefficient  is  used 
in  civil  engineering  practice.  Accordingly,  the  coefficient 
of  permeability  is  defined  as  the  average  filter  velocity 
through  a  porous  medium,  when  the  hydraulic  gradient  (h/L)  is 
unity.  It  is  obvious  that  such  a  definition  renders  the 
permeability  dependent  on  the  properties  of  the  porous  medium 
as  well  as  upon  those  of  the  percolating  fluid.  When  one 
deals  with  many  different  fluids,  it  is  necessary  to 
dissociate  the  two  influences  -  of  the  properties  of  the 
porous  medium  and  of  those  of  the  fluid  -  on  the  permeability, 
so  as  to  make  the  permeability  dependent  solely  upon  the 
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characteristics  of  the  porous  medium.  Nutting  (1930)  was  the 
first  to  set 

k  =  K  V\i  U) 

where,  p.  =  the  viscosity  of  the  fluid,  and 

K1  =  the  1  specific  permeability1  of  the  porous  medium. 
Incorporating  this  into  Darcy Ss  law,  equation  (1)  becomes: 

Q  =  KJ/p..A.h/L  (3) 

It  is  well-known  that  when  two  ends  of  a  continuous 
column  of  liquid  are  maintained  at  two  different  levels,  the 
pressure  difference  (^P)  between  the  two  ends  is  related  to 
the  hydraulic  head  difference  (h)  by  the  expression 

^  P  =  h.y  (4-) 

where  y  is  the  unit  weight  of  the  liquid.  Since  the  unit 
weight  of  the  liquid  is  again  a  property  of  the  liquid,  and 
since  the  phrase  "hydraulic  head1  may  suggest  that  the 
percolating  liquid  is  water,  it  is  preferable  to  formulate 
the  equation  in  terms  of  the  pressure  difference,  instead  of 
hydraulic  head  difference;  thus  we  let 

K*  =  K.y  (5) 

and  Eq.  (3)  becomes 

Q  =  K/|i.A.  (aP)/l  (S) 

It  will  be  noticed  the  final  coefficient,  K,  in  Eqa 
(6)  has  the  units  of  a  distance  squared,  and  that  it  is 
independent  of  the  physical  properties  of  the  percolating 
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2.2.  Limitation  of  Darcy's  law: 

Since  Darcy's  law  is  based  on  experimental  observations 
it  is  to  be  expected  that  it  will  be  valid  over  a  certain 
range  of  conditions  only.  As  the  law  has  come  to  be 
recognized  as  the  basis  of  all  analyses  of  flow  problems 
initiated  by  pressure  potentials ,  the  range  of  validity  of 
the  law  is  of  particular  interest.  To  begin  with,  since  the 
average  filter  velocity  is  set  as  proportional  to  the  first 
power  of  the  pressure  gradient  causing  flow,  the  assumption 
that  the  friction  losses  are  proportional  to  the  steady-state 
velocity  is  implied.  In  the  realm  of  hydraulics,  such  a  flow 
is  implied.  In  the  realm  of  hydraulics,  such  a  flow  is 
referred  to  as  'laminar  flow1.  Thus,  Darcy's  law  is  valid 
for  cases  of  steady,  laminar  flow  only.  The  common  way  of 
describing  this  type  of  flow  is  by  means  of  the  Reynolds 
number,  which  is  the  ratio  of  the  inertia  forces  to  the 
viscous  forces;  in  the  case  of  pipes  flowing  full,  the 
Reynolds  number  is  2000  or  less  if  the  flow  is  laminar. 

Attempts  to  establish  a  similar  critical  value  for  the 
cases  of  flow  through  soils  have  proved  to  be  futile,  since 
the  'critical*  value  is  found  to  vary  erratically  between 
0.1  and  75  (Scheidegger  1957) .  Now  it  is  generally  accepted 
that  the  flow  through  soils  is  essentially  laminar  if  the 
Reynolds  number  is  less  than  unity.  Incidentally,  the  vast 
difference  between  the  critical  Reynolds  numbers  for  the 
case  of  flow  through  pipes  and  that  of  flow  through  soils, 
indicates  that  the  two  flow  systems  are  not  dynamically  similar 
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2.2.2.  Molecular  effects : 

Darcy's  law,  as  well  as  Poiseuille's  law,  assume 
full  viscous  resistence  to  flow  and  so  the  velocity  of 
flow  must  be  zero  at  the  boundaries  of  the  flow  region, 
namely,  the  walls  of  the  conducting  medium.  Actually, 
due  to  the  molecular  nature  of  the  percolating  fluid, 
there  is  a  “hopping  of  the  molecule1  effect  at  the  walls 
of  the  container.  This  gives  rise  to  “slip  flow1  (also 
called  Knudsen  flow)  and  is  more  pronounced  in  the  case  of 
flow  of  gases,  when  the  distance  between  the  walls  of  the 
container  is  of  the  same  order  of  magnitude  as  the  molecular 
free  path  of  the  flowing  gas.  Further,  the  permeability  of 
the  same  medium  (in  length  squared  units)  for  flow  of 
liquids  and  of  gases  are  different,  with  the  liquid 
permeabilities  being  lower  (Scheidegger  1957) . 

2.2.3.  Other  effects: 

Besides  increased  friction  losses  resulting  from  high 
flow  rates  and  increased  quantities  of  flow  resulting  from 
molecular  slip,  there  are  some  more  anomalies  that  have  been 
observed  by  numerous  workers.  One  such  anomaly  results  from 
the  boundary  effect,  which  causes  a  reduction  in  the  velocity 
of  the  fluid  stream  near  the  walls  of  the  container  and  near 
the  center  of  the  bed,  with  the  maximum  velocity  somewhere 
in  between;  this  is  at  variance  with  the  assumption  of  an 
'average  uniform  filter  velocity'.  Another  series  of  anomalies 


are  ascribed  to  ionic  effects,  wherein  the  permeability  of  the 
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same  medium  is  found  to  vary  with  the  acidity  or  alkalinity 
of  the  percolating  fluid  (Quirk  and  Schofield  1955) .  Finally, 
as  explained  in  a  subsequent  chapter  of  this  dissertation, 
whenever  a  fluid  flows  through  a  capillary  passage,  an 
electric  potential  called  the  streaming  potential  will  be 
generated  across  the  ends  of  the  capillary;  under  the 
influence  of  this  induced  electric  potential,  a  small 
quantity  of  fluid  may  flow  in  the  reverse  direction,  closely 
analogous  to  electro-osmotic  flow. 

Thus,  Darcy 11  s  law  breaks  down  whenever  the  flow  rate 
is  very  high,  or  when  molecular  slip  flow  is  present,  or 
when  the  pores  are  small  enough  to  give  rise  to  streaming 
potentials;  above  all,  the  law  is  valid  only  under  conditions 
of  steady  state  laminar  flow;  the  pressure  gradient  as  well 
as  the  pore  geometry  must  be  independent  of  time.  In  fissured 
clays,  percolation  is  almost  always  associated  with  the 
leaching  of  the  fines  and  these  fine  particles  may  be 
redeposited  in  a  different  place;  also,  entrapped  air  may  go 
into  or  come  out  of  solution.  These  and  similar  factors 
give  rise  to  unpredictable  changes  in  the  pore  geometry.  In 
spite  of  these  limitations,  the  range  of  validity  of  the  law 
is  sufficiently  large  to  include  many  problems  of  practical 


interest . 
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2.3.  Flow  through  saturated  soils. 

2.3.1.  Introduction: 

Natural  masses  of  soils  well  below  the  range  of 
fluctuation  of  the  ground  water  table  are  essentially 
saturated.  Although  it  is  probable  that  the  zone  over 
which  the  water  table  fluctuates  is  not  saturated,  it  is 
commonly  assumed  that  such  zones  are  also  fully  saturated. 

In  a  similar  manner,  the  soil  beneath  an  impounding 
reservoir  is  also  assumed  to  be  fully  saturated.  Thus, 
the  seepage  of  water  beneath  dams  is  a  problem  in  flow 
through  saturated  soils;  the  flow  of  water  into  open  cuts 
and  trenches  that  extend  to  a  level  below  that  of  the  ground 
water  table  and  the  yield  from  natural  wells  are  problems 
of  the  same  category;  so  is  filtration. 

2.3.2.  Application  of  Darcy's  law: 

In  general,  the  problem  may  involve  flow  in  one 
direction  only  (filtration) ,  two  directions  (seepage  under 
long  dams)  or  three  directions  (flow  into  a  well  that  does 
not  completely  penetrate  the  aquifer)  .  Uni-dimensional  and 
two-dimensional  flow  systems  are  relatively  easy  to  solve, 
and  the  various  methods  of  solution  are  briefly  outlined  in 
Appendix  A  of  this  dissertation.  Solutions  of  three- 
dimensional  flow  systems  are  extremely  complicated,  and 
comparatively  only  few  specific  solutions  have  been  obtained 
for  particular  problems  (cf.  Muskat,  1937) .  Analytical 
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solutions  are  the  most  convenient  of  all  methods  for  certain 
classes  of  uni-dimensional  flow  problems;  two-dimensional 
flow  problems  are  somewhat  difficult  to  solve  analytically 
(except  under  favourable  boundary  and  initial  conditions) 
and  recourse  may  have  to  be  taken  to  analog  solutions  and 
graphical  solutions  (Schmid  1958 ,  Casagrande,  A.  1937) . 

2.3.3.  Influence  of  capillary  fringe: 

It  is  well  known  that,  in  most  cases,  the  region  just 
above  a  water  table  is  essentially  saturated  with  water  (see 
Terzaghi  1959,  pp.  303),  especially  so,  if  the  water  table 
has  dropped  to  its  present  level  by  the  process  of  drainage. 
Experiments  with  electric  analogs  by  Childs  (1945)  and  with 
mechanical  models  by  Swertzendruber  and  Kirkham  (1956) 
indicate  that  this  capillary  fringe  (saturated  portion  above 
the  water  table)  may  increase  or  decrease  the  flow  through 
the  soil,  depending  upon  the  particular  conditions.  Thus, 
if  the  flow  is  essentially  horizontal,  capillary  fringe  adds 
more  area  to  the  flow  region,  thus  reducing  the  overall 
resistance  to  flow;  this  gives  rise  to  an  increase  in  the 
quantity  of  flow  for  the  same  pressure  gradient.  On  the 
other  hand,  the  presence  of  a  capillary  fringe  in  a  situation 
where  the  flow  is  mainly  in  the  vertical  direction,  causes 
an  increase  in  the  resistance  to  flow  tending  to  reduce  the 


total  flow. 
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Swertzendruber  and  Kirkham  (1956)  have  reported  from 
theoretical  considerations  that  the  contribution  by  the 
fringe  to  the  total  flow  may  reach  a  maximum  value  of  170 
per  cent  of  the  flow  beneath  the  phreatic  surface,  when  the 
length  of  the  flow  region  (parallel  to  the  overall  direction 
of  flow)  is  ten  times  the  thickness  of  the  flow  region 
beneath  the  phreatic  line.  The  fringe  component  of  flow 
rapidly  falls  off  with  increasing  thickness  of  the  flow 
region  beneath  the  phreatic  surface,  till  it  attains  a 
value  of  13  per  cent  for  a  length  -  thickness  ratio  of  two. 

Such  substantial  increases  of  flow  through  capillary  fringes 
may  have  important  implications  in  problems  involving  seepage 
through  earth  dams,  wherein  the  current  practice  is  to  ignore 
the  effects  of  capillarity. 

2.4.  Darcy*s  coefficient  of  permeability: 

The  various  methods  of  solution  of  the  Laplacian 
equation  (outlined  in  section  2.4.2)  require  that  the 
permeability  of  the  soil  be  known  in  advance.  There  are 
many  experimental  methods  available  for  the  determination  of 
this  quantity,  and  they  are  briefly  reviewed  elsewhere  in 
this  dissertation  (see  Appendix  B) . 

Having  realized  that  the  experimental  determination  of 
the  permeability  is  susceptible  to  many  errors  and  difficulties 
in  interpretation,  numerous  investigators  have  attempted  to 
correlate  the  permeability  to  other  physical  properties  of 
the  soil.  Such  attempts  have  never  proved  to  be  uniquely 
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successful  and  so  a  detailed  consideration  of  the  various 
theories  is  not  warranted.  A  brief  outline  of  these 
theories  is  given  in  Appendix  B;  for  greater  details,  the 
reader  is  referred  to  Scheidegger 1  s  monograph  (1957)  . 

2.5.  Flow  through  partially  saturated  soils. 

2.5.1.  Practical  cases  of  occurrence: 

Highway  and  airport  engineers  are  mainly  concerned 
with  flow  of  water  through  partially  saturated  soils  only, 
since  most  of  the  artificial  structures  like  roads, 
embankments  and  earth  dams  are  constructed  at  optimum 
moisture  content  of  the  soil  for  maximum  densif ication;  as 
is  well  known,  soils  compacted  under  such  conditions  have 
10  per  cent  or  more  of  their  voids  filled  with  air.  Movement 
of  water  through  such  soils  is  not  strictly  comparable  to 
that  through  saturated  soils. 

Considerable  amount  of  research  has  been  devoted 
to  the  task  of  evolving  a  theoretical  explanation  for  the 
phenomenon  of  flow  through  partially  saturated  soils,  and 
many  theories  have  been  put  forward  by  the  various 
investigators.  In  the  following  sections,  the  various 
theories  of  unsaturated  moisture  flow  are  briefly  discussed. 

2.6.  Theories  of  flow  through  partially  saturated  soils; 

2.6.1.  Nature  of  the  problem: 

In  comparison  to  the  case  of  flow  through  saturated 

media,  the  problem  of  flow  through  partially  saturated  media 
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is  beset  with  many  more  complicating  factors.  The  foremost 
among  these  is  the  dependence  of  the  permeability  on  the 
moisture  content  of  the  medium;  the  permeability  has  been 
found  to  depend  both  upon  the  potential  gradient  as  well 
as  upon  the  moisture  content  of  the  medium  (Klute,  1952 , 

Russel  and  Spangler  1941) .  In  contrast,  in  such  realms 
as  flow  of  electricity  or  of  heat,  the  rate  of  flow  depends 
solely  upon  the  potential  gradient  causing  flow.  Another 
complicating  factor  is  the  presence  of  numerous  three-phase 
interfaces  (soil-water-air)  within  the  region  of  flow, 
giving  rise  to  a  series  of  interfacial  electrical  and 
physico-chemical  phenomena,  such  as  differential  heat  of 
wetting  and  of  adsorption,  microplastic  swelling  and  so  on. 

And,  the  entire  process  is  not  thermodynamically  reversible; 
that  is,  the  rate  of  drying  of  a  wet  soil  is  different  from 
the  rate  of  wetting  of  a  dry  soil,  thus  introducing  hysteresis 
effects.  These  and  many  similar  difficulties  have  prompted 
many  investigators  to  try  various  hypotheses  to  explain  the 
phenomenon . 

These  various  hypotheses  may  be  broadly  classified 
into  three  groups:  In  the  first  group,  the  laws  of  viscosity 
and  surface  tension  have  been  used  in  conjunction  with  the 
expression  for  the  pressure  difference  across  a  curved 
meniscus  to  deduce  the  distribution  of  water  above  a  water 
table  (see,  for  example.  Keen,  1931,  Schofield,  1935) .  Such 
an  approach  gives  reliable  results  in  relatively  coarse-grained 


-  ...  i.  •.  .  .  . : 

...  ■.  i ...  '■  .:  1  'j:  '>  .'  ■  .. 

■  •  '■  ;  ...  .  ; 

■  •  -  ’ . ...  j.  .  ' 

x  x  }  :  .  '•  :  '  ;o 

,  .  .  r 

;/  ' . 

:  .  .  '  ''  . .  . : 

:  •  '■  :  j  .....  '■  :  ...  , 

x  ■  )  .  .  .  . 

.  .  .  ■  :  . 

X  ....  .  :  ... 

■  •  ■  -  , 

• ■  ■  v  , 

„  ■ 

'  ■  ' . f  , :  .  '  .  . 

. 

•  :  '•  •  •  •  ,  .'  '•  ... ..."  .;  u'1,1 

x  '■  :  :  :  '■  ' 

.  •  -  >  ■  '■  "  ,:'i  ■:  '  '.I.  .,  '•  ■:  ■  ...  ; 

.  '  •  '  '..CDS,  sons  1 

■  ■  '.  ">  /  '/  .  ...  .,  .  ;■  ,  1  .;  ;  ;;v 

■  •  '•  x  ■  x  ,  , 


soils  where  physico-chemical  reactions  are  negligible.  In 
the  second  group,  the  analysis  is  based  on  energy 
relationships,  and,  since  the  thermodynamical  energy 
relationships  must  hold  irrepective  of  the  mechanism  at 
work,  this  approach  is  fundamentally  sound  and  theoretically 
valid  for  all  types  of  soils.  In  the  third  group,  the 
moisture  flow  is  treated  as  a  diffusion  phenomenon  with  the 
rate  of  diffusion  depending  upon  the  moisture  concentration 
gradient.  As  is  usually  the  case,  all  these  hypotheses  are 
based  on  certain  simplifying  assumptions  -  both  implicit  and 
explicit  -  and  an  attempt  will  be  made  in  the  following 
sections  to  analyze  these  assumptions.' 

2.6.2.  Viscosity  and  surface  tension  theories: 

These  theories  are  all  based  on  the  validity  of 
Darcy ts  law  for  the  case  of  flow  through  partially  saturated 
soils.  Gardner  (1951)  was  presumably  the  first  to  suggest 
such  a  validity  and  now  it  is  generally  accepted  that 
Darcy *s  law  is  valid  "at  any  given  degree  of  saturation 
and  at  flow  rates  low  enough  for  laminar  flow  to  obtain  in 
unsaturated  soils"  (Russel  and  Klute,  1954) .  However,  in 
contrast  to  the  assumed  constancy  of  the  permeability  of 
any  given  saturated  medium,  the  permeability  of  a  partially 
saturated  medium  is  a  rapidly  varying  function  of  the 
moisture  content.  This  is  "due  to  the  fact  that  the 
effective  cross-section  available  for  flow  of  fluid  is 


. 

'  1. 

> 5  : 

:r  ''.at 

f J :  •  • 

..  .  ■  ■ 

;.7rC! 

.  ■ .  - 

..7 

7 

■  : 

(ta 

,  ^ 

.  •  ;> 

' 

:  : 

. 

•.  '  r  r  v 

3'  ' 

•• 

\ 

.  '  '• 

.  '■ 

• 

'  :  l 

:  :  :: 

n;  :: . 

' 

. 

/.  ,  '■  .  .  .■  .  o  , : 

'?  : :  '  '•  , .  ■  •*  '  '•  .■  '  ,  n, '  • j  •.  - '■> 

l  O.'-  «• i ;  Dj:;h >a 

... 

.  :  ''  ;  ■  <  "J 

r  .  :  7  '  1 

••  j  '  '.i.  ■;  .  i  '  .  . 

■  v  '■  ’  '■/'  :7  :  ;  ■  -v.'.  .  ■  •  /  .  '  j. 

1 

■■  ;;  T  •;  ;h~  . 

.  '  ,  .  ■  '  ::  >  ■' 

'■  '  '■  i-.: : 

11 '  -:r: :  ■. :  ; .  ;  v; ;  .  v.;.:r';  . 

■  .• 

..  ‘  ;  /  '  i  •  •  ■; ;j::>  s\ ;; : 


reduced  in  an  unsaturated  medium  as  compared  to  that  of  a 

saturated  medium"  „  .  .  (and  the  conductivity  in  an 

unsaturated  medium)  will  decrease  with  decreasing  moisture 

content  of  the  medium"  (Klute,  1952) . 

Having  established  the  validity  of  Darcy *s  law  to 

the  present  problem,  the  authors  of  these  theories  combine 

the  law  with  the  fundamental  law  of  conservation  of  matter 

to  obtain  a  final  equation  of  the  form  (Gardner,  1958)  : 

v  <£>] 

where  e  is  the  water  content  on  a  volumetric  basis 

k  (e)  is  the  so-called  capillary  conductivity,  assumed 
to  be  a  function  of  (©)  and 

cj>  is  the  potential  function  per  unit  of  mass  responsible 
for  flow. 

The  quantity,  k  (e)  has  been  variously  described  as 
the  capillary  conductivity,  coefficient  of  aqueous 
conductivity,  unsaturated  permeability,  permeability, 
conductivity  and  so  on  (see  Gardner,  1957,  Klute,  1952, 

Russel  and  Spangler  1941,  Klute  et  al .  1956)  .  The  writer 
favours  the  use  of  the  term  "unsaturated  permeability"  since 
it  is  brief  and  descriptive. 

As  already  mentioned,  the  unsaturated  permeability 
is  a  variable  function  of  the  (volumetric)  moisture  content 
of  the  medium  and  possibly  of  other  factors  as  well. 

Although  Gardner  and  Gardner  (1951)  express  the  opinion  to 
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this  effect,  Gardner  (1956)  assumes  that  the  unsaturated 
permeability  is  a  single-valued  function  of  the  moisture 
content  alone . 

Regarding  the  potential  function,  (f),  many 
investigators  (including  Gardner)  have  assumed  that  the 
total  potential,  0,  is  the  sum  of  capillary  potential,  4)  , 
and  the  potential  of  the  gravitational  force  field,  gz, 
all  potentials  being  expressed  for  one  unit  mass  of  the 
substance.  Thus 


where  z  is  the  vertical  distance  to  the  location  of  the 
unit  mass  from  the  datum,  and 
g  is  the  acceleration  due  to  gravity. 

Expressing  all  potentials  in  terms  of  the  gravitational 
constant,  the  capillary  potential  becomes  the  "suction" 

head  of  soil  water,  S,  and  eg.  (8)  becomes 

_Lvcb  -  ^7  (  -  S  +  *) 

9 


Substituting  equation  (9)  in  the  basic  equation,  Gardner 
arrives  at  the  following  equation  of  flow: 

^  =  _  -7.  (10) 

At  this  stage  of  the  derivation,  he  introduces  the 
further  assumption  that  the  moisture  content  and  the  suction 
head  are  related  by  a  single-valued  function.  Experiments 
conducted  by  numerous  investigators  (Schofield  1935,  Klute 
1952,  Croney  et  al„  1948,  1958,  Klute  et  al.  1956,  to  name 
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but  a  few)  have  established  that  this  is  not  so  and  that 


there  are  hysteresis  effects  to  be  observed  between 
wetting  and  drying  processes.  Thus,  the  assumption  of  a 
reversible  relationship  between  moisture  content  and 
suction  cannot  be  justified.  However,  the  solution  of  the 
non-linear  differential  equation  (10)  will  be  extremely 
difficult,  if  not  altogether  impossible,  without  this 
simplification.  With  this  assumption,  k  (9)  becomes  k  (S) , 
that  is  the  unsaturated  permeability  becomes  a  function  of 
the  suction  head  of  the  soil  water. 


The  nature  of  the  functions  k  (0)  and  k  (S)  is  not 


well  known  (Gardner  1956) .  However,  for  uni -directional 
steady-state  flow  conditions,  equation  (10)  reduces  to 


(ID 


where  q  is  the  volume  flux. 

Solutions  of  this  equation  have  been  reported  for  particular 
forms  of  the  function  k  (©)  =  k  (S)  .  As  reported  by 

Gardner  (1958) ,  Richards  assumed  that  k  (S)  =  a  S  +  b.  Wind 


-3  /2 

assumed  k  (S)  =  aS  ,  Remson  and  Fox  assumed  that  k  (S) 
-1 

aS  ,  while  Gardner  himself  set 


k  (S)  “  (Sh  +  b) 


(12) 


and  gave  solutions  for  various  values  of  the  exponent,  n. 

In  each  case  where  Gardner  tried  to  fit  the  solution  of  the 


flow  equation  to  the  experimental  data  furnished  by  previous 
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investigators,  he  was  able  to  obtain  a  good  agreement  after 
manipulating  the  constants  a,  b  and  n.  In  the  writer's 
opinion,  this  cannot  be  entertained  as  conclusive  proof  of 
the  soundness  of  the  theoretical  analysis  unless  and  until 
one  is  able  to  obtain  the  values  of  the  k  (S)  function 
beforehand  and  check  the  predicted  relationship  against  the 
experimentally  observed  data.  To  the  writer's  knowledge, 
no  such  rigorous  verification  has  been  reported  so  far. 

Working  along  similar  lines,  but  recognizing  that 
the  moisture  content  -  soil  suction  and  unsaturated 
permeability  -  soil  suction  relationships  are  not  perfectly 
reversible.  Miller  and  Miller  (1955)  have  developed  a 
'physical  and  mathematical  theory  of  capillary  flow1.  They 
argue  that  'the  liquid  filled  spaces  -  and  therefore  the 
water  content  and  the  capillary  conductivity  as  well  -  depend 
upon  the  tension  and  upon  the  previous  history  thereof, 
but  not  upon  the  time  rate  at  which  successive  portions  of 
that  history  occurred,  the  order  of  events  being  preserved1. 
They  report  that  the  final  differential  equations  obtained 
by  eliminating  the  flow  velocity  from  the  generalized  form 
of  Darcy's  law  and  the  pertinent  law  of  conservation  of 
matter,  is  identical  with  the  equation  derived  by  Richards 
in  1931 As  an  improvement  over  the  earlier  formulation, 
the  present  authors  have  appended  a  complete  list  of  "reduced 


1 


The  writer  has  not  been  able  to  obtain  a  copy  of  the 
paper  published  by  Richards. 
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variables"  from  considerations  of  dynamical  similitude, 
thus  enabling  the  comparison  between  two  geometrically 
similar  media.  This  involves  that  fictitious  quantity 
called  the  "characteristic  pore  size"  and  the  exact 
definition  of  this  characteristic  size  has  not  been  given. 

Klute  and  Wilkinson  (1958)  have  reported 
experimental  data  supporting  this  "similar  media  concept" 
of  Miller  and  Miller.  They  determined  the  unsaturated 
permeability  as  a  function  of  the  moisture  content  for 
five  different  size  fractions  of  lplainfield  sand1  and 
plotted  the  data  in  terms  of  the  reduced  variables  suggested 
by  Miller  and  Miller.  When  so  plotted,  the  data  seemed  to 
coalesce,  in  accordance  with  the  theory,  into  a  single 
curve  "within  limits  of  experimental  error." 

2.6.3.  Diffusion  theories: 

The  application  of  the  potential  theory  -  the 
surface  tension  and  viscous  resistance  theories  discussed 
in  an  earlier  section  -  has  been  restricted  to  the  case  of 
steady-state  flow  conditions  only,  since  the  differential 
equation  describing  non-steady  flow  conditions  is  non-linear 
with  attendent  difficulties  of  mathematical  solution.  In 
fine-grained  soils,  steady  flow  conditions  are  seldom 
realized  due  to  the  extremely  low  permeability  of  the  soils. 
For  dealing  with  such  problems,  some  other  approach  is 
clearly  needed  (Childs  1936) .  One  possible  method  of  approach 
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is  to  view  the  phenomenon  as  one  of  diffusion  of  moisture 
initiated  by  the  moisture  concentration  gradient.  It  may 
be  noted  that  the  flow  equation  (Eq.  (10))  reduces  to  the 
non-linear  diffusion  equation,  if  the  gravity  term  (  ^  ) 
is  omitted. 


2. 6. 3.1.  Childs1  theory  assuming  constant  diffusivity: 

Childs  (1936)  assumed  that  the  mass  of  water 
passing  a  unit  cross-section  of  a  soil  mass  in  unit  time 
is  directly  proportional  to  the  moisture  concentration 
gradient,  expressed  by  an  equation  of  the  form: 


M  M 

at  Sx 


(13) 


5Q 


where,  tfr  =  the  mass  rate  of  flow  of  water  in  the  x  - 


direction, 

D  =  the  diffusion  coefficient,  and, 

©  =  the  moisture  concentration,  "expressed  as  a 

fraction  of  the  dry  matter"  (the  ratio  of  the 
mass  of  water  to  the  mass  of  dry  matter  in  a 
given  volume)  . 

Combining  this  equation  with  continuity  equation,  and 
eliminating  the  mass  flux,  Childs  arrived  at  the  final 
equation: 

U4) 


dt 


d  be 

G  dx^ 


where,  s’  equals  the  dry  unit  weight  of  the  soil.  Assuming 
that  the  ratio  of  the  diffusion  coefficient  to  the  unit 


weight  of  the  soil  is  independent  of  depth  and  time,  Childs 
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obtained  a  solution  of  the  above  equation  in  the  form  of  a 
Fourier  series.  By  the  very  nature  of  the  solution,  the 
moisture  profile  is  obtained  as  a  portion  of  an  infinite 
periodic  function,  the  exact  shape  of  which  is  dependent 
on  the  initial  and  boundary  conditions  of  the  particular 
problem . 

The  nature  and  soundness  of  the  assumptions  made 
in  the  above  formulation  are  worthy  of  discussion,  if  only 
to  point  out  the  inherent  weaknesses  of  the  arguments 
advanced.  In  the  first  place,  it  is  implicitly  assumed 
that  "the  equilibrium  profile,  which  the  soil  moisture 
would  exhibit  if  given  sufficiently  long  time  of  undisturbed 
rainless  conditions  .  .  .  would  indicate  constant  moisture 

content  irrespective  of  depth"  (Nicholson  and  Childs,  1936) . 
Secondly,  the  ratio  of  the  diffusion  coefficient  to  the  dry 
unit  weight  of  the  soil  is  assumed  to  be  independent  of 
depth  and  time.  The  admissibility  of  the  former  assumption 
cannot  be  verified  by  direct  experiments,  since  the  necessary 
conditions  can  seldom  be  fulfilled.  The  latter  assumption  has 
been  proved  to  be  unwarranted  (Klute  1952,  Philip  1958) :  in 
the  course  of  many  infiltration  experiments,  it  has  been 
observed  that  there  is  a  distinct  wetting  front,  across  which 
there  is  a  rather  sharp  moisture  gradient  (Bodman  and  Coleman 
1944,  Coleman  and  Bodman  1944,  and  others) .  The  simple 
theory  assuming  a  constant  diffusivity  fails  to  predict  any 
such  wetting  front.  Klute  (loc.  cit.)  has  shown  that  a 
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wetting  front  is  indeed  predicted,  if  the  diffusivity  is 
treated  as  a  variable  function  of  the  moisture  concentration. 
In  a  more  recent  paper,  Childs  and  George  (1948)  state  a 
form  of  Darcy's  law,  wherein  the  dependence  of  the 
diffusivity  on  the  moisture  content  is  implied. 

2. 6. 3. 2.  Diffusion  theory  with  concentration-dependent 


diffusivity : 


In  an  attempt  to  develop  a  numerical  method  of 


solution  for  the  flow  equation  (Eq.  10 )  ,  Klute  was  able  to 
derive  the  diffusion  type  of  equation  explicitly  from  the 
continuity  equation  and  the  Darcy *s  law  (Klute  1951) „  The 
final  equation  so  obtained  is  of  the  form: 


wherein,  in  addition  to  the  symbols  already  defined, 

D (p)  is  the  diffusivity  function,  assumed  to  be  a  function 
of  the  pressure,  p,  in  the  soil  water, 
k(p)  is  the  unsaturated  permeability  of  the  soil,  again 
assumed  to  be  a  function  of  the  pressure, 

P  is  the  density  of  water,  and 

p.  is  the  viscosity  of  the  water  (see  Klute  et  al.  1956)  . 
It  will  be  noticed  that  if  the  last  term  on  the  right  side 
of  the  equation  is  dropped,  and  if  the  diffusivity  is 
treated  as  a  constant,  Eq.  (15)  reduces  to  the  same  form  as 


postulated  by  Childs. 
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As  equation  (15)  is  again  non-linear,  analytical 
solutions  are  very  difficult;  in  order  to  overcome  this 
difficulty,  Klute  (1951)  and  Philip  (1955)  have  developed 
numerical  methods  of  solving  such  equations  for  horizontal 
flow  systems  of  semi-infinite  extent.  Although  the 
methods  are  very  tedious,  they  do  reveal  the  existence  of 
a  locally  sharp  moisture  gradient,  that  is,  a  wetting  front, 
in  problems  involving  infiltration  of  water  into  an 
initially  dry  soil.  Gardner  (1959)  has  developed  a 
numerical  method  of  solving  the  diffusion  type  flow  equation 
for  the  case  of  isothermal  evaporation  of  water  from  soils; 
he  has  also  indicated  a  transformation  of  the  variables  which 
enables  the  inclusion  of  the  gravity  term  into  the  analysis. 
On  the  strength  of  the  results  of  drying  experiments  of 
artificially  prepared  soil  columns,  Gardner  states  that 
the  experimental  cumulative  evaporation  from  soil  columns 
agreed  well  with  the  theoretical  analysis. 

2.6.40  Theories  based  on  energy  relationships: 

2. 6 .4.1.  Preliminary  consideration: 

With  the  introduction  of  the  term  capillary  potential 
into  soil  moisture  studies  by  Buckingham  (1907)  and  with  the 
elaboration  of  the  scope  of  this  quantity  by  Gardner  and 
others,  energy  considerations  in  unsaturated  flow  problems 
have  received  ever-increasing  attention.  It  is  a  matter  of 
common  observation  that  free,  pure  water  will  spontaneously 
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enter  a  partially  saturated  soil;  defining  the  datum  of 
energy  as  that  of  the  free  water,  it  is  at  once  apparent 
that  the  energy  of  the  soil  water  in  a  partially  saturated 
soil  is  negative  in  character.  Buckingham  defined  the 
capillary  potential  as  the  difference  in  free  energy 
between  soil  water  and  free  water,  expressed  as  the 
equivalent  height  of  a  column  of  water  in  centimetres. 
Schofield  (1935)  introduced  the  logarithmic  scale  and  set 

PF  =  1O^10H  ('S) 

where  H  is  the  height  in  centimetres  of  the  equivalent 
water  column.  The  symbol  in  analogy  with  the  symbol 

used  in  Sorensen's  acidity  scale,  signifies  the  logarithm 
of  the  height  of  the  water  column  to  give  the  amount  of 
negative  “suction®  head  in  question. 

Schofield  gives  the  following  equations  for 
determining  the  1?F  of  soil  water  as  a  function  of  the  relative 
humidity  of  the  soil  atmosphere  or  the  depression  of  the 
freezing  point  of  water  in  soils: 

£»  =  6.5  +  log.  _  (2  -  log  h)  for  water  at  20°C  (\7) 

and  =4.1  +  logn_t  for  water  at  freezing  point.  (18) 

In  these  equations,  h  is  the  relative  humidity  of  the  soil 
atmosphere  in  per  cent  and  t  is  the  depression  of  the  freezing 
point  of  soil  water  in  degrees  Centigrade.  Croney  et  al . 

(1948)  and  Dhawan  et  al.  (1956)  have  shown  that  the  first  of 
these  equations  can  be  readily  derived  as  a  consequence  of 
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Kelvin tt s  fundamental  equation  describing  the  relationship 
between  an  adsorbate  and  any  porous  adsorbent.  The  second 
equation,  relating  the  change  in  P^,  with  the  freezing  point 
depression  is  again  a  simple  consequence  of  the  classical 
laws  of  thermodynamics  (see  Appendix  C;  also  Edlefson  and 
Anderson  1943) .  Thus,  the  energy  relationships  invoked 
in  this  theory  may  be  regarded  as  established  on  sound 
theoretical  grounds  of  proven  validity. 

2. 6. 4. 2.  Principle  of  practical  application: 

Having  seen  the  soundness  of  the  underlying 

framework  of  the  theory,  let  us  see  how  it  can  be  applied 

to  the  problem  of  movement  of  water  through  soils.  It  is, 

of  course,  known  that  two  bodies  at  different  energy  levels 

cannot  remain  in  equilibrium.  In  the  same  way,  if  two 

sections  of  a  soil  mass  are  at  different  energy  levels  (or 

Pp  levels) ,  a  state  of  equilibrium  can  never  exist; 

consequently,  some  movement  must  take  place  in  order  to 

equalize  the  energies  of  the  co-existing  phases.  Thus, 

soil  water  will  always  migrate  from  regions  of  low  P^  to 

.  * 

regions  of  higher  Pp  until  the  two  values  are  equal,  due 
consideration  being  given  to  the  potential  energy  of  the 
various  sections  of  the  mass  due  to  their  positions  in  the 
gravitational  force  field.  It  is  interesting  to  note  that 
water  may  actually  migrate  from  regions  of  low  moisture 


A  low  value  of  P signifies  a  low  numerical  value  of  the 
free  energy;  but,  as  the  free  energy  of  soil  water  is 
negative,  the  absolute  magnitude  of  free  energy  in  this 
case  is  greater  than  if  the  Pp  were  higher. 
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content  to  regions  of  higher  moisture  content  (the  moisture 
contents  being  expressed  on  a  weight  basis) ;  thus,  the 
migration  of  water  in  soils  is  governed  by  the  gradient 
and  not  the  moisture  content  gradient.  This  is  true 
irrespective  of  the  particular  mechanism  at  work,  that  is, 
whether  the  migration  is  in  the  vapour  phase  under  a  vapour 
pressure  gradient  or  in  the  liquid  phase  under  a  pressure 
gradient . 

If  the  at  any  point  in  a  soil  mass,  initially  in 

r 

a  state  of  moisture  equilibrium  is  changed,  the  soil  water 
will  redistribute  itself  until  a  new  state  of  equilibrium 
of  energies  is  reached.  However,  under  field  conditions,  a 
state  of  equilibrium  could  hardly  be  set  up,  since  the 
initial  and  boundary  conditions  are  perpetually  varying, 
such  variations  being  caused  by  infiltration  of  rain  water, 
evaporation  into  atmosphere,  and  similar  factors.  But,  if 
the  soil  mass  is  covered  by  an  impervious  blanket  -  as  is 
done  in  the  case  of  highway  and  airport  pavements,  an 
equilibrium  distribution  of  moisture  content  beneath  the 
pavement  is  obtainable,  and  such  an  equilibrium  distribution 
can  be  approximately  determined  from  the  energy  relationship 
so  far  seen. 

Such  determinations  have  been  attempted  in  the  Road 
Research  Laboratory  England  (Croney  and  Coleman  1948, 

Croney  et  al.  1950,  Croney  1952,  Croney  et  al.  1958)  with 
success.  Two  assumptions  have  been  made  in  these  attempts; 
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(1)  that  the  surfacing  is  impervious  to  water  in  the  liquid 
and  vapour  forms;  and,  (2)  that  the  moisture  transfer  between 
the  soil  and  the  atmosphere,  which  occurs  at  the  boundary  of 
the  covered  area,  has  a  negligible  effect  on  the  moisture 
distribution  beneath  the  pavement.  The  first  of  these 
assumptions  is  readily  fulfilled  by  most  of  the  flexible 
and  rigid  pavements  normally  used  in  civil  engineering 
practice;  the  second  assumptions  is  not  so  easily  satisfied 
in  many  cases,  and  thus  the  analysis  can  be  considered  only 
as  a  first  approximation. 

In  incompressible  soils,  the  overburden  pressure  is 
entirely  carried  by  the  soil  grains  and  the  soil  beneath 
the  ground  water  table  is  essentially  saturated  with  water. 

At  any  distance,  x,  above  the  water  table,  the  1?  is  given 
by  log^Q  x,  and  the  corresponding  equilibrium  moisture 
content  can  be  readily  obtained  from  experimentally  determined 
sorption  isotherms.  Due  to  the  presence  of  hysteresis 
effects,  two  equilibrium  profiles  will  be  obtained,  depending 
on  the  initial  moisture  conditions  of  the  soil.  If  the 
moisture  profile  at  the  time  of  paving  the  surface  is  known, 
a  single  equilibrium  profile  can  be  easily  approximated. 

Such  approximations  have  proved  to  be  successful  in  many 
cases  in  England  for  comparatively  shallow  water  table 
conditions  (depth  not  more  than  about  10  feet)  (Coleman  1959) . 

The  case  of  compressible  soils,  like  heavy  clays,  is 
complicated  due  to  the  possibility  of  presence  of  pore  water 
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pressures  from  other  causes.  Also,  a  part  of  the  overburden 
pressure  may  be  carried  by  the  pore  water,  and  thus,  the 
suction  at  any  level  is  dependent  on  the  position  of  the 
water  table  as  well  as  the  overburden  pressure,  which 
depends  on  the  bulk  density  of  the  soil  and  hence  on  the 
moisture  profile.  If  a  fraction,  a,  of  the  overburden 
pressure,  P,  is  carried  by  the  pore  water,  then  the 
following  equation  can  be  written: 

a.p  +  s  =  u  (;i9); 

where,  IS5  is  the  suction  in  pressure  units  and  *u5  is  the 
pore  water  pressure  determined  by  the  position  of  the  water 
table.  (Note  that  both  u  and  S  are  negative  for  all  planes 
above  the  water  table) .  A  successive  approximation  technique 
may  be  used  for  obtaining  the  equilibrium  suction  values  at 
the  different  levels:  assuming  a  tentative  moisture  profile, 
the  bulk  densities  at  various  elevations  may  be  computed 
and  substituted  in  Eq.  (19)  to  obtain  the  first  set  of 
suction  values;  these  values  may  then  be  used  in  conjunction 
with  the  corresponding  sorption  isotherms  to  obtain  a  closer 
approximation  to  the  moisture  profile  and  the  whole  procedure 
may  be  repeated  till  the  desired  degree  of  refinement  is 
obtained . 

There  seems  to  be  no  reliable  method  for  obtaining 
the  value  of  the  co-efficient  'a1  in  Eq.  (19) ;  in  compressible 
soils,  volume  changes  occur  with  changes  in  moisture  content 
and  thus,  the  shrinkage  characteristics  of  the  soil  may  be 
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utilizable  in  determining  the  value  of  'a1.  Work  along  these 
lines  is  reportedly  under  way  in  the  Road  Research  Laboratory 
England  (Croney  1952) .  Further,  efforts  to  correlate  the 
suction  properties  with  the  index  properties  of  the  soil 
have  also  been  reported  „  Although  the  application  of  the 
concept  to  the  problem  of  migration  of  water  in  soils  has 
had  a  respectable  amount  of  success,  much  work  remains  to 
be  done?  the  investigations  so  far  have  been  concerned  with 
shallow  water  table  conditions  only;  with  the  water  table 
at  greater  depths,  continuous  air  voids  are  distinctly 
possible,  thus  promoting  active  transfer  of  water  in  the 
vapour  phase.  To  the  writer^s  knowledge,  the  E>F  concept 
has  not  been  applied  to  movement  of  water  in  the  vapour  phase 
Thus,  in  order  to  predict  the  direction  and  magnitude 
of  movement  of  water  under  a  given  set  of  conditions,  one 
must  know  the  relationship  between  moisture  content  and  F 

F 

over  the  full  range  of  moisture  contents  for  the  soil. 

Knowing  this  relationship  and  the  initial  and  boundary 
conditions  of  the  particular  problem,  it  is  possible  to 
determine  the  equilibrium  distribution  of  moisture  within  the 
mass,  that  is,  the  distribution  that  will  be  attained,  if 
the  conditions  remained  constant  for  a  sufficient  period. 
Hence,  a  considerable  amount  of  research  has  been  devoted  to 
the  task  of  perfecting  the  techniques  for  obtaining  the 
so-called  sorption  isotherms — the  moisture  content-suction 


curves — of  soils. 


L  T  \ 

.  * 

1 

'  : 

' 

• 

° 

'  . 

:  . 

.  . 

’ 

~  fii) : 

V  : 

;  .v:  •  -V... 

:  c-  . 

.  : 

/ 

. 

V 

V 

.  :  a 

'• 

r 

'  -  i  '  j* 

V 

1 

1 

■  ■ 

'  f< 

:/  1 

a  ’  : 

' 

' 

\  ' 

: 

7 

* 

rx.  ' ' 

:  - .  :  : 

j 

'  •  : :  '•  r  \ 

■ 

- 

Oil  : :  1 

...  r.- 

[ 

•  : 

'■  V 

:  ~J.  '  ; 

'  i  .  ■' ... 

• 

"  ■  O'..* 

’ 

> 

: ;  . 

J'.  A  J  .  . 

< 

. 

.  '  . 

::r; 

:: 

- :  :  i  ■ 

. 

: :  s  £■ 

. 

V 

■ 

: X  ■ 

. :  . : 

] 

. 

*•  •  ••  , 

'  1 

:  . : '  '•  i  j 

:  /  •:  :  .  . 

'■ 

i:  r.: 

A 

.  .  . 

2. 6. 4. 3.  Methods  of  obtaining  sorption  isotherms: 

The  entire  range  of  moisture  content  of  a  soil — 

from  saturation  to  oven  dryness — is  covered  by  a  Pp  range 

of  0  to  approximately  7.  No  single  method  can  possibly  be 

relied  upon  to  give  the  moisture  content  -  P  relationship 

F 

of  a  soil  over  such  a  vast  range  of  values.  Fortunately, 
there  are  a  number  of  methods  available,  each  one  being 
suitable  for  a  particular  range  of  suction  values.  These 
various  methods  are  listed  below  (after  Croney  et  al „  1952) : 
Direct  methods- 
Suction  Plate 
Tensiometer 
Centrifuge 
Pressure  Membrane 
Indirect  Methods- 

Freezing  point  depression 
Vacuum  desiccator 
Sorption  balance 
Electrical  resistance  gauges 
The  suction  plate  and  the  tensiometer  methods 
employ  similar  principles:  the  soil  is  allowed  to  come  into 
equilibrium  with  a  large  quantity  of  free  water  contained  in 
a  free  water-manometer  system;  as  some  of  the  free  water  tends 


Pp  range  0  to  3 
Pp  range  0  to  3 
Pp  range  3  to  4.5 
Pp  range  0  to  4 


Pp  range  3  to  4 
Pp  range  4.5  to  7 
Pp  range  4.5  to  7 
Pp  range  3  to  7 


to  enter  the  soil  to  bring  about  a  state  of  equilibrium,  the 
pressure  in  the  free  water-manometer  system  decreases  by  an 
amount  equal  to  the  pressure  deficiency  in  the  soil  or  the 
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soil  suction,  which  is  indicated  by  the  manometer.  The  range 
of  utility  of  these  two  instruments  is  restricted  by  the  fact 
that  they  operate  in  a  gravitational  field  and  can  never 
exceed  one  atmosphere  of  pressure  deficiency  or  suction. 

The  centrifuge  method  involves  supplementing  the  gravitational 
force  field  with  a  centrifugal  force  field  and  thus  the 
range  of  utility  is  considerably  increased.  In  the  pressure 
membrane  apparatus,  one  face  of  the  specimen  inside  a 
pressure  chamber  is  maintained  at  atmospheric  pressure  in 
equilibrium  with  a  body  of  free  water  through  permeable 
membrane;  the  chamber  pressure  is  raised  to  some  definite 
value  above  that  of  the  atmosphere  and  the  soil  specimen  is 
allowed  to  reach  equilibrium  with  the  external  body  of  water. 
If  the  pressure  in  the  chamber  is  now  suddenly  reduced  to 
atmospheric,  and,  if  the  soil  specimen  is  removed  from  the 
chamber  before  any  reimbibition  of  water  can  occur,  the 
moisture  content  of  the  specimen  will  correspond  to  a 
particular  value  of  suction  which  is  numerically  equal  to 
the  chamber  pressure  applied  initially.  A  number  of  trials 
may  be  made  with  various  chamber  pressures  in  order  to  define 
the  sorption  isotherm;  recent  improvements  of  the  apparatus 
enable  one  to  follow  the  changes  in  moisture  content  of  the 
specimen  without  having  to  remove  the  specimen  from  the 
pressure  chamber. 

Of  the  indirect  methods,  the  method  relying  on  the 
depression  of  the  freezing  point  of  soil  water  is  based  on  the 
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fact  that  the  freezing  point  of  water  decreases  with  decreasing 
free  energy  (see  Appendix  C;  also  Edlefson  and  Anderson  1943) . 
The  reason  for  this  depression  is  not  fully  known,  but  that 
does  not  invalidate  the  thermodynamical  arguments  that  are 
used  for  deriving  the  expression  relating  the  freezing  point 
depression  to  the  free  energy  or  the  ]?F  of  soil  water  (Croney 
1952) .  As  freezing  is  essentially  a  drying  process,  the 
results  obtained  by  this  method  define  the  desorption  curve 
only.  The  vacuum  desiccator  method  is  based  on  the  principle 
that  the  free  energies  of  several  co-existing  phases  in 
equilibrium  are  equal  to  one  another.  In  the  test,  the  soil 
specimen  is  placed  in  a  desiccator  of  known  humidity  (obtained 
by  using  sulphuric  acid  of  suitable  concentration)  and  the 
moisture  content  of  the  specimen  is  determined  after 
equilibrium  conditions  are  reached.  From  the  known  relative 
humidity  of  the  atmosphere  inside  the  desiccator — which  is 
the  same  as  that  of  the  soil  air — the  free  energy  of  the 
soil  air  can  be  computed,  which,  in  its  turn,  is  equal  to  that 
of  the  soil  water.  The  time  required  for  reaching  equilibrium 
is  radically  reduced  by  using  a  vacuum  inside  the  desiccator; 
in  this  way,  the  time  will  be  about  three  days  for  a  2-gm. 
specimen,  which  may  require  about  two  months  under  atmospheric 
pressure  conditions  (Croney  1952) . 

The  principle  of  the  sorption  balance  method  is  the 
same  as  that  of  the  vacuum  desiccator  method,  in  that  the 
relation  between  the  free  energy  of  a  liquid  and  the  pressure 
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in  the  vapour  phase  in  equilibrium  with  it,  is  made  use  of. 
While  the  humidity  is  controlled  and  measured  by  the 
concentration  of  the  acid  in  the  latter  case,  the  required 
humidity  is  directly  obtained  in  the  former  case  by  allowing 
water  to  evaporate  into  a  previously  evacuated  container; 
the  value  of  the  relative  humidity  is  deduced  from  the 
measured  vapour  pressure.  A  known  weight  of  the  soil —  about 
2  gms0 — is  placed  on  the  pan  of  a  special  type  of  balance 
inside  a  chamber;  enough  water  vapour  is  admitted  into  this 
chamber  to  give  the  desired  suction  value  (as  computed  from 
the  vapour  pressure  data) ,  and  the  changes  in  moisture 
content  of  the  specimen  due  to  changes  in  the  vapour  pressure 
are  continuously  recorded  by  the  special  balance .  Absorption 
and  desorption  curves  can  be  readily  obtained  by  increasing 
or  decreasing  the  vapour  pressure. 

And,  finally,  the  electrical  resistance  method  of 
determining  the  E»F  aims  at  deducing  the  free  energy  of  the 
soil  water  from  the  noted  change  in  the  free  energy  of  a 
porous  plug  with  which  it  is  in  equilibrium;  the  porous  plug 
is  provided  with  two  electrodes  embedded  in  it,  the  resistance 
between  which  increases  with  increasing  suction.  Due  to 
difficulties  in  calibration,  the  suction-moisture  content 
relationship  of  the  porous  plug  and  the  moisture  content- 
electrical  resistance  relationships  are  separately  determined 
from  which  the  suction-electrical  resistance  relationship  is 
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obtained  by  interpolation.  This  method  has  not  been  found  to 
be  useful  for  laboratory  work  (Croney  et  al.  1951) . 

2.6„404.  Hysteresis  effects: 

All  experimental  data  assembled  so  far  have  one 
unique  feature  in  common:  the  presence  of  hysteresis  between 
the  absorption  and  desorption  isotherms;  for  the  same  value 
of  moisture  content,  the  corresponding  to  the  desorption 

r 

(drying)  process  is  higher  than  the  corresponding  to  the 
absorption  (wetting)  process.  Thus,  although  Buckingham 
(1907)  assumed  that  a  soil  exerts  a  definite  equilibrium 
suction  depending  solely  on  its  moisture  content,  this 
“equilibrium  suction1  is  found  to  depend  very  much  on 
whether  the  soil  is  being  wetted  or  dried  (Schofield  1935) . 
The  reason  for  this  marked  hysteresis  effect  is  generally 
believed  to  be  that  the  cross-sectional  areas  of  individual 
pores  of  the  soil  vary  widely  from  point  to  point.  Miller 
and  Miller  state  "A  given  pore  will  empty  suddenly  when 
the  tension  is  just  sufficient  to  “suck  in'  the  air-water 
interface  through  its  largest  entry  port.  To  refill  the 
pore,  the  tension  must  be  reduced  until  the  surface  tension 
of  the  soil  water  can  pull  the  interface  past  the  largest 
internal  cross-section  of  the  pore,  at  which  point  the 
refilling  process  will  again  occur  suddenly.  ...With  finer 
media,  the  individual  jerks  (in  emptying  and  filling  of  pores) 
will  not  be  observed  macroscopically ,  but  their  smoothed-out 
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effect  will  still  be  present  as  hysteresis."  Whatever  be 
the  cause „  the  presence  of  hysteresis  is  an  established 
fact,  and  this  constitutes  one  of  the  grounds  for  questioning 
the  validity  of  the  approach  of  the  potential  theory  (see 
section  2.72.  above)  which  tacitly  assumes  a  single-valued 
relationship  between  moisture  content  and  the  capillary 
potential . 

2.6.5.  Comparison  of  the  various  theories: 

The  potential  theory — based  on  surface  tension, 
viscous  resistance  and  the  hydrodynamic  aspects  of  the 
problem — can  be  made  sufficiently  general  in  its  formulation 
to  permit  solutions  of  most  problems  of  practical  interest 
(Muskat  1937) .  However,  solution  of  problems  involving 
transient  flow  conditions  is  very  difficult,  on  account 
of  the  non-linearity  of  the  differential  equation.  In  such 
a  case,  it  is  helpful  to  treat  the  problem  as  one  involving 
diffusion  of  water  governed  by  a  concentration-dependent 
diffusivity,  and  resort  to  numerical  methods  of  solution. 

Both  these  theories — potential  theory  and  diffusion 
theory — are  based  on  an  indispensable  assumption  of  a  single¬ 
valued  relationship  between  moisture  content  and  soil  suction 
or  the  unsaturated  permeability  of  the  soil.  The  "hysteresis 
functions"  concept  postulated  by  Miller  and  Miller  (1955) 
promises  to  be  a  definite  improvement  over  this  questionable 
assumption,  but  no  experimental  method  of  evaluating  these 
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"hysteresis  functions"  has  so  far  been  developed.  Philip 
(1958)  claims  that  the  presence  of  hysteresis  need  not 
necessarily  invalidate  the  theoretical  arguments  and  that 
the  effect  of  hysteresis  over  small  suction  ranges  normally 
encountered  in  practice „  is  only  small. 

The  1?  concept  is  theoretically  applicable  to  all 

r 

types  of  soil  moisture  problems.  The  reports  of  practical 
applications  of  the  concept  to  deal  with  equilibrium 
distribution  of  water  beneath  impervious  pavements  in 
England  are  very  encouraging.  There  seems  to  be  a  vast 
scope  for  further  research  in  this  field;  in  view  of  the 
fact  that  the  thermodynamic  relationships  are  universally 
true,  it  should  be  possible  to  evolve  a  generalized 
theoretical  analysis  to  solve  problems  involving  polyphase 
transfer  of  water  through  soils,  besides  providing  adequate 
explanations  for  such  controversial  phenomena  as  cohesion 
and  swelling  of  clays. 
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CHAPTER  III 


FLOW  OF  WATER  UNDER  ELECTRIC  POTENTIALS 


3.1.  Introduction : 

The  subject  matter  of  this  chapter  is  the  complex 
phenomenon  of  electro-osmosis  of  water  through  fine-grained 
soils.  Although  electro-osmosis  is  not  generally  directly 
responsible  for  the  accumulation  of  water  beneath 
impervious  pavements,  a  brief  outline  of  the  basic 
phenomenon  is  included  herein,  since  there  is  some  evidence 
to  believe  that  thermo -osmosis  of  water  through  soils  is 
only  a  special  form  of  electro-osmosis  (Winterkorn  1947) . 

In  this  dissertation,  the  pertinent  physico-chemical 
and  electro-kinetic  phenomena  occurring  at  the  solid  water 
interfaces  are  described,  as  these  are  essential  for  an 
understanding  of  the  principles  of  electro-osmosis.  A 
comprehensive  review  of  the  literature  on  electro-osmosis 
has  already  been  published  (L.  Casagrande,  1953) .  Only 
the  historical  developments  are  traced  herein  for  the  sake 
of  maintaining  the  continuity  of  presentation;  the  reader 
is  referred  to  Dr.  Casagrande ls  review  for  more  detailed 


information. 


3.2.  Physico-chemical  phenomena  at  interfaces: 

3.2 .1.  Free  energy  of  surfaces: 

The  fundamental  property  of  liquid  surfaces  is 
that  they  tend  to  contract  to  the  smallest  possible  area, 
consistent  with  the  external  conditions  of  forces  acting 
on  the  liquid.  In  text -books  on  physical  chemistry  and 
surface  chemistry ,  this  tendency  is  usually  explained 
by  the  hypothesis  that  every  surface  molecule  of  the 
liquid  is  subject  to  a  strong  inward  attraction  by  the 
interior  molecules  and  that  the  surface  molecules  are 
continually  moving  inwards  more  rapidly  than  others  move 
outwards  to  take  their  places;  the  contraction  of  the 
surface  continues  until  the  maximum  possible  number  of 
molecules  are  in  the  interior ,  that  is,  until  the  surface 
is  the  smallest  possible  for  a  given  volume  (see,  for 
example,  Adam  1941) . 

On  account  of  the  foregoing,  it  is  evident  that 
work  must  be  done  to  extend  the  surface,  that  is,  to  bring 
molecules  from  the  interior  to  the  surface  against  the 
inward  attractive  forces.  Since  the  molecules  have  definite 
sizes,  there  will  be  a  definite  number  of  them  in  the  surface 
and  so,  the  work  done  in  extending  the  surface  will  also  be 
definite.  This  definite  quantity  of  work  is  termed  the  free 
energy  of  the  surface.  Although  the  foregoing  explanation 
is  concerned  with  liquid  surfaces,  it  is  obvious  that 
similar  reasoning  applies  to  solid  surfaces  also.  This 
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tendency  of  surfaces  to  contract  spontaneously  gives  rise 
to  the  important  phenomenon  of  adsorption. 

3.2.2.  Adsorption: 

Since  any  surface  tends  to  have  as  small  an  area 
as  possible,  it  follows  that  the  free  energy  of  the  surface 
also  tends  to  be  as  small  as  possible  (Preece  1947) .  In  a 
pure  substance,  minimum  free  energy  can  be  attained  only 
by  having  a  minimum  surface  area;  however,  if  more  than  one 
substance  is  present,  an  alternative  means  of  reducing  the 
surface  free  energy  is  available.  Of  the  various  species 
of  molecules  present,  the  ones  having  the  smallest  intrinsic 
energy  may  accumulate  at  the  surface,  with  molecules  having 
greater  intrinsic  energies  in  the  interior  (Adam  1947) . 

This  phenomenon  -  of  preponderance  of  one  type  of  molecules 
at  the  surface  -  is  called  adsorption;  an  atom  or  a  molecule 
of  the  adsorbate  saturates  some  of  the  unbalanced  forces 
of  the  surface,  thereby  decreasing  the  free  energy  of  the 
surface.  In  soil-water  systems,  dissociated  ions  correspond 
to  the  ‘molecules1,  while  the  solid  surface  becomes  the 
‘adsorbent*  (Preece,  1947) . 

3.2.3.  Formation  of  the  double-layer: 

The  presence  of  adsorbed  ions  on  the  soil-water 
interfaces  gives  rise  to  the  formation  of  a  layer  of 
counter-ions,  a  short  distance  away  from  the  interface.  This 
formation  of  a  layer  of  counter-ions  is  a  direct  consequence 
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of  the  attractive  and  repulsive  forces  exerted  by  the 
individual  ions.  The  adsorbed  layer  of  ions  and  the  layer 
of  counter-ions  are  collectively  termed  as  the  electric 
double  layer.  The  presence  of  a  layer  of  anions  (negative 
ions)  in  close  proximity  to  a  counter  layer  of  cations 
(positive  ions)  necessitates  that  a  potential  difference 
act  across  the  double  layer ,  opposing  the  electrostatic 
forces  of  attraction  and  repulsion:  this  potential  difference 
across  the  electric  double  layer  is  called  the  zeta  potential 
or  the  electrokinetic  potential. 

3.3.  Electro-kinetic  phenomena: 

The  presence  of  the  zeta  potential  across  the  double 
layer  maintains  the  same  in  equilibrium,  as  long  as  the 
solid-liquid  interface  is  not  disturbed.  If  such  a 
disturbance  occurs,  a  series  of  electrokinetic  effects  are 
produced,  depending  upon  the  nature  of  the  disturbance. 

Thus,  if  a  liquid  is  forced  through  a  capillary  by 
hydrostatic  pressure,  a  potential  difference  is  set  up 
between  the  ends  of  the  capillary,  opposing  flow  (Streaming 
potential) ;  if  charged  particles  are  allowed  to  move  through 
a  liquid  by  mechanical  means,  a  potential  gradient  is  set  up, 
opposing  the  motion  of  the  particles  (Sedimentation  potential) 
If  the  disturbing  force  is  of  electrical  origin  the  converse 
of  these  effects  occurs  (Butler,  1951) .  Thus,  an  electric 
field  applied  axially  to  a  capillary  containing  liquid  causes 
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the  liquid  to  flow  (electro-osmosis)  7  an  electric  field 

applied  to  a  suspension  containing  charged  particles 

causes  the  particles  to  move  (Cataphoresis  or  electrophoresis) . 

These  various  effects  were  observed  and  reported 
by  Quincke ,  Dorn,  Reuss  and  Reuss  respectively  (Butler,  1951) . 
Helmholtz  furnished  "the  classical  theory  for  all  the  primary 
electrokinetic  effects"  assuming  a  rigid  double  layer,  that 
is,  one  in  which  "there  is  a  plane  with  a  uniform  distribution 
of  charge  of  density  teI  per  unit  area,  separated  by  a 
distance  1dI  from  a  similar  plane  with  charge  density  ’-e1 
(Butler,  loc.  cit.) . 

3.3.2.  Structure  of  the  double  layer  and  consequences: 

Helmholtz's  assumption  of  a  rigid  double  layer  has 
been  criticised  on  the  grounds  that  the  electrical  forces 
causing  the  formation  of  the  double  layer  and  the  osmotic 
forces  tending  to  maintain  homogeneity  of  the  solution 
oppose  each  other  and  thus  there  can  be  no  abrupt  change  in 
the  concentration  of  any  kind  of  ions  in  the  vicinity  of 
the  double  layer,  but  "merely  a  gradual  increase  of 
concentration  of  ions  of  one  sign  and  a  rarefaction  of 
ions  of  the  other  sign" .  A  few  modified  theories  of  the 
structure  of  the  double  layer  have  been  put  forward  (see 
Butler  1951,  Chapter  II) .  However,  the  simple  Helmholtz's 
theory  has  been  found  to  be  sufficiently  accurate  for  soil 
moisture  studies  involving  moisture  contents  well  above  the 
plastic  limit  of  the  soil,  if  the  applied  potential  gradient 
is  small  (Winterkorn  1947) . 
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According  to  this  theory ,  as  already  mentioned,  the 
double  layer  can  be  visualized  as  a  rigid  inner  layer  of 
water  molecules  and  anions  on  the  surface  of  the  soil  grain, 
and  a  movable  outer  layer  of  water  molecules  and  cations. 
Under  the  influence  of  an  externally  applied  electric 
potential  gradient,  the  cations  of  the  outer  layer  will 
migrate  towards  the  cathode,  shearing  the  double  layer  and 
carrying  the  water  molecules  by  viscous  drag  (Thompson  1942) . 
This  effect  is  termed  "electro-osmosis" . 

3.4.  Application  of  electro-osmosis  to  soil-water  systems: 

The  electro-osmotic  velocity  of  flow  is  given  by 
the  expression: 

_  D  £  E 

V  -  4  h  p.1 

where  v  is  the  electro-osmotic  velocity  through  a  capillary 
tube , 

D  is  the  dielectric  constant  of  water, 

is  the  electro-kinetic  (zeta)  potential, 

E  is  the  applied  electric  potential, 
p.  is  the  viscosity  of  water,  and 
1  is  the  length  of  the  capillary. 

(For  a  derivation  of  this  equation  from  Helmholtz's  theory, 
see  Winterkorn  1947,  Bull  and  Gortner  1931,  Leo  Casagrande 
1949)  ,  This  equation  may  be  reformulated  (Winterkorn  1947) 
similar  to  Darcy's  law,  yielding 

i  E 

v  =  ke  •  j 

wherein  ke  is  the  composite  electro-osmotic  permeability.  If 
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one  assumes  that  the  soil  mass  may  be  represented  by  a 
bundle  of  circular  capillaries  of  the  same  length  and  cross 
section/  the  rate  of  electro-osmotic  discharge  may  be  obtained 
as 


where ,  qg  =  the  rate  of  discharge  of  water  per  unit  area, 

r  =  the  radius  of  the  capillary,  and 

b =  the  number  of  capillaries  per  unit  area. 

r  ^ 

Treating  the  quantity  (  '  )  as  a  two  dimensional 

expression  of  the  porosity  of  the  soil,  n#  the  above  equation 
reduces  to 


It  may  be  inferred  from  the  above  equation  that  the 
rate  of  electro-osmotic  discharge  is  proportional  to  the 
porosity  of  the  soil,  the  electro-osmotic  permeability  and 
the  applied  electric  potential  gradient.  The  electro-osmotic 
permeability,  in  its  turn,  is  a  function  of  the  zeta 
potential  across  the  double  layer,  the  dielectric  constant 
of  the  water  and  its  viscosity.  These  three  quantities  are 
influenced  to  a  greater  or  less  extent  by  the  electrolyte 
concentration,  by  the  vicinity  of  adsorptive  force  fields, 
and  by  the  temperature  of  the  system  (see  Winterkorn  1947) „ 

In  particular,  with  high  concentrations  of  electrolytes, 

"the  zeta  potential  can  drop  to  zero  .  .  „  or  the  charges  in 

the  double  layer  may  become  reversed  so  that  the  direction 
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of  flow  will  also  be  reversed"  (Casagrande  1951) .  As  a 
result  of  a  series  of  experiments,  Dr.  Casagrande  gives  the 
value  of  0.5  x  10^  cm. /sec.  per  unit  potential  gradient  for 
the  electro-osmotic  permeability  of  most  soils;  however, 
Winterkorn  (loc.  cit.)  has  reported  experimental  results 
suggesting  that  the  value  varies  with  the  moisture  content 
of  the  soil  and  the  applied  electric  potential  gradient. 

The  reliability  of  these  results  has  been  questioned  by 
Casagrande  (pp.  18,  1953) . 

3.4.2.  Pore  pressures  induced  by  electro-osmosis: 

According  to  the  Helmholtz *s  theory,  the  movable 
part  of  the  double  layer  carries  the  molecules  of  the  bulk 
water  by  viscous  drag  and  as  such  no  stresses  should  be 
developed  in  the  pore  water.  Tests  by  Casagrande  (1951) 
reveal  that  "electro-osmosis  does  create  in  compressible 
fine  grained  soils  tension  in  the  pore  water,  and 
compressive  forces  and  a  corresponding  amount  of  consolidation 
in  the  grain  skeleton  ..."  However,  Casagrande  states 
that  this  question  has  not  been  specifically  investigated 
or  discussed  in  the  literature. 

3.5.  Uses  of  electro-osmosis  in  soil  mechanics. 

3.5.1.  Temporary  stabilization  with  osmotic  drainage: 

Many  practical  cases,  where  electro-osmotic  drainage 
of  “treacherous1  soils  has  been  employed  during  construction 
periods  have  been  described  (Casagrande  1951,  Graham  1951, 
Hardy  1953,  Loughney  and  Hickey  1959) .  In  all  these  cases. 
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the  actual  decrease  in  water  content  of  the  soil  due  to 
electro-osmosis  was  found  to  be  relatively  small;  but  even 
this  small  decrease  in  water  content  was  found  to  be  ample 
for  the  purpose  of  stabilizing  the  construction  site. 

Thus,  it  is  necessary,  before  making  a  final  choice  of  the 
drainage  procedure,  to  assess  the  feasibility  of  using 
electric  drainage  -  by  means  of  laboratory  tests  to  determine 
the  osmotic  permeability,  natural  water  content  -  strength 
relationship,  electrical  resistivity  and  similar  properties 
of  the  soil.  In  some  cases,  full  scale  field  tests  also 
have  been  employed  (Casagrande  1951) . 

3.5.2.  Reduction  of  skin  friction  of  cathodic  piles: 

At  the  beginning  of  electro-osmotic  drainage  of  a 
soil,  before  steady-state  conditions  are  established,  the 
rate  of  discharge  at  the  cathode  is  somewhat  greater  than 
the  rate  at  which  water  is  supplied  to  the  cathodic  region, 
thus  resulting  in  a  partial  drying  of  the  soil  in  the 
immediate  vicinity  of  the  cathode.  (In  due  course,  however, 
this  condition  may  reverse,  resulting  in  an  excess  of  water 
at  the  cathodic  region  depending  upon  the  drainage 
facilities) .  During  this  short  period  of  desiccation, 
shrinkage  cracks  may  develop  in  the  soil  in  the  immediate 
vicinity  of  the  cathode,  thus  reducing  its  coefficient  of 
skin  friction.  Begemann  (1953,  also  included  in  Casagrande ts 
review)  explored  the  possibility  of  reducing  the  skin  friction 
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of  piles  in  clayey  soil  by  subjecting  them  to  an  electric 
potential  for  a  short  period  of  time  (not  exceeding  a  few 
minutes) ;  he  found  that  the  friction  of  the  cathodic  piles 
showed  a  marked  reduction  due  to  this  treatment  and  that 
prolonged  treatment  caused  a  recovery  of  the  friction.  No 
large  scale  practical  application  of  this  technique  has 
been  reported  so  far. 

3.5.3.  Permanent  hardening  of  clays: 

Since  electro-osmotic  drainage  of  water  involves 
migration  of  ions  through  the  soil,  it  is  possible  to 
introduce  particular  types  of  cations  at  the  anode  which 
will  move  through  the  soil  and  be  fixed  on  the  soil  surface 
by  the  mechanism  of  ion  exchange,  thus  altering  the 
physico-chemical  properties  of  the  soil.  The  introduction 
of  the  desired  cations  at  the  anode  may  be  achieved  by 
using  anodes  containing  those  cations.  Thus,  when  aluminum 
anodes  are  employed, Al-cations  are  stripped  away  from  the 
anode  and  are  either  adsorbed  on  the  clay  surface  or  form 
aluminum  hydroxide  on  the  surface  (Shukla,  1953) .  Both  of 
these  results  cause  an  irreversible  hardening  of  the  soil, 
the  former  on  account  of  the  physico-chemical  modification 
of  the  soil  constituent,  and  the  latter,  presumably  due  to 
the  binding  action  of  aluminum  hydroxide. 

This  method  has  been  successfully  employed  for 
increasing  the  bearing  capacity  of  piles  (see  papers  by 
Casagrande,  Spangler  and  King,  and  others,  reviewed  by 
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Casagrande) .  In  all  cases,  it  has  been  found  that  only  a 
small  portion  of  the  soil  surrounding  the  aluminum  anodic 
pile  has  been  altered  by  the  treatment.  However,  Murayama 
and  Mise  (1953)  have  reported  evidence  indicating  that  the 
location  of  the  hardened  zone  can  be  controlled  artificially 
by  suitably  altering  the  salinity  of  the  soil. 

3.5.4.  Detrimental  effects  of  electro-osmosis: 

These  undesirable  effects  are  the  result  of  accidental 
discharge  of  electric  current  in  a  soil  mass;  such  accidental 
discharge  of  current  may  take  place  due  to  improper  insulation 
of  underground  electric  cables  or  of  similar  electrical 
appliances;  notable  among  the  latter  case  are  the  switch 
boxes  of  railway  station  yards  using  some  system  of  track- 
circuiting.  Depending  upon  the  nature  of  the  soil  involved, 
these  accidental  discharges  may  cause  either  an  accumulation 
of  water  in  the  soil  resulting  in  loss  of  strength,  or  a 
desiccation  of  the  soil,  causing  shrinkage  cracks. 

3.6.  Outlook  for  the  use  of  electro-osmosis: 

Winterkorn  (1947)  has  expressed  the  hope  that  by 
using  perforated  anodes,  a  constant  supply  of  some  liquid 
dielectric  may  be  maintained  at  the  anode,  to  replace  the 
moving  water  and  to  assure  a  continuous  removal  of  water 
without  having  to  resort  to  high  potential  gradients.  If 
such  an  arrangement  proves  feasible,  it  may  be  utilized  in 
introducing  the  so-called  moisture  barriers  into  subgrades 
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and  eliminate  accumulation  of  water  beneath  pavements  and 
frost  damage. 

The  actual  decrease  in  water  content  of  a  soil 
subjected  to  electro-osmosis  is  very  small  and  thus^ 
electro-osmosis  alone  is  not  economically  feasible  for 
large  scale  consolidation  of  soil  masses.  Casagrande  (1949) 
has  expressed  the  view  that  it  may  be  possible  to  introduce 
a  few  electrodes  at  the  bottom  of  the  deposit  and  rely 
upon  osmotic  action  to  mobilize  the  full  weight  of  the 
overlying  soil  and  water  for  the  purpose  of  consolidation. 
The  possibilities  of  this  technique  have  not  apparently 
been  explored  as  yet. 
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CHAPTER  IV 


FLOW  OF  WATER  UNDER  THERMAL  GRADIENTS 


4.1.  Introduction: 

Ground  temperature  investigations  show  that  the 
temperature  in  the  region  immediately  below  the  surface 
varies  with  depth  and  with  the  daily  and  monthly 
fluctuations  of  the  surface  temperature ,  and  that  below 
a  certain  critical  depth,  the  temperature  is  sensibly 
independent  of  the  surface  temperature  fluctuations 
(Crawford  and  Legget,  1955) .  Above  this  critical  depth, 
the  daily  fluctuations  of  the  surface  temperature  give 
rise  to  a  reversal  of  the  direction  of  the  temperature 
gradient  every  day;  in  countries  where  the  monthly  average 
temperature  varies  considerably  from  the  annual  mean 
temperature,  the  monthly  fluctuations  of  the  surface 
temperature  also  give  rise  to  a  reversal  of  the  temperature 
gradient . 

If  substantial  quantities  of  water  may  be  shown  to 
migrate  in  response  to  the  prevalent  temperature  gradient, 
the  implications  of  the  above  noted  thermal  gradients  become 
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obvious:  assuming,  in  accordance  with  the  potential  theory, 
that  the  quantity  of  water  moved  in  proportional  to  the 
potential  gradient  causing  flow  and  the  duration  of 
application  of  the  gradient,  a  net  accumulation  of  water 
on  one  side  will  be  expected,  if  the  durations  of  the  two 
opposing  temperature  gradients  are  not  equal.  Field  data 
reported  by  Benkelman  and  Williams  reveal  that  the  durations 
of  the  temperature  gradients  in  a  subgrade  are  indeed  not 
equal.  A  typical  set  of  data  reported  by  these  authors  is 
reproduced  here  as  an  example;  as  will  be  seen  from  this 
plot,  the  temperature  at  the  top  of  the  base  is  higher  than 
that  at  the  bottom  of  the  base,  only  for  about  9  hours  per 
day.  Thus,  there  is  ample  justification  for  inquiring 
into  the  mechanism  of  moisture  migration  in  soils  due  to 
thermal  gradients  and  for  ascertaining  the  extent  to  which 
such  a  migration  may  account  for  the  observed  accumulation 
of  water  beneath  impervious  pavements. 

In  the  following  sections,  the  available  literature 
on  this  subject  is  reviewed  in  some  considerable  detail. 

The  various  quantitative  methods  of  analyses  are  briefly 
described  and  evaluated.  Finally,  the  present  state  of 
our  knowledge  of  the  mechanisms  is  summarized. 

* 


"A  co-operative  study  of  Structural  Design  of  Non-rigid 
Pavements",  Highway  Res.  Board  SR  46,  1959. 


' 

.  '  . .  ■  :  J  / 

, 

v 

'  : 

l 


|- 


\  '■ 

•  •  •  v  -  . 

i  '  " 

■ 

'  x 

* 

'  ■ 

:  ’  :  :  ■ 


■  \  .  ^  : 


59 


EASTERN  STANDARD  TIME- 5- OURS 

“Plate  1.  TEMPEV  ATJt:E  AT  \7,  RASE  COURSE  CECTIOM  [l»  OAV  AVERAOe] 


(  9 


4.2.  Review  of  Literature 

* 

4.2.1.  Historical 

In  1907,  Lippman  observed  that  when  a  temperature 
gradient  was  applied  across  a  membrane  of  gelatine  separating 
two  sections  of  a  body  of  water,  the  water  moved  from  the 
cold  to  the  hot  side,  and  that  the  observed  rate  of  flow 
was  about  50  mg.  per  minute  for  a  temperature  difference 
of  80°C ;  he  called  this  the  thermo-osmosis  of  liquids 
through  porous  membranes.  Subsequently,  Auburt,  in  1912, 
found  that  the  direction  of  flow  depended  upon  the  nature 
of  the  membrane:  gelatine  and  pig *3  bladder  gave  rise  to 
flow  from  the  cold  to  the  hot  side,  while  parchment  paper 
and  viscose  gave  rise  to  flow  in  the  reverse  direction.  He 
also  found  that  porous  plugs  of  inorganic  materials  gave 
rise  to  ho  flow;  that  the  thermo -osmotic  flow  depended 
upon  the  presence  of  water-soluble  substances  in  the 
membrane;  and  that  the  flow  was  not  caused  by  the  diffusion 
of  these  substances  into  the  water  on  either  side,  but  the 
flow  took  place  only  when  these  substance  were  still  present 
in  the  membrane.  These  findings  led  Auburt  to  postulate 
that  the  phenomenon  was  electrical  in  origin  and  related 
to  electro-osmosis. 

In  1941,  considerable  quantities  of  flow  of  water, 
carbon  tetrachloride  and  other  liquids  through  sintered  glass 
plates  under  a  thermal  gradient  was  reported  by  Derjaquin 
and  Sidorenkov.  This  is  at  variance  with  Auburt “s  findings, 
who  could  detect  no  flow  through  inorganic  materials. 

* 

The  material  for  this  section  was  largely  drawn  from  the 
paper  published  by  Hutchison  et  al .  (1948) . 
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Hutchison  et  al .  (1948)  conducted  a  series  of 

experiments  to  check  the  data  reported  by  Derjaquin  and 
Sidorenkov  and  found  that  there  was  no  flow  through  a 
sintered  glass  porous  plate,  if  the  plate  and  the  liquid 
were  chemically  pure 7  flow  occurred  only  when  electrolytes 
were  present  in  the  porous  plate.  They  suggested  that 
the  movement  observed  by  the  Russian  investigators  might 
be  due  to  simple  thermal  expansion  of  the  liquid  or  to 
the  presence  of  impurities.  These  authors  derived  an 
equation  for  the  thermo-osmotic  pressure  in  terms  of  the 
heat  of  transfer,  molar  volume  of  the  substance  and  the 
mean  temperature;  Winterkorn  (1955)  rederived  this  equation 
from  a  somewhat  different  line  of  reasoning,  and  utilized 
it  in  a  Darcy  type  of  equation  for  use  in  soil  moisture 
studies . 

No  further  work  on  pure  thermo-osmosis  has  been 
reported.  Thus,  the  available  data  reveal  that  thermo-osmosis 
does  not  occur  through  mineral  (inorganic)  porous  substances, 
except  under  the  presence  of  electrolytes  in  the  membrane. 

As  mentioned  in  an  earlier  chapter,  soil  grains  usually 
have  an  ionosphere  surrounding  them,  and  so  thermo -osmosis 
should  be  possible  in  soils,  even  though  the  solid  constituent 
are  purely  inorganic  in  nature. 

4.2.2.  Thermo-osmosis  in  soils: 

The  effect  of  thermal  gradients  on  soil  moisture 
migration  was  reportedly  studied  by  Bouyoucos  as  early  as 
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* 

1915  .  His  extensive  experiments  revealed  that  there  is  a 
flow  of  moisture  from  the  hot  to  the  cold  side  in  a  closed 
soil  column;  that  there  is  an  optimum  initial  moisture 
content  for  maximum  thermo -osmotic  flow;  and  that  this 
optimum  moisture  content  falls  in  the  vicinity  of  the 
plastic  limit  of  the  soil.  Working  with  soil  columns 
separated  into  sections  by  air  gaps,  Bouyoucos  observed 
that  the  net  transfer  of  moisture  from  the  hot  to  the  cold 
end  was  considerably  reduced.  From  these  findings,  he 
concluded  that  the  flow  took  place  on  account  of  the 
difference  in  capillary  potential  between  the  two  ends, 
and  that  the  flow  was  mainly  in  the  liquid  phase. 

At  about  the  same  time,  Lebedeff  in  U.S.S.R. 
observed  substantial  upward  migration  of  water  in  field 
soils;  after  accounting  for  precipitation  and  evaporation, 
he  computed  that  a  volume  of  water  equivalent  to  66  mm. 
depth  per  unit  area  had  migrated  upwards  from  the  ground 
water  table  during  the  winter  1914-15  (Lebedeff  1927) . 

He  ascribed  this  migration  to  the  phenomenon  of  evaporation 
from  the  water  table  and  condensation  at  the  surface 
layers  and  implied  that  the  flow  was  the  result  of  diffusion 
of  water  vapour  under  a  vapour  pressure  gradient. 


The  original  paper  published  by  Bouyoucos  was  not 
available  to  the  writer  at  the  time  of  preparing  this 
dissertation;  Bouyoucos "  work  has  been  quoted  by  many 
of  the  subsequent  investigators  (see  Winterkorn  1947, , 
Rollins  et  al.  1954) . 
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In  1939 ,  Smith  encountered  the  thermo-osmotic 
phenomenon  when  he  was  interested  in  the  determination  of 
thermal  conductivities  of  moist  soils.  When  the  initial 
moisture  content  of  the  sample  was  such  as  to  give  75% 
relative  humidity  for  the  soil  atmosphere^  the  specimen 
appeared  to  be  of  uniform  consistency  after  the  application 
of  a  thermal  gradient;  when  the  relative  humidity  was 
raised  to  99. 3%*  there  appeared  a  distinct  line  of  colour 
change,  indicating  a  heterogeneous  consistency  in  the 
final  specimen.  Smith  investigated  the  matter  further  by 
determining  the  final  moisture  contents  at  three  locations 
along  the  specimen  for  various  soils  and  obtained  results 
very  similar  to  those  obtained  by  Bouyoucos  earlier.  He 
reasoned  that  the  changes  in  surface  tension  of  the  water 
induced  by  the  applied  thermal  gradients  were  too  small 
to  account  for  the  transfer,  and  that  vapour  diffusion 
under  the  imposed  vapour  pressure  gradients  was  too  slow 
a  process  to  account  for  the  rapidity  with  which  equilibrium 
was  obtained.  He  postulated  that  the  phenomenon  involved 
a  mass  convection  of  the  water  vapour  across  the  pores 
(Smith  1939) . 

Subsequent  investigations  (Smith  1940)  led  Smith 
to  alter  his  earlier  ideas  and  postulate  that  the  actual 
mechanism  involved  evaporation  and  recondensation  of 
water  in  the  soil  pores  nearest  to  the  warm  end,  which 
caused  a  disturbance  in  the  equilibrium  of  the  capillary -he Id 
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water  in  the  adjoining  pores,  thus  causing  migration  to 
the  cold  side . 

Winterkorn  (1947)  lists  four  main  theoretical 
concepts  that  may  be  applied  to  thermo -osmotic  phenomenon 
as  follows: 

1.  diffusion  of  water  vapour  through  the  pore  space  under 
a  concentration  or  partial  pressure  potential; 

2.  diffusion  of  water  in  solid  solution*  similar  to  the 
case  observed  for  hydrophilic  membranes  separating  chambers 
of  different  water  vapour  pressure; 

3.  movement  as  capillary  water  due  to  the  difference  in 
surface  tension  at  different  temperatures;  and 

4.  flow  in  the  film  phase  along  the  internal  surface  of  the 
soil  due  to  change  in  water  affinity  with  change  in 
temperature „ 

After  discarding  the  first  and  the  third  of  these 
concepts  as  of  insignificant  importance*  he  proposes  that 
concept  (4)  should  be  applicable  to  soil-water  systems  and 
states"  ....  the  exchangeable  cations  possess  a  greater 
activity  (or  are  more  dissociated)  at  the  cold  than  at  the 
warm  side;  they  cannot  move  to  the  warm  side  because  they 
are  held  by  the  negative  charges  of  the  mineral  surfaces. 

The  only  way  to  decrease  this  concentration  potential  is 
by  movement  of  water  from  the  warm  to  the  cold  side" .  Based 
on  this  qualitative  explanation*  he  developed  a  quantitative 
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method  of  computing  the  amount  of  thermo-osmotic  movement 

in  soils;  the  thermo-osmotic  permeability  computed  in 

—  8 

this  manner  had  a  value  of  3.2  x  lO  cm. /sec.  for  a 

temperature  gradient  of  1°C  per  cm.,  while  the 

-7 

experimental  value  was  found  to  be  2.8  x  10  cm. /sec.  for 

the  same  conditions.  Such  a  correspondence  was  considered 
to  be  close.  In  the  same  paper,  Winterkorn  made  a  brief 
mention  about  the  detection  of  electric  potentials  between 
the  hot  and  cold  sides  of  the  specimen  and  stated  "the 
available  evidence  indicates  that  thermo-osmosis  is  a 
special  form  of  electro-osmosis"  (compare  this  with  Auburt5s 
findings  mentioned  earlier) . 

Jones  and  Kohnke  (1952)  studied  the  relation 
between  water  vapour  diffusion  and  soil  moisture  tension  (&„) 

r 

by  applying  a  sharp  temperature  difference  for  8  hours  at 
the  middle  of  soil  columns  contained  in  6 -inch  test  tubes 
and  determining  the  final  moisture  profiles.  From  these 
experiments,  these  authors  concluded  that  the  movement  of 
water  vapour  in  the  soils  studied  increased  rapidly  with 
the  initial  and,  after  a  certain  level,  the  rate  of 

r 

movement  dropped  off  sharply;  and  that  the  volume  of 
unsaturated  pores,  and  not  their  size,  governed  the  transfer 
of  water  vapour.  It  is  to  be  noted  that  these  authors  used 
a  temperature  gradient  as  a  means  of  establishing  a  vapour 
pressure  gradient  in  the  specimens;  as  such,  their 
investigations  are  not  primarily  concerned  with  thermo-osmosis 
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nevertheless,  their  conclusion  regarding  the  relative 
influence  of  the  overall  porosity  and  of  the  pore  size 
is  of  much  interest:  if  the  porosity  alone  is  the  governing 
factor,  soils  having  widely  different  pore  sizes  but  the 
same  porosity  will  be  expected  to  have  the  same  thermo- 
osmotic  permeability  (compare  this  with  the  similar  condition 
obtaining  in  the  case  of  electro-osmosis) .  Unfortunately, 
this  factor  has  not  been  systematically  investigated  as  yet, 
Gurr,  Marshall  and  Hutton  (1952)  reported  a  series 
of  experiments  on  thermo-osmosis  of  water  through  soils 
with  widely  different  initial  moisture  contents.  These 
authors  attempted  to  assess  the  extent  of  movement  in  the 
liquid  and  vapour  phases  by  tracing  the  movement  of  small 
quantities  of  salt  added  to  the  soil  prior  to  the  application 
of  the  thermal  gradient.  All  samples  were  subjected  to  a 
temperature  gradient  for  a  constant  period  of  5  days,  and 
then  the  salt  content  and  water  content  of  the  various 
sections  were  determined.  In  all  except  the  "driest  and 
wettest"  columns,  there  was  a  net  transfer  of  water  to  the 
cold  side  and  of  salt  to  the  hot  side.  The  authors 
concluded  that  the  thermo -osmotic  transport  of  water  occurred 
in  the  vapour  phase  from  the  hot  to  the  cold  side  of  the 
specimen  and  that  the  resulting  accumulation  of  water  at  the 
cold  end  caused  a  pressure  gradient  in  the  liquid  phase  in 
the  reverse  direction,  giving  rise  to  a  return  flow  of  water 
in  the  liquid  phase  to  the  hot  end. 
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Taylor  and  Cavazza  (1954)  reported  results  of  a 
series  of  tests  on  thermo -osmosis  of  water  through  a  silty 
loam  soil:  they  introduced  air  gaps  varying  from  1  mm. 
to  10  mm.  widths  between  adjacent  sections  of  the  specimen 
and  applied  temperature  gradients  varying  between  1.5  to 
1.7°C/cm.  They  found  that  a  stable  (equilibrium)  distribution 
of  moisture  could  be  produced  in  a  continuous  specimen 
after  a  sufficient  period  of  application  of  the  temperature 
gradient  while  such  a  dynamic  equilibrium  state  was  not 
obtained  during  the  same  period  in  the  segmented  specimen. 

They  reasoned  that  the  return  flow  in  the  liquid  phase  was 
responsible  for  the  establishment  of  steady  conditions  in 
the  former  case,  and  the  absence  of  it  caused  a  continuous 
movement  of  vapour  towards  the  cold  side  in  the  latter  case. 
The  writer  feels  that  had  the  experiments  been  continued 
for  longer  periods,  a  steady  equilibrium  condition  would 
have  been  produced  in  the  segmented  column  also,  since 
the  moisture  contents  of  the  various  segments  would  tend  to 
attain  equilibrium  with  the  prevalent  vapour  pressure  and 
temperature  conditions. 

Rollins,  Spangler  and  Kirkham  (1954)  constructed 
an  apparatus  in  such  a  way  that  the  water  accumulating  at 
the  cold  end  could  be  transferred  to  the  hot  end  of  the 
specimen  through  an  external  capillary  tube?  a  drop  of 
mercury  placed  in  the  external  capillary  tube  served  as 
a  means  of  indicating  the  direction  of  flow;  the  rate  of 
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movement  of  this  drop  could  also  be  taken  as  an  indication 
of  the  rate  of  thermo -osmotic  flow.  They  observed  that 
there  was  no  tendency  for  the  drop  to  move  from  the  warm 
side  to  the  cold  side  of  the  specimen  and  concluded  that 
liquid  flow  from  the  hot  end  to  the  cold  end  could  not  be 
of  significance.  With  open  systems  (water  movement  through 
the  external  capillary  tube  permitted)  they  observed  that 
the  accumulation  of  water  at  the  cold  end  was  not  appreciable, 
while  there  was  considerable  accumulation  in  a  closed, 
non-circulating  system.  Thus,  their  experiments  confirmed 
the  earlier  hypotheses  of  Gurr  et  al.  (1952)  and  of  Taylor 
and  Cavazza  (1954) . 

The  effect  of  density  and  air  content  of  the  samples 
on  the  thermo -osmotic  flow  rates  were  also  investigated. 

As  would  be  expected  from  a  qualitative  analysis,  the  flow 
rate  was  found  to  increase  with  increasing  air  contents 
and  decreasing  densities. 

The  influence  of  boundary  temperatures  and  of  the 
temperature  gradient  on  the  flow  rate  was  studied.  While 
the  flow  rate  increased  with  increasing  temperature 
gradients,  the  variation  was  not  linear.  Such  a  trend  needs 
further  verification  and  corroboration;  if  the  non-linear 
relationship  between  the  flow  rate  and  the  temperature 
gradient  is  definitely  established,  use  of  potential  theory 
to  describe  the  thermo-osmotic  phenomena  may  be  open  to 
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question,  since  the  theory  assumes  a  constant  proportionality 
factor  between  the  flux  and  the  potential  gradient. 

Hutcheon  (1955)  postulated  a  qualitative  line  of 
reasoning  to  account  for  the  existence  of  an  optimum 
moisture  content  for  maximum  thermo -osmotic  flow.  This 
line  of  argument  is  based  on  the  hypotheses  of  Gurr  et  al  „ , 
and  of  Taylor  and  Cavazza,  that  a  return  flow  of  liquid 
may  take  place  from  the  cold  to  hot  end  under  certain  conditions 
of  moisture  content. 

Swenson  and  Sereda  (1955)  described  the  construction 
of  an  apparatus  in  which  movement  of  water  through  soils, 
concrete  and  similar  porous  materials  could  be  studied 
keeping  the  moisture  content  of  the  specimen  constant;  this 
was  achieved  by  maintaining  suitable  suction  pressures 
across  porous  plates  in  contact  with  the  ends  of  the  sample. 
Results  of  preliminary  experiments  bn  Ottawa  day  confirmed 
the  findings  of  earlier  workers.  The  authors  state  "  .„. 
at  this  stage  it  is  not  at  all  evident  what  contribution  is  made 
by  the  existence  of  the  interface  between  the  sample  and  the 
porous  plate.  This  interface  might  even  make  the  main 
contribution  to  the  observed  flow.  Until  this  is  assessed 
by  further  experimental  work,  this  investigation  must  remain 
as  a  proposed  method  to  study  the  phenomenon  of  moisture 
movement  due  to  a  temperature  gradient" . 
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Winterkorn  (1955a)  presents  experimental  and 
theoretical  evidence  to  show  the  strongly  electric 
character  of  the  moist  soils  and  argues  that  if  electrical 
forces  could  produce  mechanical  effects  in  soil-water 
systems  (as  evidenced,  for  example,  by  electro-osmotic 
drainage  of  soils) ,  thermal  gradients  should  be  logically 
expected  to  have  electrical  consequences.  In  support 
of  this  deductive  reasoning,  he  cites  the  data  obtained 
by  Gowda  and  later  by  Williams,  both  of  whom  studied  the 
generation  of  electric  potentials  during  thermo -osmotic 
flow  through  natural  and  several  homoionic  modifications 
of  "New  Jersey  Hagerstown  soils",  H-bentonite,  H-kaolinite 
and  H-grundite.  Fairly  high  temperature  gradients  (2084 
to  4.03°C/cm.)  were  used  and  the  tests  were  run  for  periods 
not  exceeding  100  minutes.  According  to  the  author, 

"since  the  rate  at  which  equilibrium  is  approached  varies 
as  a  function  of  the  distance  from  the  equilibrium,  such 
data  (obtained  with  longer  durations  of  testing)  cannot 
serve  as  a  good  basis  for  the  calculation  of  thermo-osmotic 
transmission  coefficients  . . . " . 

In  all  experiments,  electric  potentials  were  measured 
and  it  was  observed  that  the  establishment  of  an  electric 
potential  was  "almost  instantaneous  with  application  of  hot 
and  cold  plates" ,  and  did  not  await  the  establishment  of  a 
uniform  temperature  gradient.  The  observed  electric 
potentials  varied  from  300  to  700  millivolts;  the  polarity 
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of  the  thermo-electric  potential  for  ‘normal1  soils  was  the 
same  as  that  of  applied  potential  to  cause  migration  of 
water  in  the  same  direction;  the  polarity  was  reversed  in 
the  case  of  "refractory1  Al-  and  Mg-  soils. 

As  for  the  mechanism  of  thermo-osmotic  transfer* 
the  movement  is  ascribed  to  vapour  flow  in  the  case  of  dry 
soils  (appreciably  drier  than  the  plastic  limit)  and  to 
liquid  flow  in  the  film  phase  in  the  presence  of  continuous 
adsorbed  films.. 

In  a  subsequent  paper  (1955b)  Winterkorn  derives 
an  equation  for  the  pressure  or  suction  associated  with 
the  transfer  of  one  gram  of  water  from  a  region  of  temperature* 
T,  to  another  at  T  +  dT;  the  final  equation*  based  on  the 
second  law  of  thermodynamics*  is  of  the  form: 


dP  = 


Q  dT 


(1) 


1  x  T' 

where  P  is  the  maximum  pressure  difference  associated  with  a 
temperature  gradient  dT/dl* 

Q  is  the  latent  heat  involved  in  the  change  of  restraint  of 
one  gram  of  water  incident  to  the  same  temperature  gradient* 
and  the  numeral  "unity1  represents  the  specific  volume  of 
water.  Substituting  this  value  in  a  Darcy  type  equation* 
Winterkorn  obtains  a  thermo-osmotic  flow  equation. 

Eq.  (1)  has  been  derived  from  a  different  view  point 
earlier  by  Hutchison  et  al .  (1948)  *  although  these  authors 

were  concerned  with  the  thermo-osmosis  through  membranes  and 


so  the  'latent  heat1  becomes  the  "heat  of  solution1  of  the 
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liquid  on  the  membrane.  But  the  use  of  this  equation  in  a 
Darcy  type  law,  as  is  done  by  Winterkown,  and  the  eventual 
formulation  of  a  mathematical  relationship  between  the 
thermo-osmotic  permeability  and  the  hydraulic  (Darcy's 
coefficient  of)  permeability  appear  to  be  gross 
oversimplifications  of  the  entire  phenomenon.  No 
experimental  support  of  this  relationship  has  been  adduced, 
either  by  Winterkorn  himself,  or  by  other  investigators. 

On  the  contrary,  data  furnished  by  Rollins  et  al  (1954) 
reveal  that  the  volume  rate  of  flow  under  thermal  gradients 
depends  upon  the  temperature  gradient,  as  well  as  upon 
the  boundary  temperatures;  in  such  case,  the  simple 
potential  theory  -  as  embodied  in  Darcy's  law  -  cannot  be 
assumed  to  hold. 

Habib  (1957)  reported  some  experiments  on  thermal 
transfer  of  moisture  in  saturated  and  unsaturated  soils; 
fully  saturated  specimens  were  tested  both  as  open  and 
closed  systems  separately.  A  continuous  flow  of  water  in 
the  former  case  and  the  establishment  of  positive  pore 
pressures  in  the  latter  case  were  observed.  The  magnitudes 
of  the  reported  pore  pressures  are  very  small,  not  exceeding 
50  gm./sq.  cm.,  but  according  to  the  author,  "the  pressure 
differences  (closed  system)  and  the  flows  (open  system)  are 
multiplied  by  three  or  four  when  the  interstitial  water 
contains  an  ionizable  salt" .  No  reasons  for  this  rather 
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significant  effect  of  salts  are  furnished  by  the  author; 
data  on  the  observed  flow  rates  have  also  not  been  reported. 
The  writer  is  inclined  to  believe  that  the  rates  could  not 
have  been  significant  and  might  be  attributable  to  flow 
under  a  surface  tension  gradient. 

Kuzmak  and  Sereda  (1958)  made  use  of  the  special 
apparatus  developed  earlier  by  Swenson  and  Sereda  (1955) 
(wherein  flow  of  water  through  porous  materials  could  be 
maintained  without  altering  the  moisture  content  of  the 
specimen)  and  studied  the  movement  of  water  under 
temperature  gradients  through  porous  plates  of  different 
textures  separated  by  gaps  of  different  widths Q  They 
found  that  a  flow  occurred  only  when  the  suction  pressure 
in  the  gap  was  increased  to  such  a  value  as  to  desaturate 
the  same,  thus  confirming  the  concept  of  vapour  flow 
under  a  thermal  gradient.  Salt  tracer  experiments  with 
a  sample  of  fine  sand  in  the  apparatus  revealed  that  the 
observed  transfer  of  salt  from  the  warm  to  the  cold  side 
could  be  completely  accounted  for  by  the  isothermal 
diffusion  of  salt.  Measured  streaming  potentials  during 
thermo -osmotic  flow  were  in  general  found  to  be  the  same 
as  the  values  observed  under  conditions  of  no  flow.  These 
two  observations  lent  support  to  the  earlier  conclusion  and 
thus  the  authors  concluded  that  vapour  flow  was  the  chief 


mechanism  involved. 
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The  writer  feels  that  a  more  sensitive  test  of  the 


concept  could  have  been  obtained  by  using  salt  solutions 
of  equal  solute  concentrations  on  both  sides,  thus  eliminating 
isothermal  diffusion. 

Hutcheon  (1958)  reports  very  extensive  experimental 
data  on  thermo-osmotic  flow  of  water  under  small  temperature 
gradients  (0.5  to  0.8°C/cm.);  in  some  of  his  experiments* 
the  cold  end  was  maintained  at  -4°C,  thus  freezing  the 
accumulated  water  and  preventing  a  return  flow  of  liquid 
from  the  cold  to  the  hot  side.  To  account  for  the  moisture 
profiles  obtained*  he  postulates  that  a  return  liquid  flow 
must  exist  even  in  the  warm  section  of  the  soil  column* 
transferring  water  from  the  interior  sections  to  the  warmest 
place,  that  is,  the  ''evaporation  sink1. 

4.2.3.  Quantitative  methods  of  analysis: 

DeVries  (1950)  assumed  the  validity  of  an  equation 
of  the  form: 

P  ‘d'H  ,  ^ 

q  =  -a  •  a  •  D  •  yw  •  p— -  •  (2) 

for  unidimensional  molecular  diffusion  of  water  vapour  through 
soils.  The  various  symbols  in  the  above  equation  were  defined 
as  follows: 

q  =  vapour  flux  in  g/cm^  sec., 

a  =  a  function  of  the  soil  structure,  dimensionless, 

a  =  volume  fraction  of  gas  filled  pores, 

2  . 

D  =  diffusion  coefficient  of  water  vapour  in  air,  cm  /sec 
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P  =  total  pressure  of  moist  air  760  mm.  Hg., 

•n  =  partial  pressure  of  water  vapour  in  mm.  Hg., 

x  =  coordinate  distance,  in  cm. 

Penman  (1940)  gave  the  empirical  value  a  =  0.66.  The  value 

of  the  diffusion  coefficient  was  assumed  to  be  given  by 

D  =4.42  x  10-4  T2.3  (3) 

P 

where  T  =  the  absolute  temperature  in  °K. 

When  this  equation  was  applied  to  predict  the 
movement  of  water  under  temperature  gradient,  it  was  found 
that  the  observed  quantity  of  flow  (or  the  rate  of  diffusion) 
was  considerably  greater  than  the  calculated  amount,  the 
ratio  of  the  two  rates  varying  from  three  to  eighteen 
(Hutcheon  1958,  Philip  and  deVries,  1957) .  This  discrepancy 
between  the  diffusion  theory  and  the  experimental  data  was 
ascribed  to  the  possible  presence  of  molecular  slip  flow, 
or  to  the  interaction  between  the  water  vapour  and  the  soil 
(Penman “s  coefficient  of  0.66  was  based  on  experiments 
with  carbon  dioxide  and  acetone,  neither  of  which  presumably 
reacted  with  the  solid  or  liquid  phase  of  the  medium;  see 
Philip  and  deVries,  1957) . 

A  very  plausible  explanation  of  the  observed 
discrepancy  is  furnished  by  the  interesting  experiments  of 
Woodside  and  Kuzmak  (1958) .  These  authors  constructed  an 
enlarged  “model1  of  the  solid  grain-pore  system:  6  inch 
diameter  marble  spheres  were  used  to  represent  the  soil 


,  .  ■  ■ 

:  .':o o'  ‘ :  w 

• 

nr.  . 

l  '  '  .  . 

. 

.  ■  '  .  ■■  ■  n 

■  * 

.  ;  .  j  '  '•  : 

7;  s 

i 

.S’.:  J  i: 

1  .  :  »;•.  ' 

. 

v 

' : : , i  .■ 

. 

*■ 

•  •  . 

j.  -  '  ■■■ 

. 

,  •  ' 

i';J-  ;  J;  u  .. 

V 

"S  :  '•  :  : 

.  •  .o'r  '  i  ':lx  j; 

. :  u  .  : 

: :  •:  '•  . 

.  .  .  o  {  :  ‘ /•:  •  ■ . 

1  ’ 

'■  ■  ; :  :  j v;  -  ;■■■.  -j  .  r  u  ■/:  . 


■-  t  :■  ...  I 


.  .  . 


\ 


.  .V:  ;■  •  : 


1 


grains,  bitumen  to  represent  the  interstitial  water  and 
silica  aerogel  to  represent  the  soil  air?  the  materials 
were  so  chosen  as  to  give  thermal  conductivity  ratios 
comparable  to  those  of  the  soil-water-air  system.  They 
found  that  the  temperature  gradient  across  an  air  filled 
pore  space  in  a  solid-air  system  was  about  six  times 
the  average  overall  temperature  gradient;  though  no 
actual  measurements  were  carried  out,  they  express  the 
opinion  that  the  ratio  will  be  less  than  six  in  a  soil- 
water-air  system  and  will  increase  with  decreasing  water 
contents.  In  the  light  of  these  findings,  these  authors 
propose  that  movement  of  water  under  a  temperature 
gradient  "involves  multiple  vapour  diffusion  and  liquid 
flow  steps  in  series"  and  that,  since  the  resistance  to 
liquid  flow  through  the  contact  wedges  is  negligible, 
the  rate  of  transfer  is  governed  by  the  rate  of  diffusion 
across,  the  pore  space.  The  validity  of  these  experiments 
and  of  the  conclusions  derived  therefrom  have  been 
seriously  questioned  later  by  deVries  and  Philip  (1959)  . 

Philip  and  deVries  (1957)  have  developed  a 
mathematical  theory  to  describe  thermo-osmosis,  recognizing 
that  the  transfer  involves  a  "series-parallel  flow"  through 
air-filled  pores  and  contact  wedges?  the  vapour  flux  is 
assumed  to  be  governed  by  the  vapour  pressure  gradient 

across  the  air-filled  pores,  and  the  liquid  flux,  on  account 
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of  the  negligibly  small  resistance,  is  assumed  to  adjust 
itself  to  the  vapour  flux.  The  final  equation  is  of  the 
form: 

q  -  (a  +  0)  D  .  -B-  .  3  (VT)  (4) 

-p-n  a  v 

where,  in  addition  to  the  symbols  already  defined. 


0  =  volumetric  moisture  content  of  the  soil, 

dP, 


dT 


=  the  rate  of  change  of  density  of  saturated 


water  vapour,  due  to  change  in  absolute 
temperature,  and 

( v T)  =  the  temperature  gradient  across  the  air-filled 

cl 

pores. 

Applying  this  equation  to  the  experimental  data 
reported  by  Rollins  et  al.  (1954)  and  Gurr  et  al .  (1952) , 

the  authors  find  that  the  agreement  is  "quite  satisfactory 
in  view  of  the  approximate  character  of  the  analysis" . 

4.3.  Present  Status  of  Our  Knowledge 
4.3.1.  Concerning  the  phenomenon: 

A  survey  of  the  available  literature  given  above, 
reveals  that  the  existence  of  a  thermo-osmotic  flow  through 
soils  has  been  definitely  established.  All  previous 
investigators  have,  in  general  observed 

(1)  that  there  is  a  transfer  of  water  from  the  hot  to  the 
cold  side  of  a  specimen; 

(2)  that  the  quantity  of  water  transported  is  a  minimum  in 
very  dry  and  very  wet  soils; 

(3)  that  there  exists  an  optimum  initial  moisture  content 
for  maximum  thermo-osmotic  flow;  and 

(4)  that,  in  fairly  moist  soils,  there  may  be  a  return  flow 
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of  water  from  the  cold  to  the  hot  side  of  the  specimen. 
These  observations  are  for  partially  saturated  soils  only. 
The  only  set  of  experimental  data  on  saturated  soils  is 
that  reported  by  Habib;  the  values  of  pore  pressures  in  a 
closed  system  observed  by  him  are  very  small  for  most 
practical  purposes;  flow  rates  in  open  systems  have  not 
been  reported. 

4.3.2.  Concerning  the  mechanism  of  thermo-osmosis: 

There  appears  to  be  a  considerable  amount  of 
controversy  as  to  the  nature  of  the  actual  mechanism 
involved  in  the  phenomenon.  In  general,  five  distinctly 
different  mechanisms  have  been  proposed  by  various  workers: 

(1)  liquid  flow  under  a  thermally  induced  capillary 
potential  (suction)  gradient  (Bouyoucos) ; 

(2)  evaporation  and  recondensation  of  water  vapour  in  the 
hot  end  which  gives  rise  to  a  liquid  flow  in  the  colder 
segments  due  to  capillary  potential  gradient  (Smith) ; 

(3)  osmotic  flow  in  the  film  phase  due  to  thermally 
induced  differential  water  affinity  of  the  soil  (Winterkorn) 

(4)  vapour  diffusion  to  the  cold  side  followed  by  a  return 
circulation  in  the  liquid  form  under  favourable  conditions 
(Gurr  et  al „ ,  Taylor  and  Cavazza,  Hutcheon  and  others) ;  and 

(5)  multiple  vapour  diffusion  and  liquid  flow  steps  in 
series  (Woodside  and  Kuzmak)  ,  and  a  series-parallel  flow  of 
vapour  and  liquid  (Philip  and  deVries) „ 
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It  is  very  difficult  to  state  definitely  which  of 
these  mechanisms  is  the  most  probable  one;  obviously  it  is 
quite  possible  for  all  of  them  to  be  involved  at  the  same 
time,  each  one  making  a  more  or  less  contribution  to  the 
total  flow.  The  critical  question  to  be  resolved  is:  what 
are  the  relative  contributions  of  each  of  these  mechanisms 
under  different  conditions  of  soil  moisture  content?  It 
is  not  possible  to  give  a  quantitative  answer  to  this 
question  at  the  present  stage  of  the  development;  numerous 
experiments,  conducted  under  widely  different  conditions, 
will  be  required  before  a  tentative  answer  could  be  proposed. 
4.3.3.  Concerning  the  quantitative  methods  of  analysis: 

The  inadequacy  of  the  simple  molecular  diffusion 
theory  has  already  been  mentioned.  The  extended  analysis 
proposed  by  Philip  and  deVries  (1957)  has  not  been 
systematically  verified  yet.  Even  so,  the  analysis  is  for 
the  case  of  dry  media  only,  where  there  is  no  continuity 
in  the  liquid  phase;  for  the  moist  media  with  liquid 
continuity,  they  propose  empirical  coefficients  to  account 
for  the  decreased  volume  of  the  pore  space  available  for 
vapour  flow,  Finally,  none  of  these  methods  of  analysis 
recognize  the  existence  of  a  return  flow  of  liquid  from 
the  cold  to  the  hot  side;  as  such,  they  can  be  used  only 
for  open  systems  where  recirculation  is  prevented. 
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4.3.4.  Concerning  the  physical  explanation  of  the  phenomenon: 

A  valid  physical  explanation  of  the  phenomenon 
will,  of  course,  depend  upon  the  nature  of  the  mechanism 
assumed  to  be  active.  As  already  seen,  Winterkorn  (1947) 
believes  that  the  flow  is  in  the  adsorbed  films  of  water 
along  the  internal  capillaries  of  the  soil;  he  ascribes 
the  cause  of  this  movement  to  the  greater  water  affinity 
of  soils  at  lower  temperatures  than  at  higher  temperatures. 

In  support  of  this  hypothesis,  he  cites  thermo -osmotic 
experiments  on  an  extracted  clay  where  consistent  electric 
potentials  were  observed  between  the  ends  of  the  specimen. 

As  already  mentioned,  when  a  liquid  is  forced 
through  a  fine  capillary  passage,  a  potential  gradient 
(streaming  potential)  is  set  up  between  the  ends  of  the 
capillary,  tending  to  oppose  the  motion.  Thus,  the 
existence  of  an  electric  potential  during  thermo-osmosis 
may  be  taken  as  a  strong  indication  of  the  presence  of  a 
flow  in  the  liquid  phase.  However,  Kuzmak  and  Sereda  (1958), 
working  with  sandy  soils,  report  that  no  streaming  potentials 
were  observed. 

The  other  physical  explanation  is  that  advanced  by 
Hutcheon  (1958) ,  based  on  the  'vapour  diffusion-recirculation 
in  the  liquid  phase1  concept.  Briefly,  his  arguments 
proceed  as  follows:  In  very  dry  soils,  the  rate  of  change 
of  relative  humidity  (and  hence  the  vapour  pressure)  of  the 
soil  air  due  to  changes  in  moisture  content  is  very  high; 
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consequently,  small  changes  in  moisture  contents  are 
sufficient  to  equalize  the  thermally-induced  vapour  pressure 
gradients.  The  total  transfer  is  hence  small,  and  a 
recirculation  in  the  liquid  phase  is  impeded  by  the  very 
low  unsaturated  permeability  of  the  soil.  As  the  initial 
moisture  content  increases,  the  rate  of  change  of  relative 
humidity  falls  off,  thus  necessitating  a  larger  change  in 
moisture  content  to  equalize  the  same  vapour  pressure 
gradient;  at  the  same  time,  the  air-filled  porosity  is 
decreasing,  thus  decreasing  the  rate  of  dissipation  of  the 
vapour  pressure  gradient.  At  the  cold  end,  the  rate  of 
condensation  and  consequently  the  rate  of  building  up  of  a 
suction  gradient  are  increasing,  due  to  the  increased 
relative  humidity  of  the  soil  atmosphere. 

Thus  the  system  now  involves  two  effects:  with 
increasing  moisture  contents,  the  rate  of  dissipation  of 
the  vapour  pressure  gradient  decreases,  while  the  rate  of 
building  up  of  a  suction  gradient  from  the  cold  to  the  hot 
side,  and  the  unsaturated  permeability  of  the  soil  are 
increasing.  It  is  apparent  that  there  must  exist  a  particular 
moisture  content  at  which  these  two  opposing  tendencies  will 
be  so  balanced,  as  to  give  a  maximum  net  transfer  of  water 
from  the  hot  to  the  cold  side.  At  very  high  moisture 
contents,  the  air  filled  pores  become  discontinuous  and  so 
vapour  diffusion  ceases.  The  only  movement,  under  these 
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circumstances ,  will  be  in  the  liquid  phase  under  the 
thermally  induced  suction  gradients ,  and  the  magnitudes 
involved  are  very  small. 

The  above  arguments  (presented  in  a  slightly 
modified  and  abbreviated  form)  appear  to  be  sound  and 
plausible;  above  all,  they  clearly  demonstrate  the 
existence  of  an  optimum  moisture  content  for  maximum 
thermo -osmotic  flow.  But,  if  the  primary  flow  is  in 
the  vapour  phase,  it  will  be  very  difficult  to  explain 
how  an  electric  potential  is  created  across  the  specimen 
during  thermo-osmosis:  it  is  not  established  -  to  the 
writer^s  knowledge  -  whether  liquid-gas  interfaces  may 
also  give  rise  to  electrokinetic  phenomena,  similar  to 
solid-liquid  interfaces;  if  they  do  not,  then,  a  vapour 
flow  relative  to  a  liquid  boundary  cannot  cause  a 
streaming  potential  to  be  created.  This  argument  presupposes 
that  a  streaming  potential  is  always  observed  during 
thermo-osmosis;  so  far,  the  data  reported  are  contradictory 
in  this  respect.  Winterkorn  and  his  co-workers  observed 
and  measured  consistent  values  of  electric  potentials  with 
clayey  soils,  while  Habib,  working  with  "mud  samples",  and 
Kuzmak  and  Sereda,  working  with  "fine  sand",  could  not  find 
any  consistent  voltage  or  polarity  across  the  samples.  Thus, 
the  question  is  still  unsettled. 
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4.3.5.  Concerning  similarities  between  electro-osmosis 

and  thermo-osmosis: 

Auburt^  findings  (see  Hutchison  et  al .  1948)  led 
him  to  postulate  that  the  phenomenon  of  thermo-osmosis  is 
electrical  in  origin  and  related  to  electro-osmosis. 
Subsequently,  Winterkorn  (1947,  1955a)  independently 
arrived  at  the  same  conclusion,  from  a  line  of  reasoning, 
entirely  different  from  the  one  employed  by  Auburt  earlier. 
He  also  found  that  an  electric  potential  was  created 
across  the  ends  of  a  specimen  of  clayey  soil  subjected  to 
thermal  gradients,  and  that  the  establishment  of  this 
electric  potential  was  almost  instantaneous.  It  is  true 
that  subsequent  investigators  (Habib  1957,  Kuzmak  and 
Sereda  1958)  have  been  unable  to  detect  any  consistent 
electric  potentials  during  thermo-osmosis,  but,  this 
need  not  necessarily  exclude  the  possibility  of  the  presence 
of  electric  potentials;  in  other  words,  owing  to  the  type 
of  soil  used  by  the  later  investigators,  the  induced 
electric  potentials  may  have  been  too  small  to  be 
measurable  by  ordinary  methods. 

Jones  and  Kohnke  (1952)  have  stated  that  the  volume 
of  air-filled  pores  and  not  the  size  of  the  pore  governs 
the  vapour  transfer  in  soils.  As  mentioned  in  an  earlier 
section,  this  condition  is  very  similar  to  the  one  obtaining 
in  electro-osmosis,  where  the  discharge  is  governed  by  the 
porosity  of  the  soil  and  not  by  the  individual  pore  size. 
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Thus,  the  available  evidence  indicates  that  it  is 
distinctly  possible  that  electro-osmosis  and  thermo-osmosis 
are  related  in  some  as  yet  unknown  manner. 
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CHAPTER  V 


EXPERIMENTS  ON  THERMO-OSMOTIC  FLOW  IN  SOIL 


5.1.  Scope  of  the  Investigations: 

Previous  investigators  have  assumed  that  in 
relatively  moist  soils,  thermo-osmotic  flow  gives  rise 
to  a  return  flow  of  liquid  water  from  the  cold  to  the 
hot  side  (Gurr  et  al.  1952,  Taylor  and  Cavazza  1954) , 
and  Hutcheon  (1955,  1958)  has  proposed  a  physical 
explanation  of  the  overall  phenomenon  accordingly. 

The  present  investigations  are  aimed  at  determining 
whether  such  a  flow  takes  place  or  not,  thus  verifying 
the  hypothesis  advanced  by  Hutcheon. 

5.2.  Experimental  Arrangement 
5.2.1.  Preliminary  considerations: 

If  the  soil  moisture,  transferred  to  the  cold  side 
by  thermo-osmosis,  tends  to  redistribute  itself  under  a 
suction  gradient,  it  is  obvious  that  the  moisture  profile 
of  the  specimen  should  be  time-dependent  during  the  early 
stages  of  the  flow.  Thus,  the  experimental  set-up  should 
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be  so  devised  as  to  subject  several  specimens  of  identical 
initial  conditions  to  varying  durations  of  application  of 
a  constant  temperature  gradient.  A  somewhat  similar 
testing  program  was  employed  by  Taylor  and  Cavazza  when 
they  studied  the  moisture  profiles  of  continuous  and 
segmented  columns  of  soil*  with  a  view  to  differentiate 
between  liquid  flow  and  vapour  flow;  however*  these 
authors  used  only  one  initial  moisture  content  that  gave 
an  equilibrium  profile  in  the  continuous  columns*  at  the 
end  of  seven  days.  In  the  present  series  of  investigations* 
six  different  initial  moisture  contents  were  employed. 

The  apparatus  used  in  this  investigation  is  broadly  similar 
to  the  one  employed  by  Hutcheon  (1958)  in  his  investigations „ 
5.2.2.  Soil  sample: 

About  85  lbs„  of  the  sample  were  obtained  by 
grinding  pure  silica  sand  so  that  the  final  material  will 
completely  pass  a  60-mesh  sieve.  Attempts  to  obtain  a 
finer  grinding  proved  to  be  uneconomical  and  time-consuming 
on  account  of  the  large  quantities  of  sample  involved  and 
of  the  extremely  high  resistance  of  the  soil. 

Subsequent  to  grinding*  the  sample  was  digested 
with  excess  dilute  hydrochloric  acid  for  a  period  of  72 
hours  so  as  to  remove  the  iron  as  ferric  chloride;  this 
was  followed  by  repeated  washing  with  distilled  water  to 
remove  the  salt.  A  vacuum  technique  was  used  to  filter  the 
distilled  water  out  of  the  sample.  After  washing*  the  sample 
was  over-dried  at  105  -  5°C. 
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5.2.3.  Soil  container: 

Commercially  made  44  mm.  "Lucite"  tubes  were 
employed  as  specimen  containers;  the  tubes  were  cut  into 
rings  2  cm.  long,  and  eight  such  rings  were  taped  together 
by  adhesive  tape  to  provide  a  container  for  each  specimen, 
thus  yielding  specimens  16  cm.  long  x  4„4  cm.  diameter. 
Necessary  holes  were  drilled  on  the  sides  of  these  rings 
for  introducing  thermo-couples. 

5.2.4.  End  plugs: 

The  end  plugs  were  made  of  open-ended  brass  rings 

that  would  fit  snugly  into  the  container  ring.  The  end 

1  " 

closest  to  the  specimen  was  covered  by  an  brass  disc 

soldered  in  position;  the  other  end  was  closed  by  a  similar 
plate  with  two  tubes  projecting  out.  Thus,  the  end  plugs 
could  be  used  as  thermal  contacts,  by  circulating  water 
at  the  desired  temperature  through  them. 

5.2.5.  Hot  and  cold  sources: 

Two  Y6  G.I.  tanks,  1*9"  x  1*9"  x  lt9",  with 
individual  capacities  of  approximately  34  gallons,  served 
to  hold  the  hot  and  cold  water  to  be  circulated  through 
the  end  plugs  of  the  specimens.  An  ordinary  laboratory 
freezing  unit,  operating  in  conjunction  with  a  sensitive 
bi-metallic  thermo-regulator  was  employed  for  maintaining 
constant  cold-end  temperature.  A  thermostatically  controlled 
(as  claimed  by  the  manufacturer)  domestic  heating  coil  was 
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used  in  the  hot  side;  soon  it  was  found  necessary  to  have 
an  auxiliary  thermo-regulator  on  this  side  as  well,  and 
so  another  bi-metallic  thermo-regulator  was  connected  in 
series  with  the  heating  coil,, 

5.2.6o  Circulation  system: 

Water  from  the  hot  and  cold  tanks  was  circulated 
through  the  respective  end  plugs  of  the  specimens  by 
means  of  two  water-throw  pumps;  all  specimens  were  arranged 
in  parallel  (see  schematic  diagram  of  the  lay-out)  »  Initially, 
the  pumps  were  so  arranged  as  to  suck  water  from  the  tank, 
force  it  through  the  battery  of  specimens  and  then  discharge 
the  same  back  into  the  tank.  Preliminary  tests  revealed  that 
the  temperature  of  the  water  was  raised  by  2  to  3°F  during 
the  passage  through  the  pump,  presumably  due  to  frictional 
effects.  In  order  to  minimize  the  temperature  difference 
between  the  tank  and  the  ends  of  the  specimen,  the  pumps 
were  rearranged  so  as  to  suck  water  from  the  battery  of 
specimens  and  deliver  it  into  the  tank.  In  this  manner, 
water  could  flow  through  the  end  plugs  of  the  specimens 
under  the  combined  influence  of  the  hydrostatic  pressure 
of  the  body  of  water  in  the  tank  and  the  suction  pressure 
exerted  by  the  pump;  this,  together  with  the  smaller  length 
of  discharge  pipe,  facilitated  an  increased  discharge  while 
eliminating  the  difference  in  temperature  between  the  tank 
and  the  ends  of  the  specimens  (Plate  1(a)) . 
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At  the  cold  end,  the  heat  generated  by  the  passage 
of  water  through  the  system  was  taken  away  by  the 
refrigerator  unit  automatically.  Since  there  was  no  such 
heat-absorbing  unit  at  the  hot  end,  a  small  length  of 
thin  copper  (—  dia.  )  tube  was  interposed  between  the 
discharge  end  of  the  pump  and  the  outlet  in  the  hot  tank; 
this  piece  of  copper  tube  was  kept  immersed  in  the  cold 
tank  so  that  the  water  flowing  through  would  be  cooled 
down;  the  temperature  was  again  raised  to  the  desired 
value  by  the  heating  element  in  the  hot  tank.  The  very 
large  volume  of  water  in  both  tanks  (about  30  gallons  each, 
compared  to  the  rate  of  pumping  which  did  not  exceed  3 
gallons  per  minute)  acted  as  a  buffer  in  minimizing  the 
effect  of  the  quantity  of  heat  on  the  temperature  of  the 
body  of  water.  In  the  final  arrangement,  the  individual 
temperatures  of  water  in  both  tanks  did  not  vary  by  more 
than  -  1 °F . ,  with  the  cold  end  temperature  varying  within 
even  narrower  limits. 

5.2.7.  Thermal  insulation : 

The  two  tanks  were  completely  insulated  by  commercially- 

II 

made,  3  thick  cotton  wool  pads.  The  battery  of  specimens, 
together  with  the  pipe  lines  were  placed  inside  a  specially 
constructed  plywood  box  (approximately  12 I0"  x  2 I0"  x  O'S") 
and  the  box  was  filled  with  "zonolite"  completely,  thus 
providing  adequate  thermal  insulation  for  the  system. 
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5.3.  Experimental  Procedures 
5.3.1.  Compaction  of  specimens: 

With  the  desire  to  achieve  a  highly  uniform  and 
reproducible  degree  of  compaction,  various  conventional 
and  novel  methods  were  tested:  these  methods  included 
(1)  dynamic  compaction  (similar  to  conventional  Proctor's 
method) ,  (2)  static  compression  from  one  end,  (similar 

to  Porter's  method) ,  (3)  static  compression  from  both 

ends,  and  (4)  static  axial  compression  aided  by  lateral 
vibratory  blows  delivered  manually  and  with  a  high 
frequency  electric  hammer.  Finally,  static  compression 
from  both  ends,  aided  by  lateral  vibratory  blows  with  an 
electric  hammer  was  chosen  as  the  technique  to  be  used; 
even  in  this  case,  the  dry  densities  of  the  various 
sections  of  the  specimen  were  found  to  vary  erratically 
over  a  range  of  -  3  pcf  about  an  average  mean  density  of 
95  pcf.  The  random  nature  of  the  variation  of  densities, 
together  with  the  fact  that  the  experimental  determinations 
involved  a  probable  error  of  -  2.5  pcf  in  the  dry  density 
values,  indicated  that  a  better  degree  of  uniformity  could 
not  be  obtained. 

Before  the  start  of  the  experiments,  about  500  gm„ 
of  the  dry  soil  was  mixed  thoroughly  with  enough  quantity 
of  distilled  water  to  give  the  desired  moisture  content 
and  was  placed  in  a  sealed  water-proofed  cardboard  cylinder. 
All  these  cylinders  were  left  in  the  humid  room  for  a  period 
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Plate  3 

COMPARISON  BETWEEN  VARIOUS  COMPACTION  TECHNIQUES 
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dynamic  compaction  with  a  falling  rammer 

STATIC  COMPRESSION 
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of  at  least  72  hours,  at  the  end  of  which  period,  the 
cylinders  were  re-opened,  the  soil  was  mixed  again  and 
the  necessary  quantity  of  the  wet  soil  to  give  the  desired 
dry  density  was  weighed  into  the  assembled  mold.  The 
remaining  wet  soil  was  used  for  determining  the  initial 
moisture  content  by  oven-drying. 

After  compaction,  the  specimen  (inside  the  'Lucite1 
container)  was  removed  from  the  mold,  covered  tightly  with 
aluminum  foil  and  was  left  in  the  humid  room  till  it  was 
ready  to  be  mounted  in  the  apparatus.  In  all,  six  different 
initial  moisture  contents  were  used  with  five  specimens 
each  for  the  first  run  (temperature  gradient  38  to  110°F.) 
and  four  specimens  each  for  the  second  run  (42  to  78°F.) . 
(The  apparatus  was  designed  to  accomodate  twenty  four 
specimens  at  any  one  time) . 

5.3.2.  Assembling  the  specimens  in  position: 

II 

A  3  layer  of  "Zonolite"  was  first  laid  inside  the 
plywood  box  and  the  brass  end  plugs  were  assembled  in  their 
respective  positions.  (The  inlet  and  outlet  tubes  of  the 
brass  plugs  were  connected  to  the  main  pipes  by  means  of 

II 

short  lengths  of  "Tygon"  tubing) .  An  8  wide  plank  was 
laid  on  top  of  the  insulation  and  was  adjusted  to  be  truly 
horizontal.  The  specimens,  with  their  ends  trimmed  plane, 
were  then  mounted  in  position  between  the  end  plugs;  the 
plugs  and  the  specimens  were  bound  together  with  adhesive 
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tape.  Necessary  thermo-couple  leads  were  then  inserted 
into  the  specimen  at  predetermined  locations  and  then 
the  box  was  completely  filled  with  " Zonolite" . 

5.3.3.  Actual  experiment: 

After  assembling  the  apparatus,  circulation  of 
water  through  the  end  plugs  was  started  and  periodic 
temperature  readings  were  taken,  both  of  the  water  in  the 
tanks  and  of  the  different  locations  along  the  various 
specimens . 

Due  to  experimental  difficulties,  the  first  specimen 
of  the  first  run  of  tests  could  not  be  taken  out  for 
moisture  profile  determination  before  an  elapsed  time  of 
77  hours.  In  view  of  the  fact  that  the  tests  were 
concerned  with  changes  in  the  moisture  profiles  at 
different  times,  it  was  decided  to  retain  this  period 
as  the  same  for  all  the  runs.  Subsequent  specimens  of 
each  run  were  taken  out  at  24  hour  intervals. 

As  soon  as  a  specimen  was  taken  out,  it  was 
clamped  between  two  flat  plates  held  together  in  an 
universal  soil  extruder.  (A  bench  vise  could  have  been 
used  in  the  place  of  the  soil  extruder,  had  the  specimen 
been  much  shorter) „  The  individual  rings  were  cut  out 
by  means  of  a  razor  blade  and  a  wire  saw  and  the  weights 
were  recorded.  The  moisture  content  of  the  soil  in  each 
ring  was  obtained  by  oven-drying  the  wet  soil  for  a 
period  of  about  15  hours?  the  weight  and  diameter  of  the 


. 

_  •  ,  (  *  *  i*  • 

■  .  :  '  :  x  '■ 

'  ''  :  i 

'  -  v  :  '  '.;  ■  ■■  '■  .f'  : 

.  v:  ;  •  i.v  i.  »J' 


. 

..  ■'  v  :  ;  .  ■  '  \t:  : :  >  • . :  • 

:  ::  •  ■ 

i 

. 

'  '■  ■  '■ ' '  :  : 


■  v:  '  :  :  .  .v  ■  !  ‘  . '  >  •  '  f.i  '■■  ■ 

:  , 


'•  .  ‘  .  . :  ;  ::  ' 

.  :  i-d  i  f 


96 


rings  were  measured  and  recorded.  All  weight  measurements 
were  carried  out  in  a  semi-automatic  balance  correct  to 
0.0001  gm.,  and  the  diameter  was  measured  with  a  specially 
made  tapered  micrometer,  correct  to  -  0.001  in. 

During  the  course  of  the  experiments,  the 
temperatures  at  various  sections  of  individual  specimens 
were  measured  at  periodic  intervals  using  iron-constanten 
thermo-couples.  With  a  view  to  minimize  the  errors  in 
the  computed  values  of  density  and  void  ratio  of  the 
various  sections  of  the  specimens,  the  number  of  thermo¬ 
couples  was  kept  as  small  as  possible.  In  general,  the 
temperatures  at  homologous  points  of  the  various  specimens 
showed  no  tendency  to  vary  regularly  with  the  distance 
along  the  apparatus;  this  will  be  seen  from  the  typical 
data  sheet  shown  on  page  98. 

5.3.3.  Computations: 

The  individual  dry  densities  of  the  various 
segments  of  the  specimen  were  used  to  compute  the 
corresponding  void  ratios,  and  then,  in  conjunction  with 
the  computed  average  moisture  content  of  the  soil  within 
the  ring,  the  degrees  of  saturation  were  computed.  The 
degree  of  saturation,  rather  than  the  moisture  content, 
was  chosen  as  the  means  of  presentation  of  the  results, 
on  account  of  the  inevitable  variations  in  dry  density 
(and  hence  void  ratio)  across  the  specimen. 


In  some  cases,  especially  in  the  drier  specimens, 
dry  densities  could  not  be  determined  for  all  of  the 
individual  segments,  since  the  dry  soil  tended  to  slough 
off;  in  all  these  cases,  the  moisture  content  of  the  soil 
in  the  warmer  sections  of  the  specimen  was  not  appreciably 
different  from  zero  per  cent;  for  these  segments,  the  void 
ratio  was  taken  as  the  average  of  the  values  determined 
for  the  colder  segments.  Finally,  the  initial  degrees  of 
saturation  for  all  the  specimens  were  computed  on  the 
basis  of  the  overall  average  void  ratio  of  the  specimen. 

A  typical  data  sheet  showing  the  various  observed 
and  computed  quantities  is  shown  on  page  99. 

Lastly,  mention  must  be  made  of  the  fact  that  the 
average  final  degree  of  saturation  of  any  one  specimen 
(obtained  by  averaging  the  individual  degrees  of 
saturation  of  the  various  segments)  was  found  to  differ 
erratically  from  the  initial  degree  of  saturation  of  the 
same  specimen.  In  view  of  the  uncertainties  involved  in 
the  calculation  of  the  initial  degree  of  saturation,  and 
of  the  experimental  errors  in  the  computation  of  the  final 
degree  of  saturation,  such  differences  were  deemed  to  be 
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FIG.  [  6  ]  TAKING  THERMO-COUPLE  READINGS 
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CHAPTER  VI 


DISCUSSION  OF  RESULTS 


6.1.  Preliminary  considerations: 

Only  one  type  of  soil  was  used  in  the  present 
series  of  tests;  the  compacted  density ,  and  hence  the 
void  ratio,  of  the  specimens  were  constant  within  limits 
of  experimental  error.  Although  the  initial  density 
conditions  of  the  various  sections  of  the  specimens  were 
not  determined,  it  is  believed  that  the  density  conditions 
remained  unaltered  during  thermo-osmotic  flow. 

The  data  obtained  from  the  tests  involve  five 
variables  as  follows: 

i)  the  initial  degree  of  saturation, 

ii)  the  final  degree  of  saturation,  or,  the  change  in 
degree  of  saturation  at  various  points  along  the  specimen, 

iii)  the  distance  along  the  specimen, 

iv)  the  elapsed  time,  and 

v)  the  temperature  gradient . 

As  only  two  temperature  gradients  were  employed,  it  is 


possible  to  eliminate  the  temperature  gradient  as  a  variable. 


by  considering  the  two  series  separately.  The  remaining 
four  variables  have  been  combined  in  different  ways  to 
obtain  four  different  sets  of  graphs  for  each  series,  each 
set  emphasizing  one  aspect  of  the  phenomenon.  These  sets 
of  graphs  will  be  taken  up  one  at  a  time  for  a  discussion 
of  their  implications. 

It  is  well  to  point  out  that  the  initial  degrees 
of  saturation  of  the  various  specimens,  originally 
prepared  for  the  same  moisture  contents,  were  not  exactly 
the  same;  the  individual  values  were  found  to  vary  within 
-  1  per  cent  (degree  of  saturation)  from  the  arithmetic 
mean  value  for  the  group. 

6.2.  Change  in  degree  of  saturation  vs.  time: 

The  first  set  of  graphs  (Plates  4  to  15  inclusive) 
show  how  the  degree  of  saturation  of  the  different  segments 
of  the  specimen  varies  with  time;  in  order  to  eliminate 
the  influence  of  minor  differences  in  the  initial  degrees 
of  saturation  of  the  specimens  having  the  same  nominal 
initial  degree  of  saturation,  the  change  in  degree  of 
saturation  has  been  expressed  as  a  percentage  of  the  initial 
degree  of  saturation;  each  curve  represents  the  changes 
taking  place  at  the  noted  distance  from  the  hot  end. 

6.2.1.  For  the  temperature  gradient  of  4.5°F./cm.: 

Graphs  4  to  9  (pp.io8-"3)  are  for  specimens  subjected 
to  the  high  temperature  gradient  of  4.5°F./cm.  In  the  driest 
group  of  specimens,  the  degree  of  saturation  for  the  first 
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4  cm.  section  decreases  to  practically  zero,  and  does  not 
vary  appreciably  with  time.  For  the  next  three  2  cm. 
sections  of  the  specimen,  the  degree  of  saturation 
decreases  with  time.  At  the  10  cm.  section,  the  degree 
of  saturation  reaches  nearly  zero  at  149  hours,  but 
undergoes  a  small  increase  in  the  next  24  hours.  In 
contrast,  the  degrees  of  saturation  of  the  soil  in  the 
last  three  2  cm.  sections  appear  to  fluctuate  considerably 
with  time.  Such  fluctuations  indicate  that  the  water 
accumulated  at  the  cold  end  must  be  flowing  backwards 
under  the  influence  of  some  potential  gradient.  The 
absence  of  any  such  fluctuations  at  the  warmer  sections 
of  the  specimen,  taken  in  conjunction  with  the  observed 
fact  that  these  sections  are  practically  dry,  suggest 
that  the  return  flow  is  in  the  liquid  phase,  and  that 
the  flow  takes  place  up  to  such  distances  behind  the 
cold  end,  as  permitted  by  the  presence  of  continuous 
water  films. 

Ample  evidence  of  the  fluctuations  in  the  degree 
of  saturation  of  the  various  sections  of  the  specimens 
is  also  provided  by  the  remaining  five  graphs  of  the  set. 

It  will  be  noticed  that,  with  increasing  degree  of 
initial  saturation,  the  zone  of  drying  moves  towards 
the  hot  end  progressively,  till  it  disappears  altogether 
with  relatively  wet  specimens  (degree  of  saturation  greater 

* 

The  "zone  of  drying"  refers  to  that  section  of  the  specimen 
where  the  decrease  in  degree  of  saturation  is  nearly  100  per 
cent . 
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than  about  20%:  this  is  illustrated  by  Plates  25-31  also) . 

It  will  .also  be  noticed  that*  as  the  initial  degree  of 
saturation  increases*  larger  and  larger  portions  of  the 
specimen  undergo  fluctuations  in  the  change  in  degree  of 
saturation.  This  is  to  be  expected*  since*  with  increasing 
degree  of  wetness  of  the  specimen*  continuous  conducting 
water  films  are  present  for  longer  distances  back  of  the 
cold  face  and  so  the  return  liquid  flow  extends  to 
greater  and  greater  distances  from  the  cold  end. 

A  final  point  worthy  of  consideration  is  that  the 
accumulation  at  the  cold  end  is  the  greatest  for  the  driest 
specimen*  on  a  percentage  basis.  However*  in  terms  of 
actual  magnitudes*  this  is  really  not  so;  comparing  the 
specimens  of  the  first  and  second  groups  (initial 
saturation  7.8%  and  16.9%  respectively)  *  it  is  easily 
checked  that  an  increase  in  the  degree  of  saturation  of 
about  260%  in  the  former  case  is  numerically  smaller  than 
a  similar  increase  of  about  200%  in  the  latter. 

6.2.2.  For  a  temperature  gradient  of  2.25°F./cm.  (Plates  10-15) 

The  general  trends  in  all  cases  of  this  series  are 
well  defined.  The  same  general  conclusions  deduced  from 
the  first  series  are  equally  applicable  to  this  series* 
thus  supporting  the  validity  of  the  conclusions.  The  main 
difference  between  the  two  sets  of  plots  lies  in  the  extent 
of  drying  zones  of  parallel  specimens  and  the  corollary 
extent  of  the  region  of  recirculation  of  water.  This  difference 
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arises  from  the  fact  that  the  vapour  pressure  gradients 
were  smaller  in  the  latter  case  than  in  the  former;  but 
the  initial  relative  humidity  and  all  other  conditions 
were  nearly  identical.  Hence,  the  rate  of  building  up  of 
a  suction  gradient  from  cold  to  the  hot  side  was  greater 
in  the  latter  case  than  in  the  former,  thus  starting  a 
return  flow  before  appreciable  portions  at  the  warmer  end 
could  become  so  dry  as  to  prevent  active  liquid  flow. 

Thus,  the  data  support  Hutcheon1^  hypothesis  (see  Section 
4.3.4.)  . 

An  inexplicable  situation  presents  itself  in  the 
case  of  the  group  of  specimens  with  an  initial  degree  of 
saturation  of  29.6%  (Plate  14)  .  The  changes  in  the  degrees 
of  saturation  of  the  various  specimens  of  the  group  except 
one  are  very  small,  within  -  15  per  cent  of  the  initial 
value.  The  second  specimen  (taken  out  after  102.5  hours) 
shows  a  reduction  of  99.1  per  cent  at  the  warmest  end  and 
an  equal  accumulation  at  the  14  cm.  section  immediately 
back  of  the  cold  end.  As  the  graphs  in  Plate  13  show, 
the  maximum  net  transfer  in  the  group  with  23.3  per  cent 
saturation  is  only  about  35  per  cent,  and  at  higher  degrees 
of  saturation,  one  should  expect  the  transfer  to  be  smaller 
still.  This  is  indeed  the  case  with  the  last  group  of 
specimens  with  37.1  per  cent  saturation  (Plate  15) .  This 
anomalous  situation  is  further  discussed  in  a  subsequent 


section. 
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6.3.  Change  in  degree  of  saturation  vs.  initial  saturation: 

The  same  data  have  been  plotted  in  a  different 
manner,  wherein  the  per  cent  change  in  the  degree  of 
saturation  has  been  plotted  as  a  function  of  the  initial 
saturation  for  various  distances  along  the  specimen  (Plates 
16-24  inclusive) .  For  the  sake  of  clarity,  the  elapsed 
time  has  been  treated  as  constant  for  any  one  series  of 
graphs  of  the  set,  although  all  the  specimens  of  any 
particular  series  were  not  sampled  simultaneously.  The 
effect  of  initial  saturation  on  the  net  change  is  brought 
out  very  clearly  by  this  series  of  graphs. 

6.3.1.  For  the  higher  temperature  gradient: 

The  most  striking  feature  of  all  the  graphs  (16-20) 
is  that  the  per  cent  change  decreases  rapidly  with  increasing 
degree  of  initial  saturation,  and  that  the  family  of 
curves  in  any  one  group  seems  to  converge  at  a  particular 
value  of  the  initial  saturation.  The  point  of  convergence 
is  not  the  same  for  the  various  time  intervals  but  ranges 
between  30  and  40  per  cent  initial  degree  of  saturation. 

This  is  explainable  as  follows:  the  moisture  profile  at  any 
time  is  governed  by  the  rate  of  vapour  diffusion  from  the 
hot  to  the  cold  side  as  well  as  by  the  rate  of  recirculation 
from  the  cold  to  the  hot  side;  these  rates  are  different 
for  specimens  with  different  degrees  of  initial  saturation, 
as  both  the  rates  are  ultimately  governed  by  the  relative 
humidity  of  the  soil  water,  the  amount  of  air  filled  pores 
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and  the  unsaturated  permeability,  all  of  which  are  in  some 
manner  related  to  the  degree  of  saturation.  Thus,  when 
two  specimens  with  different  degrees  of  saturation  are 
taken  out  after  a  certain  elapsed  time,  vapour  diffusion 
may  be  dominant  in  one  case,  while  recirculation  may  be 
dominant  in  the  other.  The  same  sequence  of  events  is  not 
likely  to  be  maintained  at  a  different  elapsed  time,  except 
probably  when  both  specimens  have  independently  attained  a 
state  of  dynamic  equilibrium.  Hence,  the  point  of 
convergence  cannot  become  stationary  unless  all  specimens 
have  attained  dynamic  equilibria,  when  the  individual  moisture 
profiles  become  independent  of  the  elapsed  time. 

The  preceding  argument  also  serves  to  explain  the 
undulating  nature  of  the  individual  members  of  each  family 
of  curves. 

6.3.2.  For  the  smaller  temperature  gradient: 

These  families  of  curves  (21-24)  need  no  special 
comment,  since  the  general  trends  are  exactly  similar  to 
the  previous  ones. 

6.4.  Moisture  profiles  at  various  times: 

In  plates  25  to  36  inclusive,  the  final  degrees  of 
saturation  at  various  distances  along  the  individual 
specimens  have  been  plotted  as  a  function  of  the  distance 
from  the  hot  end.  In  order  to  bring  out  the  variation  of 
the  moisture  profile  with  time,  the  data  obtained  from  various 
specimens  having  approximately  the  same  initial  degree  of 
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saturation,,  but  sampled  after  different  elapsed  times, 
have  been  grouped  together.  Besides  showing  that  there  is 
an  accumulation  of  water  at  the  cold  end  in  the  case  of 
"dry"  specimens,  these  graphs  also  reveal  that  a  recirculation 
from  the  cold  to  the  hot  side  must  be  taking  place. 

6.4.1.  For  the  higher  temperature  gradient: 

The  general  features  of  these  families  of  curves 
(Plates  25-30)  are  as  would  be  expected  in  the  light  of 
the  previous  discussions.  In  the  "dry"  specimens,  the 
zone  of  drying  progressively  extends  towards  the  cold  end, 
while  the  moisture  profiles  of  the  colder  sections  show 
marked  fluctuations;  this  is  the  case  in  the  first  two 
families  of  curves  (initial  degrees  of  saturation  7.8  and 
16.9  per  cent  respectively) .  In  the  third  family  of 
curves,  the  extent  of  the  zone  of  drying  also  undergoes 
fluctuations,  as  is  evident  by  comparing  the  profiles 
corresponding  to  elapsed  times  of  150  and  174  hours.  With 
further  increase  in  initial  degree  of  saturation,  the 
zone  of  fluctuation  extends  progressively  towards  the  hot 
end,  the  zone  of  drying  retreats  in  the  same  direction 
and  finally  disappears.  A  final  point  observable  is  that 
the  range  of  fluctuation  at  any  section  also  decreases 
with  increasing  degrees  of  initial  saturation. 

6.4.2.  For  the  smaller  temperature  gradient: 

The  general  features  of  these  curves  (31-36)  are 
identical  with  the  ones  described  in  the  preceding  section. 
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The  only  anamoly  is  observed  in  the  family  of  curves 

corresponding  to  the  initial  degree  of  saturation  of 

29.3  per  cent  (Plate  35) .  As  is  seen  from  these  curves, 

the  profiles  after  elapsed  times  of  79,  126  and  150 

hours  are  not  materially  different  from  one  another, 

while  the  profile  after  102.5  hours  shows  a  marked 

drying  at  the  hot  end  and  a  considerably  greater  accumulation 

of  water  at  the  cold  end;  the  intermediate  sections  have 

profiles  similar  to  those  of  the  rest  of  the  specimens  of 

this  group.  As  mentioned  in  an  earlier  section  (cf.  section 

6.2.2),  this  situation  is  inexplicable. 

That  such  a  state  of  instability  does  exist  is 
brought  out  by  comparing  the  graph  of  the  similar  group 
of  specimens  (initial  degree  of  saturation  30.9  per  cent), 
subjected  to  a  temperature  gradient  of  4.5°F./cm.  (see 
Plate  29) .  In  this  case  also,  the  profile  at  an  elapsed 
time  of  103  hours  shows  larger  fluctuations  than  the 
profiles  at  other  times.  As  both  groups  of  specimens 
were  not  tested  over  a  sufficiently  long  period  of  time, 
we  are  not  in  a  position  to  state  definitely  that  such 
an  instability  will  occur  after  an  elapsed  time  of  204 
hours  also;  but  such  a  possibility  cannot  be  excluded 
altogether „ 

If  it  is  established  that  such  a  large  accumulation 
does  take  place  periodically,  it  will  be  impossible  to 
account  for  it  from  the  concept  of  vapour  flow  alone. 
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because,  in  this  range  of  values  of  saturation,  the  air-filled 
porosity  and  the  amount  of  continuous  air  voids  are  both 
small  and  it  is  inconveivable  that  such  a  large  quantity  of 
water  could  have  been  transported  in  the  vapour  phase . 

6.5.  Moisture  gradients  vs.  initial  saturation: 

The  variation  in  the  moisture  concentration  gradient 
between  the  cold  and  hot  ends  of  the  specimen  with  the 
initial  degree  of  saturation,  at  various  elapsed  times  is 
illustrated  by  Plates  37  and  38.  The  existence  of  an 
optimum  initial  degree  of  saturation  for  maximum  gradient 
is  quite  apparent  from  these  graphs. 

Of  necessity,  the  moisture  concentration  gradient 
has  been  assumed  to  be  uniform  across  the  length  of  the 
specimen  and,  accordingly,  the  ordinates  of  these  graphs 
represent  the  difference  in  the  final  degree  of  saturation 
between  the  cold  and  the  hot  ends  of  the  specimen.  It  is 
quite  obvious,  as  the  preceding  sets  of  graphs  have  shown, 
that  the  gradient  may  be  much  steeper  over  some  portions 
of  the  specimen,  than  the  overall  gradient.  However,  as 
the  portions  having  such  steep  gradients  are  not  the  same 
for  different  specimens,  it  is  not  possible  to  compare  the 
individual  gradients  with  gradients  over  homologous 
sections  of  other  specimens. 

6.5.1.  For  the  higher  temperature  gradient: 

The  first  of  these  graphs  (Plate  37)  represents  the 
data  obtained  with  a  temperature  gradient  of  4„5°F./cm„ 
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The  'optimum1  initial  degree  of  saturation  is  found  to 
vary  between  15  and  20  per  cent.  The  influence  of  the 
initial  degree  of  saturation  in  governing  the  net  transfer 
of  water  from  the  hot  to  the  cold  side  is  once  again 
apparent?  in  particular,  the  difference  in  degree  of 
saturation  between  the  ends  is  never  greater  than  10 
per  cent,  for  all  specimens  with  an  initial  degree  of 
saturation  greater  than  about  28  per  cent. 

6.5.2.  For  the  smaller  temperature  gradient: 

In  this  case  (Plate  38) ,  the  optimum  initial 
degree  of  saturation  is  seen  to  range  between  13.5  and 
17.2  per  cent.  The  anamolous  existence  of  two  optimum 
degrees  of  saturation  for  the  specimens  subjected  to  a 
temperature  gradient  for  100-103  hours  has  already  been 
discussed.  The  lower  optimum  value  has  been  used  in  the 
computation  of  the  range  of  variation  mentioned  above. 

6.5.3.  Optimum  degree  of  saturation  vs.  time: 

In  view  of  the  fact  that,  in  both  the  series  of 
tests,  the  optimum  initial  degree  of  saturation  varies 
with  time,  it  is  of  interest  to  ascertain  whether  the 
variation  shows  any  regular  trend  or  not.  For  this 
purpose,  the  optimum  initial  degrees  of  saturation, 
obtained  from  the  above  discussed  curves,  have  been 
plotted  against  time  (Plate  39) .  As  will  be  seen  from 
these  two  curves  (one  each  for  each  temperature  gradient) , 


F  I 

'•  ,  : :  : 

1 

.  .  . 

\  V  /  .  .. 

,  .  . 

’  "... 

. 


- 

.  : . 

, 

*  , 

T 

■  .  ;  .  :  ■  /  ! 

-i  . 


z 

o 

< 

z 

D 

h 

< 

V) 


H 

Z 


D 


CL 

O 


152 


153 


the  variation  in  the  optimum  value  with  time  does  follow 
a  regular  pattern. 

While  limited  data  forbid  us  to  reach  any  definite 
conclusions  of  general  validity,  it  seems  reasonable  to 
hypothesize  that  this  variation  in  the  optimum  initial 
degree  of  saturation  is  also  due  to  the  transient  nature 
of  the  flow?  when  all  the  specimens  have  reached  a  state 
of  dynamic  equilibrium,  this  variation  will  likely 
disappear . 

6.6.  Moisture  gradient  vs.  time: 

The  same  data  have  been  plotted  in  a  different 
manner,  in  order  to  bring  out  the  time -dependent  nature 
of  the  phenomenon:  the  elapsed  time  has  been  chosen  as 
the  independent  variable?  the  data  obtained  from  specimens 
having  nearly  the  same  initial  degree  of  saturation  have 
been  grouped  together.  Both  sets  of  graphs  of  this 
series  (Plates  40  and  41)  show  clearly  that  the  moisture 
concentration  gradient  varies  with  time  and  that  it  attains 
maxima  and  minima?  of  particular  interest  is  the  fact  that 
the  abscissae  of  the  points  of  inflexion  of  the  various 
members  of  the  family  of  curves  are  not  equal.  Thus,  it 
is  seen  that  the  conditions  in  the  various  specimens  are 
different  from  one  another  at  any  one  time?  while  the 
conditions  in  one  specimen  are  such  as  to  give  a  maximum 
moisture  gradient  from  the  cold  to  the  hot  side,  the 
conditions  in  another  specimen,  for  the  same  elapsed  time, 
are  tending  towards  a  minimum  value  of  the  gradient  and  so  on. 
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The  final  point  worthy  of  mention  is  that  the  gradient 
assumes  small  negative  values  in  the  wetter  specimens 
after  certain  elapsed  times,  thus  indicating  an  accumulation 
of  water  at  the  hot  end  also.  As  will  be  seen  from  the 
graphs,  this  accumulation  is  of  a  temporary  nature  only. 

Even  so,  such  an  accumulation  at  the  hot  end  could  not 
have  taken  place,  unless  there  was  a  return  flow  of  water 
towards  the  hot  end. 

6.7.  Concluding  remarks: 

The  data  obtained  from  these  tests  show  very 
clearly  that  water  does  flow  back  from  the  cold  to  the 
hot  side,  the  extent  of  the  flow  region  depending  upon 
the  distance  up  to  which  continuous  water  films  are  present. 

As  regards  the  primary  mechanism  involved  in  the  thermo- 
osmotic  flow,  it  appears  probable  that  it  is  comprised  of 
a  vapour  flow  across  the  empty  pores  followed  by  liquid 
flow  through  the  contact  meniscii.  The  following 
considerations  lend  support  to  this  tentative  conclusion: 

(1)  It  has  been  shown  that  the  temperature  does 

not  appreciably  affect  the  t?F  of  the  soil  water  at  unaltered 
water  contents  (Croney  and  Coleman,  1948) ;  so,  the  flow 
cannot  be  wholly  accounted  for  by  the  thermally-induced 
gradient. 

(2)  Except  in  very  dry  soils,  the  soil  moisture  will 
be  distributed  as  held  water  around  contact  wedges  and  so 
flow  cannot  be  exclusively  in  the  vapour  phase  since 
continuous  air  voids  are  not  present. 
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As  regards  the  existence  of  an  optimum  initial 

degree  of  saturation  for  maximum  thermo-osmotic  flow,  the 

data  do  reveal  that  such  is  indeed  the  case,  in  accordance 

with  the  findings  of  the  earlier  workers  in  the  field. 

But,  there  seems  to  be  no  general  measure  of  agreement 

between  the  reported  order  of  magnitudes  of  this  optimum 

value:  Bouyoucos  reportedly  states  that  the  value  is 

close  to  the  plastic  limit  of  the  soil,  whereas  the  data 

of  Habib  (1957)  suggest  that  the  optimum  value  is  far 

below  the  plastic  limit  of  the  soil.  Winterkorn  (1947, 

1955a)  states  that  the  value  lies  in  the  vicinity  of  the 

Proctor  optimum  moisture  content  for  maximum  densif ication, 

while  Hutcheon  (1958)  reports  that  the  value  lies  between 

50  and  60  per  cent  of  the  3/10  atmosphere  percentage  for 
* 

the  soil. 

For  the  soil  used  in  the  present  series  of  tests,  the 
optimum  initial  degree  of  saturation  is  found  to  range  from 
15  to  20  per  cent  for  the  higher  temperature  gradient  of 
4.5°F./cm.,  and  from  13.5  to  17.2  per  cent  for  the  smaller 
temperature  gradient  of  2.25°F./cm.;  based  on  an  average 
void  ratio  of  0.72  for  all  the  specimens  tested,  it  is 
seen  that  the  air-filled  porosities  range  from  32  to  38 
per  cent  of  the  volume  of  the  soil.  Habib  reports  that 
the  air-filled  porosity  in  his  experiments  was  26%.  There 
are  two  probable  reasons  for  this  discrepancy  between  the 

*  .  , 

The  3/10  atmosphere  percentage  is  approximately  equal 

to  the  field  moisture  equivalent  of  the  soil. 


two  sets  of  data:  in  the  first  place,  the  soils  used  in 
the  two  cases  are  vastly  different:  Habib  used  a  plateau 
loam  of  low  plasticity  (w1  =  34%,  w^  =  18%,  I  =  16%)  , 
while  the  present  experiments  were  conducted  using  a 
nonplastic,  closely -graded  sample  of  fine  silica  sand. 
Secondly,  since  the  present  series  of  experiments  were 
primarily  concerned  with  the  transient  stages  of  thermo- 
osmotic  flow,  the  specimens  had  not  attained  a  state  of 
dynamic  equilibrium,  while  it  is  possible  that  Habib's 
results  are  based  on  data  obtained  from  specimens  after  a 
stable  distribution  of  moisture  was  established.  It  is 
believed  that  the  results  of  the  present  series  of  tests 
would  have  been  lowered  somewhat,  if  the  specimens  had 
been  allowed  to  reach  a  state  of  equilibrium. 

For  the  soil  tested,  the  optimum  moisture  content 
for  maximum  density  determined  by  the  Standard  Proctor  test 
is  4.85  per  cent,  corresponding  to  a  degree  of  saturation 
of  17.3  per  cent.  It  will  be  noticed  that  this  value  is 
within  the  range  of  variation  of  the  optimum  initial  degree 
of  saturation  for  maximum  thermo-osmotic  transfer  for  the 
same  soil.  Thus,  the  experimental  data  are  in  accordance 
with  Winterkorn's  contention  that  the  optimum  moisture 
content  for  maximum  thermo-osmotic  flow  falls  close  to  the 
Proctor  optimum  moisture  content  (Winterkorn  1947) . 


I 
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CHAPTER  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 


7.1.  General  conclusions: 

1.  The  tp  concept  seems  to  be  very  valuable  in 
dealing  with  problems  involving  accumulation  of  water 
beneath  impervious  pavements;  successful  applications 
of  the  concept  have  so  far  been  restricted  to  shallow 
water  table  conditions  and  to  cases  involving  liquid 
flow  only. 

2.  Under  favourable  conditions  of  soil  moisture, 
considerable  quantities  of  water  may  migrate  under  the 
influence  of  a  temperature  gradient;  the  flow  is  from 
the  hot  to  the  cold  side. 

3.  For  a  given  soil,  there  exists  an  optimum 
moisture  content  at  which  thermally-induced  transfer 
of  water  from  the  hot  to  the  cold  side  is  a  maximum. 

4.  The  overall  phenomenon  of  thermo-osmosis  through 
soils  is  still  imperfectly  understood;  many  controversial 
hypotheses  have  been  put  forward  by  the  various  investigators 
to  explain  the  observed  phenomena. 
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7.2.  Conclusions  from  the  present  series  of  experiments: 

1.  The  experimental  data  indicate  that  water  flowing 
from  the  hot  to  the  cold  side  under  a  temperature  gradient 
tends  to  redistribute  itself  in  the  liquid  phase;  the 
distance  up  to  which  such  a  recirculation  takes  place 
depends  upon  the  presence  of  continuous  conducting  water 
films „  The  zone  of  fluctuation  (where  the  thermo -osmotic 
flow  and  the  opposing  liquid  flow  take  place  at  the  same 
time)  increases  with  increasing  initial  degree  of  saturation 
of  the  specimen. 

2.  The  existence  of  an  optimum  initial  moisture 
content  (or  initial  degree  of  saturation)  for  maximum  net 
transfer  of  water  from  the  cold  to  the  hot  side  has  been 
verified.  This  optimum  degree  of  saturation  appears  to 
vary  slightly  with  the  magnitude  of  the  temperature  gradient. 
A  unique  optimum  value  for  any  given  temperature  gradient 
can  be  obtained  only  when  the  individual  specimens  have 
reached  a  state  of  dynamic  equilibrium. 

3.  For  the  soil  tested,  the  optimum  initial  degree  of 

saturation  is  found  to  vary  between  15  and  19.5  per  cent 
for  a  temperature  gradient  of  4.5°F./cm.  (38-110°F .)  ,  and 

between  13 0 5  and  17.2  per  cent  for  a  temperature  gradient 
of  2.25°F./cm.  (42-78°F.) ;  this  variation  may  be  ascribed 

to  the  fact  that  the  specimens  had  not  reached  a  state  of 
dynamic  equilibrium. 
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4.  There  exists  a  particular  value  of  initial  degree 
of  saturation  beyond  which  the  net  transfer  of  water  is 
negligible;  for  the  soil  tested,  this  value  is  about  28 
per  cent. 

5.  In  specimens  having  an  initial  degree  of  saturation 
greater  than  30  per  cent,  the  experiments  reveal  that  a 
temporary  accumulation  of  water  may  take  place  at  the  hot 
end  also;  this  condition  appears  to  be  a  direct  consequence 
of  the  existence  of  a  return  flow  of  water  from  the  cold 

to  the  hot  end . 

6.  The  data  indicate  that  the  mechanism  of  thermo¬ 
osmosis  involves  vapour  flow  across  air-filled  pores  and 
liquid  flow  through  the  contact  wedges  of  held  water. 

Upon  the  establishment  of  a  suction  gradient  of  sufficient 
magnitude,  the  accumulated  water  at  the  cold  end  tends  to 
redistribute  itself  by  flowing  towards  the  hot  end,  as 
far  as  continuous  water  films  are  present.  In  very  moist 
specimens,  a  state  of  dynamic  equilibrium  may  never  be 
reached.  Thus,  the  experimental  data  are  in  accordance 
with  the  findings  of  earlier  workers  like  Gurr  et  al., 

Taylor  and  Cavazza,  and  Hutcheon. 

7.3.  Recommendations  for  further  investigations: 

1.  Further  experimental  work  in  thermo-osmosis  may 
be  so  designed  as  to  measure  the  of  the  soil  water  at 
various  sections  of  the  specimens  as  a  function  of  time 
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and  of  the  initial  degree  of  saturation.  Such  experiments 

may  prove  to  be  invaluable  in  many  ways:  they  may  throw 

more  light  on  the  actual  mechanism  involved  in  thermo- 

osmotic  flow;  and,  they  may  also  pave  the  way  for  the  eventual 

application  of  the  P  concept  to  the  case  of  moisture 

F 

migration  in  the  vapour  phase. 

Experience  shows  that,  with  fairly  high  temperature 
gradients,  the  soil  close  to  the  warm  end  becomes  very  dry, 
especially  in  specimens  having  low  initial  moisture  contents. 
The  I?  of  soil  water  in  such  sections  is  likely  to  be  very 
high,  depending  upon  the  soil  type.  No  suitable  techniques 
for  in-place  measurement  of  F  beyond  3  are  available.  So, 
new  techniques  may  have  to  be  developed,  or  high  temperature 
gradients  must  be  avoided. 

2.  It  is  also  desirable  to  conduct  these  experiments 
using  different  types  of  soils,  and  ascertain  whether  there 
is  any  correlation  between  the  factors  involved  in  thermo¬ 
osmosis  on  the  one  hand,  and  the  properties  of  the  soils 
on  the  other  hand.  The  former  may  include  the  following: 
maximum  thermo-osmotic  transfer,  thermo -osmotic  permeability, 
optimum  initial  degree  of  saturation  for  maximum  flow,  and 
the  minimum  initial  degree  of  saturation  for  negligible  flow; 
the  latter  may  include:  limits  of  consistency,  optimum 
moisture  contents  for  maximum  densif ication  according  to 
the  different  specified  compaction  procedures,  base  exchange 
capacity,  hygroscopic  moisture  content  and  type  of  soil 
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3.  The  present  series  of  experiments  were  conducted 
with  closed  systems.  In  the  field,  the  subgrade  is 
usually  in  communication  with  a  water  table  below;  if 
the  surfacing  is  freely  permeable,  the  field  conditions 
will  approximate  to  an  open  system.  The  case  of  a  subgrade 
overlain  by  an  impervious  pavement  belongs  to  an  intermediate 
state  of  conditions,  if  we  assume  that  lateral  movement  of 
water  from  the  center  of  the  pavement  towards  the  shoulders 
is  negligible.  Thus,  it  will  be  of  great  practical  interest 
to  investigate  the  extent  of  accumulation  at  the  cold  end 

in  a  laboratory  model  simulating  these  field  conditions 0  A 
simple  testing  arrangement  for  this  purpose  would  seem  to 
be  to  replace  the  brass  end  plug  at  the  hot  end  by  a  plug 
made  of  porous  bronze  or  similar  material;  piezometer  tips 
located  close  to  the  hot  end  may  serve  to  indicate  the 
pressure  of  water  at  entry,  which  may  be  utilized  in  making 
allowance  for  the  flow  under  hydrostatic  pressure  alone. 

4.  In  view  of  the  fact  that  it  is  extremely  difficult, 
if  not  altogether  impossible,  to  prepare  two  identical 
specimens  of  the  same  soil,  it  is  desirable  that  a  suitable 
technique  for  in-place  measurements  of  moisture  contents 

of  laboratory  soil  samples  be  evolved „  Such  a  technique 
will  eliminate  the  need  for  a  large  number  of  specimens 
having  the  same  initial  moisture  content  and  the  apparatus 
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can  be  simplified  considerably.  More  important,  however, 

is  the  fact  that  it  will  be  possible  to  dispense  with 

the  need  of  removing  the  specimen  from  the  apparatus  and 

of  determining  the  moisture  content  by  the  conventional 

method  of  oven  drying.  By  removing  the  specimen  from 

the  apparatus ,  we  alter  the  temperature  conditions  across 

the  specimen  thus  upsetting  the  vapour  equilibrium 

conditions  at  the  different  sections  of  the  specimen; 

thus  liquid  water  may  pass  into  the  vapour  phase  at  the 

cold  end,  while  the  water  vapour  may  condense  into  the 

liquid  phase  at  the  warm  end.  Since  the  relative  humidity 

and  the  temperature  conditions  are  different  for  different 

sections  of  the  same  specimen,  this  difficulty  cannot  be  overcome 

by  employing  a  constant  temperature  room  or  room  with 

controlled  humidity  conditions  for  the  location  of  the 

apparatus.  The  need  for  an  acceptable  method  of  in-place 

measurement  of  the  moisture  content  cannot  possibly  be 

overemphasized . 
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APPENDIX  A 


HYDRODYNAMICS  OF  SOIL  MOISTURE 

A.  1.  Basic  equations  of  classical  hydrodynamics: 

Every  system  involving  flow  of  fluids  must  obey  three  basic 
equations;  the  equation  of  continuity,  the  equation  of  state,  and  the 
Newtonian  equation  of  motion.  The  equation  of  continuity  is  the  same 
as  the  well-known  principle  of  conservation  of  matter,  restated  in 
algebraic  terms.  The  equation  of  state  characterises  the  fluid  under 
consideration  and  the  thermodynamic  nature  of  flow:  if  fluid  is  a  liquid, 
the  liquid  may  be  incompressible  or  compressible;  if  the  fluid  is  a  gas, 
the  flow  may  be  either  adiabatic  or  isothermal.  An  algebraic  statement 
of  the  particular  nature  of  the  flow  and  of  the  fluid  constitutes  the  equation 
of  state.  The  equation  of  motion  specifies  the  nature  of  interaction  be¬ 
tween  the  fluid  and  the  flow  region  or  medium. 

The  law  of  conservation  of  matter  states  that  the  fluid  mass 
in  any  closed  system  can  be  neither  created  nor  destroyed.  For  anal¬ 
ytical  purposes,  this  may  be  restated  as  follows:  'The  net  excess  of 
mass  flux,  per  unit  of  time,  into  or  out  of  any  infinitesimal  volume 
element  in  the  fluid  system  is  exactly  equal  to  the  change  per  unit  of 
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time  of  the  fluid  density  in  that  element,  multiplied  by  the  free  volume 
of  that  element'  (Muskat  1937).  This  statement  takes  the  form  of  the 
equation 


di> 


V  (M)  -  +  ^  - 


2>y 


n 


3P 

3b 


(1) 


where  v  is  the  velocity  vector  at  any  point  with  components  v  ,  v  , 

x  y 

v  ,  in  the  three  co-ordinate  directions ,  P  is  the  fluid  density,  and  n 
z 

is  the  porosity  of  the  medium.  If  the  motion  is  steady  with  respect  to 
time,  the  right  side  of  equation  (1)  vanishes  and  the  equation  reduces 


to 


div  (v)  =  0  (2) 

The  equation  of  state  will^,  in  general,  be  a  quantitative 

statement  of  the  relationship  between  the  density  p,  the  pressure  p, 
and  the  absolute  temperature  T  of  the  fluid,  in  the  form  of  an  equation: 


0  (P, P,T)  =  0 


(3) 


For  example,  if  one  is  dealing  with  strictly  incompressible  liquids, 
the  statement  that  the  density  is  constant  is  an  equation  of  state.  In 
the  case  of  ideal  gases,  the  gas  law,  pv  =  RT  is  an  equation  of  state. 

In  addition  to  the  equation  characterising  the  physical  nature 
of  the  fluid,  v^e  must  also  have  an  equation  which  will  characterise  the 
thermodynamic  nature  of  the  flow  system,  so  that  one  of  the  variables 
may  be  eliminated.  For  homogeneous  fluids  in  viscous  flow,  Muskat 
gives  the  following  equation:  ^ 


e  -  e0  c  £  ) 


p  Pf  (  P  ~  1 


(4) 


1. 


The  equation  given  by  Muskat  is  of  a  slightly  different  form;  the 


form  used  here  is  taken  from  Scheidegger1  s  monograph  (19  57). 
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where  P0is  the  density  of  the  fluid  under  a  pressure  Pa 
p  is  the  density  under  a  pressure  p 
m  &  fif.  are  constants . 

The  particular  set  of  conditions  may  be  signified  as  follows: 


Liquids:  m  =  o 

Incompressible  liquids .  -  o 

Compressible  liquids . ^  o 


Gases:  m  <to  ,  Pf  -  o 

Isothermal  process . m  =  1 

Adiabatic  process  . m  -  c  /  c 

v  p 

where  cv  and  c^  are  the  spectific  heats  of  the  gas  under  constant 
volume  and  constant  pressure  respectively. 

For  incompressible  liquids,  equation  (4)  reduces  to 
and  substitution  in  the  equation  of  continuity  yields 

div  (  v  )  =  o  (5) 

which  now  involves  the  velocity  components  only. 

Having  formulated  the  fundamental  equation  of  continuity  and 
of  state,  it  now  remains  to  characterise  the  flow  dynamically,  and 
state  explicitly  how  the  fluid  reacts  to  press  re  gradients  and  external 
forces.  This  involves  the  statement  of  the  Newtonian  law  of  motion  in 
hydrodynamical  terms.  Newton's  law  states  that  the  force  acting  on 
any  body  is  equal  to  the  mass  of  the  body  times  the  acceleration  of 
the  body  in  the  direction  of  the  force.  Considering  the  case  of  flow  of 
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fluids,  the  external  forces  acting  on  the  mass  of  fluid  of  unit  volume 
may  be  separated  into  three  parts:  the  pressure  gradients  having 
components  in  the  three  co-ordinate  directions;  the  'body  forces'  due 
to  gravity  and  other  causes,  having  components  F  ,  F  ,  F  ;  the 

X  y  z 

viscous  retarding  forces  due  to  internal  resistance  experienced  by 
the  fluid. 

In  text-books  on  hydro  dr  ynamics,  it  has  been  shown  that 
the  final  equations  relating  all  these  forces  to  the  mass  and  acceleration 
of  the  body  may  be  written  as 
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defined  previously, 
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and  \x 

is  the  co-efficient  of  viscosity  of  the  fluid. 

Equations  (6)  are  known  as  the  Navier -Stoker s '  equations, 
which  are  the  hydrodynamical  equivalent  of  the  Newtonian  law  of 
motion.  (For  a  rigorous  derivation  of  these  equations,  see  Lamb  1879, 
Vallentine  1959).  The  terms  x,  ^  y,  and^z  represent  the  accel¬ 
eration  of  the  unit  volume  of  the  fluid  in  the  three  co-ordinate  directions 
the  velocity  of  the  fluid  element  changes  during  any  interval  of  time 
at  the  position  it  occupied  originally;  and,  during  this  interval,  the 
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element  itself  has  moved  to  a  new  position,  which  involves  an 
additional  change  in  velocity,  thus  giving  rise  to  the  total  time  der¬ 
ivative  of  the  velocity.  The  first  two  terms  on  the  right  side  of  the 
equations  need  no  explanation.  The  third  term  is  an  expression  for 
the  viscous  resistance  arising  from  the  fact  that  the  three  velocity 
components  have  gradients  in  all  the  three  directions  at  any  instant; 
the  fourth  term  is  also  a  part  of  the  expression  for  the  viscous 
resistance,  and  has  to  be  included  because  the  gradients  of  the  vel¬ 
ocity  components  of  any  element  depends  on  time  as  well  as  the 
position  of  the  element  in  space. 

Thus,  we  have  five  equations  -  one  equation  of  continuity, 
one  equation  of  state  and  three  equations  of  motion  -  and  five  unknown 
variables  -  p,  P,vx,  Vy,  vz.  (The  variable,  temperature,  can  be 
eliminated  between  the  equation  of  state  and  the  thermodynamic  equation 
characterising  the  flow.  )  These  five  equations  should  be,  and  indeed 
are,  sufficient  to  solve  for  the  five  unknowns.  However,  the  mathem¬ 
atical  difficulties,  arising  from  the  fact  that  the  Navier -Stokes  equations 
are  not  linear,  exclude  the  possibility  of  obtaining  a  general  solution 
for  all  but  the  simplest  of  problems  with  relatively  simple  geometry 
of  the  flow  region  (Muskat  1937,  Vallentine  1959).  Thus,  for  a  straight 
capillary  tube  of  circular  cross-section,  the  Navier -Stokes  equations 

can  be  solved  exactly,  yielding  the  well-known  Poiseuille's  law: 
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where,  'Q'  is  the  volume  rate  of  flow  through  the  capillary, 

'aP1  the  Pressure  drop  across  the  capillary, 

'h1  is  the  length  and  'a'  is  the  radius  of  the  capillary. 

A  few  approximate  solutions  of  the  Navier -Stokes  equations 

have  been  obtained  for  specific  problems:  Thus,  Happel  and  Byrne 

(1954)  obtained  an  expression  for  the  velocity  field  produced  by  the 

presence  of  a  sphere  at  the  axis  of  a  circular  tube  through  which  the 

fluid  is  flowing.  The  expression  was  based  on  the  classical  equations 

of  hydrodynamics,  but  the  inertia  term  -  was  dropped  out. 

Di- 

Using  this  expression,  the  authors  derived  a  simple  equation  for 

the  pressure  drop  across  dilute  systems  of  spherical  assemblages, 

assuming  the  force  exerted  on  each  sphere  to  be  the  same  as  the  force 

that  would  be  exerted  if  the  sphere  were  located  at  the  axis  of  the  tube. 

Their  final  equation  is  of  the  form: 

Ap  =  9  ( \~  e)  ■  -BlhJr  ,  (8) 

rs  cA 


or,  transposing, 

a^  P  s 

v  =  _ _ _ 

9  p-  (1  -  e)  L 


(9) 


In  these  equations,  is  the  pressure  drop  across  the  assemblage 

of  spheres, 

v  is  the  macroscopic  velocity  in  the  tube  approach¬ 
ing  the  assemblages  of  spheres, 
e  is  the  fractional  void  volume, 

L  is  the  length  of  the  assemblage  in  the  direction 


of  flow, 
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and  a  is  the  radius  of  the  spheres. 

It  may  be  noticed  that  in  the  case  of  incompressible  liquids, 
the  average  velocity  through  the  tube  must  be  the  same  as  the  average 
macroscopic  velocity  through  the  assemblage,  from  considerations  of 
the  continuity  of  flow. 

Unfortunately,  no  experimental  evidence  has  been  adduced 
by  the  authors  to  support  the  exactness  of  their  final  equation,  and, 
as  such,  the  validity  of  their  assumptions  is  not  known.  The  qualific¬ 
ation  of  the  spherical  assemblage  -  as  dilute  systems  -  is  not  clearly 
defined.  Presumably,  this  refers  to  the  case  where  the  total  volume 
occupied  by  the  spheres  is  negligible  in  comparison  with  the  volume 
of  the  fluid  involved  in  the  flow;  if  this  were  so,  the  usefulness  of 
this  analysis  is  severely  restricted  as  far  as  practical  problems  are 
concerned,  since  the  conditions  obtaining  in  the  field  are,  in  most 
cases,  just  the  reverse  of  what  has  apparently  been  assumed  in  the 
analysis.  And,  finally,  the  analysis  is  only  for  the  case  of  a  system 
of  spheres  of  equal  radius,  and  the  limited  extent  of  applicability  of 
the  analysis  to  soil-water  systems  is  quite  obvious. 

From  the  foregoing,  it  will  be  obvious  that  an  analysis  of 
the  problem  from  the  standpoint  of  classical  hydrodynamics  is  beset 
with  many  mathematical  difficulties  the  circumvention  of  which  in¬ 
volves  a  large  number  of  questionable  assumptions.  Hence,  it  is  not 
surprising  that  attempts  should  have  been  made  to  solve  the  problem 


from  a  semi-empirical  approach.  In  this  semi-empiracal  approach, 
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the  first  two  equations  of  classical  hydrodynamics  -  namely,  the 
equation  of  continuity  and  the  equation  of  state  -  have  to  be  retained, 
since  these  two  equations  represent  the  fundamental  laws  that  must 
be  obeyed  by  all  systems  regardless  of  the  nature  of  the  approach. 
So,  the  difference  between  the  classical  and  the  semi-empirical 
approaches  should  lie  in  the  formulation  of  the  equation  of  motion. 
While  the  classical  approach  adopts  the  Navier -Stoke s  equation,  the 
semi-empirical  approach  adopts  the  experimental  equation  given  by 
Darcy  as  the  starting  point.* 


A.Z.  Semi-empirical  approach  based  on  Darcy's  law: 

Since  Darcy's  law  was  obtained  from  experiments  with  uni¬ 
directional  flow  of  water  through  filter  beds,  a  generalized  version 
of  the  law,  as  applicable  to  three  dimensional  flow  systems  involving 
any  type  of  fluid  has  to  be  formulated.  For  this  purpose,  the  fun¬ 
damental  equation  of  Darcy  is  taken  in  its  alternative  form,  namely, 

V  =  A  .5*5  (10) 

p  ax 

where,  besides  other  quantities  as  already  defined,  K  is  the  so-called 

Z 


co -efficient  of  permeability  with  the  units  of  (distance).  In  this  form, 
the  co-efficient  of  permeability  is  a  property  of  the  porous  medium 
only,  and  is  independent  of  the  properties  of  the  percolating  fluid. 

(cf.  Chap.  II  of  this  dissertation.  ) 

*  The  original  authorship  of  this  semi-empirical  approach 
based  on  Darcy's  law  is  not  clearly  known;  the  analysis  to  be  presented 
in  the  following  pages  is  largely  based  on  the  work  of  Muskat  (cf.  Muskat 
1937). 
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As  a  first  step,  it  is  assumed  that  'the  resultant  velocity 
at  any  point  may  be  resolved  into  three  components,  parallel  to  the 
three  co-ordinate  axes,  each  reacting  to  the  pressure  gradients  in¬ 
dependently  of  the  others'  (Muskat  1937).  Then  Darcy's  law  may  be 
rewritten  as  follows: 

v*=  Kx  ^  H 


v  =  Ky  /  p  . 
y  Dy 


Vz  =  Kz  /  h  .  _g£ 
dz 


(ID 


Equations  (11)  imply  that  the  medium  is  strictly  anisotropic.  However, 
solutions  of  practical  problems  are  always  obtained  on  the  assumption 
that  the  medium  is  isotropic  and  the  final  answers  are  modified  suit¬ 
ably  to  account  for  the  anisotropy  of  the  medium.*  Hence,  the  medium 
will  be  treated  as  isotropic  in  the  following  analysis,  on  the  under¬ 
standing  that  a  transformed  geometry  of  the  porous  medium  is  to  be 
used  in  the  case  of  anisotropic  porous  media. 

It  is  evident  that  when  the  velocity  of  flow  has  a  component 
in  the  vertical  direction,  the  effect  of  the  acceleration  due  to  gravity 
will  have  to  be  taken  into  account.  In  addition,  if  there  are  adsorptive 
force  fields  present,  their  influence  on  the  flow  should  also  be  taken 
into  account. 


*As  an  example,  the  conventional  method  of  sketching  flownets 
may  be  cited.  When  the  medium  through  which  flow  is  taking  place  is 
anisotropic,  the  procedure  involves  the  modification  of  the  scale  of  the 
diagram  so  as  to  obtain  a  transformed  section  that  is  isotropic,  sketching 

the  flow  net  on  this  transformed  section,  and  then  transforming  the  flow 
net  back  to  the  original  scale. 
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These  forces  are  collectively  termed  as  body  forces,  * 
and,  in  the  present  case,  they  include  all  forces  except  the  ones  due 
to  fluid  pressure  difference.  If  we  assume  that  the  body  forces  have 
components  F  F  ,  and  F  in  the  three  co-ordinate  directions, 

y  Z  7 

equations  (11)  may  be  modified  as: 


v„  = 


=  K  /  VL  •  (  -  F  ) 

3x 


Y 


v 


=  K/»  ■  -Fy) 

=  KI*  •  lH-Fz> 


(12) 


If,  further,  we  assume  that  the  body  forces  can  be  repres 
ented  by  a  potential  V  ,  such  that 


Fx  +  Fy  +  F 


-  -  (  QV  -f  3  V  -j-  3  V  \ 

'  -Sv.  rs-  f 


ax  2>y 


;i3) 


3z 


equations  (12)  may  be  written  as 

VX  =  K  /  v  •  |r  (p  +  v) 
Vy  =  K  /  VI  (P  +  V) 


V  =  K  /  p  -2-  (p  +  V) 

^  oz 


Setting  0  =  K  /  p.  .  (p  +  V), 
equations  (14)  become 


(14) 


(15) 


*  It  may  be  mentioned  in  passing  that  this  grouping  of  so 
many  forces  into  one  single  factor  makes  it  very  difficult  to  take  explicit 
account  of  the  individual  components  of  the  body  forces  and  of  the 
effects  physical  changes  on  those  components. 
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(16) 


or  v  =  Vcj)  (17) 

Muskat  calls  '0'  as  the  velocity  potential,  since  the  grad¬ 
ient  of  (j  in  any  one  of  the  co-ordinate  directions  gives  the  velocity 
of  the  fluid  in  that  direction.  Equations  (16),  or  its  vector  equivalent 
(equation  (17))  is  the  final  form  of  generalised  Darcy’s  law  and  may 
be  considered  as  the  macroscopic  equivalent  for  the  Navier -Stokes 
equations  of  classical  hydrodynamics.  (Muskat  1937) 

Now  that  the  dynamical  basis  for  the  analysis  of  problems 
involving  viscous  flow  of  fluids  through  porous  media  have  been 
formulated  in  the  form  of  eq.  (16),  it  is  pertinent  to  examine  the 
assumptions  that  have  been  made  in  the  analysis.  Basically,  it  is 
assumed  that  Darcy’s  law  is  valid  for  the  case  under  consideration. 

In  an  earlier  section  of  this  report,  the  limitations  of  the  law  have 
been  enumerated  (cf.  Chap.  II).  The  generalised  form  of  the  law  is 
subject  to  the  same  limitations.  It  has  been  assumed  that  the  three 
component  velocities  of  any  element  of  fluid  react  to  the  respective 
pressure  gradients  in  the  three  co-ordinate  directions  in  exactly  the 
same  way  as  postulated  in  the  fundamental  Darcy's  law.  The  use  of 
a  velocity  potential  function  presupposes  that  the  body  forces  can  be 
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expressed  in  the  form  of  equation  (13).  In  the  case  of  gravity,  the 
vector  of  the  gravitational  force  field  is  always  directed  vertically 
downwards  and  its  magnitude  is  solely  governed  by  the  mass  of  the 
element  considered  and  its  position  in  the  force  field  with  respec  to 
some  arbitrarily  chosen  datum.  The  magnitudes  and  directions  of 
the  vectors  of  the  remaining  components  of  the  body  forces  -  such  as 
those  due  to  chemical  and  electrical  potentials  -  are  dependent  upon 
the  specific  conditions  of  the  system  and  cannot  be  uniquely  defined 
once  and  for  all. 

After  the  dynamical  equation  of  motion  has  been  formulated 
(Eq.  17),  it  is  only  necessary  to  combine  it  with  the  equations  of 
continuity  and  of  state  to  obtain  the  final  differential  equation  des¬ 
cribing  the  flow.  The  continuity  equation  is 

div*  (  Pv  )  =  -n-|f  (Eq.  i ) 

and  the  generalised  form  of  the  equation  of  state  is 

P  »  P0  (£.  V"  O’-PO  ,  (Eq.  4) 

r  o 

For  steady  flow  conditions,  the  right  side  of  equation  (1)  vanishes; 
for  liquids,  m  =  O  in  equation  (4),  and  if  the  liquid  is  incompressible, 

=  O,  thus  reducing  equation  (4)  to 

P  =  Po  (18) 

Substituting  this  in  Eq.  (1),  we  obtain 


div  (v)  =  O 


(19) 
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and,  substituting  for  v  from  equation  (17), 
div  (vqf>)  =  O 

or,  =  O  (2  0) 

2 

In  equation  (20),  v  is  the  Laplacian  operator,  given  by 
3^  +  3^  +  q2 

9x2  3y2  dz2  .  Similar  expressions  can  be  set  up  for  the  flow 

of  compressible  liquids  and  gases.  In  the  case  of  incompressible 
liquids,  the  time -dependent  quantities  disappear,  and  as  such,  no  non¬ 
steady  states  within  the  system  can  be  present,  unless  the  boundary 
conditions  themselves  vary  with  time.  In  contrast,  the  equations 
describing  the  flow  of  compressible  liquids  and  gases  involve  the 
time-derivative  and  so  their  solutions  are  much  more  complicated. 
However,  the  time  derivatives  do  vanish  under  steady  flow  conditions 
in  this  case  also,  and  the  final  equations  are  similar  to  the  Laplacian 
equation  (20),  the  difference  being  only  in  the  form  of  the  dependent 

r  ■* 

variable.  For  practical  purposes,  it  whll  always  suffice  to  treat 
steady  state  liquid-flow  systems  as  problems  involving  flow  of  in¬ 
compressible  liquids  alone,  in  view  of  the  very  low  compressibilities 
of  normal  liquids  (Muskat  1937). 

A.  3.  Methods  of  solution: 

The  solution  of  equation  (20)  can  be  obtained  in  a  number 
of  ways:  If  the  geometry  of  the  flow  region  is  simple  enough  to  be 
expressed  analytically,  a  purely  analytical  solution  is  possible, 
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employing  the  principles  of  conformal  mapping,  Green's  functions  or 
the  method  of  successive  relaxation  (for  details,  see  Muskat  1937, 

Klute  1952,  Scheidegger  1957,  Schmid  1958).  Such  solutions  require 
a  very  high  degree  of  mathematical  skill,  and,  in  view  of  the  un¬ 
certainties  involved  in  the  description  of  the  porous  medium,  it  is 
doubtful  whether  such  mathematical  precision  is  warranted  or  not. 

A  simpler  method  of  solution  is  suggested  by  the  fact  that 
Eq.  (20)  is  similar  to  many  equations  describing  physical  processes 
such  as  flow  of  heat  and  of  electricity;  solutions  depending  on  the 
analogy  between  the  flow  of  water  and  that  of  heat  or  of  electricity 
are  called  analog  solutions;  these  methods  involve  the  experimental 
determination  of  such  physical  quantities  as  magnetic  flux,  electric 
current,  electric  potential,  or  temperature,  which,  by  their  analog 
interpretation,  give  the  velocity  and  force  potential  of  the  seepage 
flow  problem.  On  account  of  the  high  degree  of  accuracy  obtainable, 
electric  analogs  are  widely  employed  in  preference  to  others. 

A  third  method  of  solution,  suggested  by  Forcheimer,  in¬ 
volves  the  solution  of  the  Laplacian  equation  by  a  graphical  procedure 
of  trial  and  error  (see  A.  Casagrande  1937;  also  Taylor  1958).  This 
method  is  based  on  the  fact  that  the  solution  of  the  Laplacian  equation 
yields  two  families  of  orthogonal  curves;  the  'flow  net'  is  sketched 
with  due  attention  to  the  boundary  and  initial  conditions  of  the  problem 
and  to  the  requirement  that  each  'field'  of  the  flow  net  must  approximate 
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a  square.  This  method  gives  reliable  results  in  a  mere  fraction 
of  the  time  required  for  analytical  solutions,  but  it  cannot  be  used 
for  solving  three-dimensional  flow  problems,  unless  there  is  radial 
symmetry  (Schmid  1958). 

Finally,  solutions  may  be  obtained  from,  mechanical  models 
of  the  flow  problem;  however,  this  method  is  very  expensive,  and 
on  account  of  the  experimental  errors  involved,  the  results  obtained 
therefrom  are  not  fully  reliable. 

Flow  systems  involving  steady  state  flow  of  liquids  with 
a  free  boundary  -  as,  for  example,  seepage  through  dams  -  are  very 
difficult  to  solve  analytically.  Analog  solutions  are  particularly 
useful  in  such  cases.  Graphical  solutions  by  the  Forcheimer  method 
requires  that  the  location  of  the  free  boundary  be  guessed  in  the  first 
place;  a  few  empirical  methods  are  available  for  this  purpose 
(see  Casagrande  1937,  Taylor  1958).  Problems  involving  flow  of 
gases  through  porous  media  are  outside  the  scope  of  this  disseration; 
for  information  on  these  problems,  the  reader  is  referred  to  the  books 
of  Muskat,  Scheidegger,  and  Carman  (1956). 


j-  .  . 

•  t  •  n  .  ...  ;  ..  j, 

<  •  .  -------  ...  .  -  ....  j.  .  . 

/ 

.  .........  x .  V. 


••  •  •  1  •.  .  •  •  •  -  -  <  .  .  :  j  ,  ..  .  ; 

•  -  '  V  •  -  ;  ^  it  .  ■  -  .  ll:  .  a  "i 


-  . 

.....'  .........  . ■ 


.  - 


;  ■  •  -  ~ 

"...  _  •  .  .  :.x,. 


t  — —  ... 

T  1  . . 


A  ■  A  -  -  •  •  -  «*  •  -  xKi  ....  .  . 


. 


'  .  ti  01  .  ......... 

.  •  J  — .  .  ....  :.  ...’  .  .  '..  .  .  :  ; 

.  '  -  •  .  . .  t  -  -  .  } 

:  -  •  • . -  d  ....  t  ...  .  . 

•  •  x  -  ••  ■  -  -  •  :  ■  -  •  ■  -  -  -  .  1-.  ; 

•  >  '  -  ‘  •  •  t  -  ....  { .  .  . 


196 


APPENDIX  B 

CONCERNING  DARCY'S  COEFFICIENT  OF  PERMEABILITY 

B.l.  Experimental  determination  of  the  saturated  permeability: 

B.1.1.  Various  methods  and  optimum  range  of  application; 

Owing  to  the  importance  of  Darcy's  law  in  most  of  the  problems 
involving  flow  of  fluids  through  porous  media,  it  is  only  natural  that 
several  methods  should  have  been  devised  for  the  determination  of  the 
Darcy's  coefficient  of  permeability  of  the  medium.  There  are,  at 
present,  six  different  methods  for  this  determination  in  the  case  of 
soils,  and  each  particular  method  gives  good  results  over  a  particular 
range  of  permeability  values.  The  various  methods  and  their  generally 
accepted  optimum  ranges  of  application  are  listed  below: 

Range  in  terms  of 

Method  permeability  in  cm/sec. 


a. 

Constant  Head  Permeability  test 

1(? 

to 

10"3 

b. 

Variable  Head  Permeability  test 

10~3 

to 

10“5 

c. 

Capillary  Permeability  test 

HT1 

to 

KT5 

d. 

Estimation  from  grain  size  analysis  data 

1C? 

to 

10“3 

e. 

Computation  from  consolidation  test  data 

1CT7 

to 

10~9 

f. 

Field  tests  including  pumping  test 
and  ground  velocity  tests 

102 

to 

10”3 

B.1.2  Critical  evaluation  of  the  merits  of  the  various  tests: 

The  constant  head  permeability  test  is  probably  the  simplest 
of  all  tests  and  the  test  arrangement  is  very  similar  to  the  one  em¬ 
ployed  by  Darcy  in  his  experiments.  The  conditions  of  steady-state 
percolation  are  definitely  established  and  thus  there  can  be  no  question 
of  the  applicability  of  Darcy's  law  for  the  evaluation  of  the  data  from 
the  test.  Some  errors  are  likely  to  result  from  the  formation  of  filter 
skins  of  fine  material  at  the  entrance  and  exit  ends  of  the  specimen. 

In  one  type  of  apparatus,  this  has  been  eliminated  by  providing  pressure 
tappings  within  the  specimen  and  beyond  the  influence  of  the  filter  skins; 
the  hydraulic  head  causing  flow  is  then  taken  as  the  pressure  difference 
between  two  such  tappings,  active  over  that  portion  of  the  specimen  be¬ 
tween  the  two  tappings.  With  relatively  impermeable  soils,  loss  of  out¬ 
flow  water  due  to  evaporation  may  introduce  considerable  errors. 

The  variable  head  permeability  test  is  by  far  the  most  versatile 
of  all  tests,  giving  good  results  over  a  large  range  of  permeability 
values.  But  the  head  causing  flow  through  the  soil  is  varying  continu¬ 
ously  with  time,  and  thus  the  flow  is  not  under  steady-state  conditions. 
The  use  of  Darcy's  law  (which  is  valid  for  steady-state  percolation 
only)  in  the  evaluation  of  the  test  data  may  seem  to  be  inherently  erron¬ 
eous.  However,  it  should  be  pointed  out  that,  in  the  analysis,  an 
infinitesimal  change  in  the  water  level  in  the  stand  pipe  is  related  to  the 
flow  through  the  soil  in  an  infir i  tesimal  interval  of  time  (as  given  by 

Darcy's  law)  in  order  to  set  up  a  differential  equation;  during  this 
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infinitesimal  interval  of  time,  one  may  assume  that  steady  state  flow 
conditions  are  prevalent.  The  values  obtained  from  a  solution  of  such  a 
differential  equation  may  not  be  reliable  if  the  permeability  varies  with 
the  pressure;  such  variations  with  pressure  have  been  observed  in  the 
case  of  flow  of  gases  (Scheidegger  1957),  but  the  writer  is  not  aware  of 
any  recorded  evidence  of  significant  variations  in  the  case  of  flow  of 
liquids.  Nonetheless,  it  is  desirable  to  keep  the  total  change  of  pressure 
head  during  any  one  experiment  as  small  as  possible. 

Both  constant  head  and  variable  head  permeability  tests  can  be 
run  with  either  undisturbed  or  remolded  specimens.  In  either  case,  it 
is  essential  to  have  complete  saturation  of  the  specimen  since  the  per¬ 
meability  has  been  found  to  be  very  sensitive  to  changes  in  the  degree 
of  saturation  (Klute,  1952).  In  undisturbed  specimens,  it  is  seldom 
possible  to  have  full  saturation  for  the  following  reasons: 

1.  Due  to  natural  resistance  to  wetting,  some  air  is  entrapped  on  the 
surface  of  the  grains  during  the  formation  of  the  deposit;  2.  If  the 
soil  is  organic,  the  pore  water  may  contain  dissolved  gases  which  will 
come  out  of  solution  when  the  overburden  pressure  is  released.  With 
remolded  samples,  a  high  degree  of  saturation  is  usually  attainable  by 
first  evacuating  the  dry  soil  and  then  saturating  it  with  distilled  de- 
aired  water. 

The  capillary  permeability  test  can  be  conducted  on  remolded 
specimens  only  since  the  soil  must  be  in  a  dry  state  at  the  start  of  the 


v-  X  J. 


- 


uz 


.  -x- 

V 


x-  .  J.ij 


199 


experiment.  Since  the  pores  of  the  specimen  cannot  be  evacuated  prior 
to  the  test  (one  end  of  the  specimen  must  be  in  free  communication  with 
the  atmosphere),  some  air  is  always  entrapped  in  the  voids  during  the 
test;  consequently,  complete  saturation  of  the  specimen  by  capillarity 
is  never  realized  in  the  test  and  the  value  of  the  permeability  obtained 
from  such  tests  should  be  interpreted  as  valid  for  the  particular  degree 
of  saturation  obtained  in  the  test.  In  the  actual  test,  the  rate  of  advance 
of  the  wetting  front  in  the  specimen  is  observed  under  two  different 
pressure  conditions  in  order  to  obtain  two  simultaneous  equation  in¬ 
volving  the  two  unknown  quantities - permeability  and  another 

(fictitious)  quantity  called  'effective  capillary  head1;  this  necessitates 
an  abrupt  change  in  the  pressure  head  causing  flow.  Depending  upon 
the  magnitude  of  the  change  in  the  pressure  head  between  the  two  stages, 
there  will  be  changes  in  the  volume  of  the  entrapped  air,  thus  giving 
rise  to  corresponding  changes  in  the  degree  of  saturation.  For  this 
reason,  considerable  skill  is  required  in  setting  up  the  test  conditions 
in  order  to  obtain  reliable  results. 

Estimation  of  the  value  of  the  permeability  of  a  soil  from  its 
grain  size  or  similar  properties  is  not  an  accurate  method,  but  the 
accuracy  so  obtained  may  be  sufficient  for  many  practical  cases.  This 
method  is  suitable  for  coarse  grained  soils  only,  that  do  not  possess  a 
marked  degree  of  "structure" 

Estimation  of  the  permeability  from  the  consolidation  character 
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istics  of  the  soil  is  primarily  based  on  the  validity  of  the  Terzaghi 
theory  of  consolidation.  It  is  well  known  that  the  theory  is  based  on 
many  simplifying  assumptions,  and  "while  the  theory  is  extremely 
successful  in  determining  the  magnitude  of  settlements,  predicted  and 
observed  settlement  rates  -  where  the  permeability  comes  into  play  - 
often  differ  considerably"  (Schmid,  1957).  Thus  the  indirect  deter¬ 
mination  of  the  permeability  from  consolidation  test  data  can  be  regarded 
as  semi- quantitative  only.  Here  again,  except  in  problems  related  with 
fundamental  research  on  permeability,  the  degree  of  accuracy  obtained 
from  this  method  is  sufficient  for  almost  all  practical  cases. 

Field  tests  for  permeability  include  the  well-known  pumping 
test  and  the  ground  water  velocity  tests.  Both  types  of  tests  are  rather 
expensive  and  as  such  are  limited  for  use  in  major  projects  only.  In 
fine-grained  soils,  a  long  time  may  be  required  for  the  attainment  of 
steady  state  conditions.  However  the  results  obtained  from  such  tests 
are  far  more  reliable  than  the  other  laboratory  methods;  moreover, 
the  field  tests  yield  the  average  permeability  of  a  large  mass  of  soil 
deposit  and  "the  results  are  representative  of  the  average  properties 
of  the  soil  in  its  natural  state",  (Taylor,  1958). 

B.Z  Various  theories  of  permeability: 

On  account  of  the  importance  of  the  permeability  in  so  many 
fields  of  science,  numerous  attempts  have  been  made  to  correlate  the 
permeability  with  other  physical  properties  of  the  porous  medium. 
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Such  attempts  (at  correlation  of  the  permeability  with  the  porosity  or  the 
pore  structure  or  the  grain  size  distribution  or  similar  properties  of 
the  porous  medium)  have  never  been  successful  (Scheidegger  1957),  and 
the  empirical  correlations  claimed  to  have  been  'found1  are,  at  best, 
valid  over  a  very  narrow  range  of  values  only* 

Having  failed  to  obtain  a  unique  correlation  between  permeab¬ 
ility  and  some  other  physical  property  of  the  porous  medium,  many  in¬ 
vestigators  attemped  to  idealize  the  pore  space  of  a  porous  medium  and 
describe  the  permeability  as  a  function  of  some  physical  parameter  of 
the  idealized  model.  The  advantage  of  this  approach  is  that,  by  choosing 
the  proper  type  of  model,  the  problem  can  be  treated  mathematically. 
One  set  of  such  models  consist  of  a  bundle  of  straight  capillary  tubes, 
arranged  in  various  ways,  and  the  mathematical  treatment  stems  from 
the  well-know  Poiseuille's  law  of  flow  through  straight  circular  capil¬ 
lary  tubes.  This  idealization  is  hardly  justifiable,  since  the  pores  of  a 
natural  porous  medium  do  not  even  remotely  resemble  straight  circular 
capillary  tubes. 

Another  type  of  models  are  based  on  the  assumption  that  the 
pores  are  analogous  to  a  series  of  a  channels.  It  will  be  recalled  that 
the  modified  coefficient  of  permeability  has  the  dimensions  of  a  length 
squared  (Eq.2.6).  This  length,  characteristic  of  the  porous  medium, 
is  taken  as  the  equivalent  of  the  'hydraulic  radius1  of  the  porous 
medium.  A  possible  measure  of  the  hydraulic  radius  of  a  porous 
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medium  will  be  the  ratio  of  the  volume  to  the  surface  area  of  the  pore 
space.  Slitcher  (1899)  considered  a  bed  of  spherical  particles  as 
equivalent  to  the  porous  medium  and  arrived  at  the  equation 


K 


N2  dZ 


96  (1-  n  ) 


CO 


where,  N  =  the  ratio  of  the  void  area  to  the  total  area  in  a  plane 
section,  when  the  plane  passes  through  the  centres  of 
adjacent  spheres, 

d  =  the  diameter  of  the  spheres,  and 
n  =  the  porosity. 

Slitcher' s  work  has  been  subsequently  refined  and  extended  by  Graton 
and  Fraser  (1935)  and  others,  but  the  general  approach  may  be  critic¬ 
ised  on  the  grounds  that  for  a  natural  porous  medium,  a  unique  value  for 
the  'diameter  of  the  sphere'  cannot  be  assigned;  further  in  a  bed  of 
spherical  particles  of  equal  diameter,  there  can  be  no  interconnected 
voids  and  consequently  no  flow  through  the  pores  (Scheidegger,  1957). 

Kozeny  replaced  the  porous  medium  by  an  assemblage  of 
channels  of  various  cross-sections  but  of  a  definite  length,  solved  the 
Navier-Sto  ke '  s  equations  simultaneously  for  all  channels  passing 
through  a  cross-section  normal  to  the  general  direction  of  flow,  and 
expressed  the  permeability  as  a  function  of  the  specific  surface  of  the 
porous  medium  (Scheidegger,  1957).  The  final  expression  obtained  by 


Kozeny  is  of  the  form 
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K  -  c.n.  3/ S2  (2) 

where,  c  =  a  dimensionless  constant,  called  the  'Kozeny  constant1, 
n  =  the  porosity  of  the  medium,  and 
s  =  the  specific  surface  of  the  porous  medium. 

The  Kozeny  equation  has  been  modified  by  a  large  number  of 
people,  foremost  among  them  being  Carman  (1937),  who  set: 

K  =  n3  /  5  scz  (1-n2)  (3) 

where  sQ  is  the  'surface  exposed  to  the  liquid  per  unit  volume  of  the 
solid  (and  not  the  porous)  material.  By  implication,  the  value  of  the 
Kozeny  constant  is  1/5.  Equation  (9),  called  the  Kozeny- Carman 
equation  has  been  extensively  used  for  the  determination  of  the  surface 
area  of  industrial  powders  (Carman  1939). 

Yet  another  approach  to  a  physical  charaterization  of  the 
permeability  may  called  the  drag  theory  of  permeability  (Scheidegger 
1957).  Here,  the  walls  of  the  porous  medium  are  assumed  to  be 
obstacles  to  an  otherwise  straight  flow  of  a  viscous  fluid.  Using 
Navier  -  Stoke1  s  equations,  the  drag  on  each  portion  of  the  walTis 
estimated  and  the  total  drag  is  set  equal  to  the  resistivity  (p/K,  as 
per  Darcy's  law)  of  the  porous  medium  to  flow  of  fluids.  This  approach 
has  been  found  to  give  reliable  results  in  highly  porous  media  such  as 
textile  fibres,  where  each  fibre  can  be  regarded  as  solitary  within  the 
fluid.  (Scheidegger,  loc.  cit.). 


The  problem  has  been  tackled  from  the  standpoint  of  statistical 
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mechanics,  assuming  the  porous  medium  to  be  intrinsically  disordered. 
Childs  and  George  (1948)  were  probably  the  first  to  imply  such  an 
approach  and  Scheidegger  (1957)  has  given  a  comprehensive  outline  of 
the  method.  For  the  details  of  this  approach,  Scheidegger' s  monograph 
may  be  consulted. 

Recently,  Schmid  (1957)  has  applied  the  concept  of  a  stationary 
boundary  layer  and  of  ineffective  porosity  to  the  problem  of  percolation 
of  water  through  soils  and  has  shown  that  there  is  a  linear  relationship 
between  permeability  and  porosity.  His  final  expression  is  of  the  form: 

k  =  (V  32  u)  .  d2E  (n-n0) 

where,  k  =  the  Darcy's  coefficient  of  permeability  in  velocity  units, 
n  =  the  total  porosity 

nQ  =  the  effective  porosity,  and  the  other  symbols  are  as 
already  defined. 

The  symbol  'D^1  is  defined  as  the  'effective  diameter  of  the  pore'  and 
the  method  of  determining  the  same  is  not  described.  It  may  be  men¬ 
tioned  that  the  experimental  results  reported  by  the  author  seem  to  agree 
well  with  the  theoretical  considerations,  in  that  the  plot  of  permeability 
against  porosity  turns  out  to  be  linear  except  when  the  total  porosity  is 
approximately  equal  to  the  effective  porosity,  where  'additional 
phenomena  occur  which  need  further  clarification'. 

Thus,  there  are  so  many  theories,  purportedly  aimed  at 
ascribing  a  physical  significance  to  the  permeability  of  porous  media. 
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The  writer  fails  to  detect  any  practical  usefulness  of  such  theories 
except  the  indirect  determination  of  the  surface  area  of  industrial 
powders.  It  is  decidedly  easier  to  determine  the  permeability  by 
direct  experiments  than  to  measure  the  so-called  physical  parameters  - 
tortuosity  factor,  effective  diameter  of  grains,  effective  pore  size  and 
so  on  -  and  use  them  in  complex  formulae  of  doubtful  reliability.  It 
would  seem  that  the  time  spent  in  formulating  such  vague  theories  may 
well  be  profitably  utilized  in  perfecting  the  technique  of  experimental 
determination  of  the  permeability  and  of  the  associated  problems  of 
sampling  and  sample  preparation. 

B.3  Factors  influencing  the  permeability  of  soils: 

As  already  mentioned,  the  nature  of  the  definition  of  the 
permeability  as  adopted  in  civil  engineering  practice  renders  the  same 
as  dependent  upon  the  properties  of  the  soil  as  well  as  those  of  the  per¬ 
colating  liquid.  The  two  most  important  properties  of  the  percolating 
liquid  that  are  of  interest,  are  the  unit  weight  and  viscosity  of  the 
liquid.  In  the  case  of  water  -  which  is  the  only  liquid  of  interest  for  a 
civil  engineer  -  the  unit  weight  may  be  regarded  as  constant  for  most 
practical  purposes,  and  so,  variations  of  its  viscosity  alone  need  be 
considered.  By  far  the  most  important  physical  agency  that  influences 
the  viscosity  is  the  temperature  of  the  liquid,  the  viscosity  decreasing 
with  increasing  temperatures.  This  suggests  a  relation  of  the  form: 
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which  has  been  found  to  be  valid  for  coarse-grained  soils.  In  clays, 
temperature  seems  to  have  a  greater  influence  on  the  viscosity  of  the 
pore  water  than  in  the  case  of  coarser  soils  (Terzaghi  and  Peck,  1958 
ed.  ,  pp.  47).  When  the  interstitial  distance  between  particles  become 
smaller  and  smaller,  the  viscosity  of  the  water  in  these  interstices  in¬ 
creases  from  that  of  the  free  water  at  the  same  temperature;  Terzaghi 
gives  the  following  empirical  equation  for  the  increase  in  viscosity  in 
narrow  slits: 

u,  ,  6.02  x  icr42  „  ,  2.42  xlC43 

•  1+ - “I -  tol  + - ^ -  <6> 

where,  is  the  viscosity  of  water  at  25  C.  ,  and 

\x  is  the  viscosity  of  water  in  a  slit  of  width  2S  at  the  same 
temperature  (reported  by  Bull  and  Moyer,  1936). 

Thus,  in  fine  grained  clay  soils  and  soil  colloids,  Eq.  (11)  cannot  be 
assumed  to  hold. 

As  regards  the  properties  of  the  porous  medium  that  influence 
the  permeability,  it  is  obvious  that  the  geometry  of  the  pore  space 
should  be  the  foremost  factor  to  be  considered.  Unfortunately,  there  is 
no  unique  definition  of  the  'pore  geometry'  of  a  porous  medium.  Such  a 
definition  should  include:  the  size,  the  shape,  the  number  and  orienta¬ 
tion  of  the  particles,  the  amount  of  sealed  pores,  and  the  shape  of  the 
continuous  voids.  Then  it  is  apparent  that  the  permeability  should  be  a 
function  of  each  one  of  these  variables.  It  may  be  recalled  that  many 
attempts  at  empirical  correlation  between  permeability  and  some  other 
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physical  property  of  the  porous  medium  have  proved  to  be  unsuccessful; 
the  various  physical  properties  -  porosity,  grain  size,  void  ratio, 
'average  pore  size1,  'average  grain  size'  and  so  on  -  are  but  individual 
manifestations  of  a  composite  property,  namely,  the  pore  geometry; 
the  failure  of  the  above-mentioned  attempts  may  be  due  to  the  fact  that 
only  one  of  the  physical  properties  was  considered  as  a  variable  in  any 
one  problem  without  paying  due  regard  to  the  fact  that  more  than  one 
variable  are  required  to  characterize  the  pore  geometry.  Such  a  comp¬ 
lete  characterization  of  the  pore  geometry  appears  impossible  since  the 
constituent  properties  are  interrelated  in  a  complex  and  as  yet  unknown 
manne  r. 

For  a  very  comprehensive  discussion  of  the  factors  influencing 
the  permeability  of  soils,  the  reader  is  referred  to  Chap.  VI  of  Taylor's 
text  book  on  Fundamentals  of  Soil  Mechanics  (1958). 
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APPENDIX  C 


THERMODYNAMICS  APPLIED  TO  SOIL  WATER  SYSTEMS 

C.l.  Fundamental  laws: 

Classical  thermodynamics  is  based  on  two  fundamental  laws 
governing  transformation  and  conservation  of  energy.  The  first  law 
of  thermodynamics  -  the  law  of  conservation  of  energy  -  may  be  stated 
quantitatively  as  follows:  The  increase  in  internal  energy  of  a  sub¬ 
stance  during  any  transformation  equals  the  amount  of  heat  taken  in, 
less  the  work  done  by  the  substance;  or  algebraically, 

de  =  dq  -  dw  (1) 

Thus,  the  first  law  specifies  quantitatively,  the  relation  that  must 
exist  during  any  interchange  of  work  and  heat. 

While  mechanical  work  can  always  be  transformed  into  heat, 
the  reverse  process  of  transforming  heat  into  mechanical  work  is  sub¬ 
ject  to  further  restrictions;  these  restrictions  are  laid  down  by  the 
second  law  of  thermodynamics:  "No  self-acting  device,  unaided  by  an 
external  agency,  can  transfer  heat  from  a  body  of  lower  temperature 
to  one  of  higher  temperature"  (Edlefson  and  Anderson,  1943).  Alter¬ 
natively,  the  law  may  be  restated  as  follows:  Work  can  never  be  pro¬ 
duced  in  a  reversible  cycle  of  changes  operating  in  surroundings,  all  at 
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the  same  temperature.  The  former  statement  specifies  the  direction 
in  which  heat  will  flow  to  do  work,  while  the  latter  statement  specifies 
the  condition  under  which  no  work  can  be  done. 

C.2.  Definitions  of  term  used: 

C.2.1.  Reversible  operations: 

A  thermodynamically  reversible  process  is  one  that  is  carried 
out  so  slowly,  that  the  system  is  in  continuous  equilibrium  with  its  sur¬ 
roundings.  The  process  can  be  made  to  proceed  in  the  opposite  direction 
by  an  infinitesimal  change  in  anyone  of  the  conditions  governing  the  state 
of  the  substance  and  the  direction  of  the  process  -  such  as  temperature 
or  pressure.  Since  the  forces  responsible  for  the  process  are  of  in¬ 
finitesimal  magnitude,  the  process  never  reaches  equilibrium.  All 
natural  processes  are  irreversible  to  a  larger  or  smaller  extent,  and 
thus  the  thermodynamically  reversible  process  is  only  an  ideal  case. 

C.2.2.  Entropy: 

The  quantity  'entropy'  is  defined  as  the  ratio  of  the  amount  of 
heat  taken  in  by  a  substance,  to  the  temperature  (in  the  Kelvin  Scale)  at 
which  it  is  taken  in;  in  symbolic  form, 

a  dT  Aq 

ds  =  ~y~  or  as  =  - (Z) 

By  convention,  '  the  amount  of  heat'  is  taken  as  positive  when  heat  flows 
into  the  substance,  and  negative  when  it  flows  out  of  the  substance  into 
the  surroundings.  It  can  be  shown  that  the  total  change  in  entropy  of  a 


substance  in  any  reversible  cycle  of  operations  is  always  zero  and  that 
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the  same  is  always  greater  than  zero  in  any  irreversible  cycle  of 
changes;  the  entropy  of  a  substance  undergoing  irreversible  processes 
always  tends  to  a  maximum  positive  value  (see  Edlefson  and  Anderson, 
1943  for  a  rigorous  proof  of  these  statements). 

C.2.3.  Free  energy: 

The  "free  energy"  of  a  substance  is  defined  as  the  difference 
between  the  total  energy  of  the  substance  and  the  amount  of  unavailable 
energy;  in  other  words,  it  is  the  amount  of  energy  readily  available 
for  doing  useful  work.  Stated  algebraically, 

f  =  e+Pv»Ts  (3) 

where,  f  is  the  free  energy  of  the  substance, 

e  is  the  internal  energy  of  the  substance, 

P  is  the  external  pressure  acting  on  the  substance 
v  is  the  specific  volume  of  the  substance, 

T  is  the  absolute  temperature,  and 
s  is  the  change  in  entropy  of  the  system. 

In  a  perfectly  reversible  process,  the  change  in  entropy  is  zero,  and  the 
free  energy  becomes  identical  with  the  total  energy  of  the  substance;  as 
the  degree  of  irreversibility  increases,  entropy  also  increases,  and  the 
free  energy  decreases.  It  is  to  be  noted  that  an  irreversible  process  is 
characterized  by  a  negative  change  in  free  energy  or  a  positive  change 
in  entropy. 

From  the  second  law  of  thermodynamics  -  which  specifies  the 
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direction  in  which  work  may  be  performed  -  and  from  the  observed 
fact  that  free  water  tends  to  move  into  a  partially  saturated  soil 
spontaneously,  we  may  conclude  that  the  free  energy  of  water  in  a 
partially  saturated  soil  is  less  than  that  of  pure  water.  By  convention, 
the  free  energy  of  pure  water  outside  the  zone  of  influence  of  any  ad¬ 
sorptive  force  field  (excepting  that  due  to  gravity)  is  taken  as  zero, 
and  this  forms  the  datum  of  free  energy  for  soil  moisture  studies;  then, 
the  free  energy  of  soil- water  in  a  partially  saturated  soil  is  always 
negative;  it  reaches  a  maximum  value  of  zero  when  the  soil  is  comp¬ 
letely  saturated.  This  negative  quantity  has  been  variously  called  soil 
suction,  pressure  deficiency,  soil  water  tension,  capillary  potential 
and  so  on;  by  common  usage,  the  negative  sign  is  usually  omitted.  In 
the  following  sections,  the  term  'free  energy'  will  be  used  consistent 
with  the  practice  in  thermodynamic  analyses. 

C.3.  Characteristics  of  free  energy: 

C.3.1.  Free  energy  of  several  phases  in  equilibrium: 

In  the  preceding  paragraph,  it  was  seen  that  the  free  energy  of 
soil  water  is  less  than  that  of  pure,  free  water  and  that  free  water  tends 
to  move  into  a  partially  saturated  soil  spontaneously.  Generalizing,  we 
may  state  that  any  system  containing  two  or  more  substances  in  co¬ 
existence,  cannot  remain  in  equilibrium  unless  the  free  energies  of  the 
individual  substances  are  equalto  each  other.  Thus,  in  a  partially 
saturated  soil  in  a  state  of  thermodynamic  equilibrium,  the  free  energy 
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of  the  soil  water  is  equal  to  that  of  the  water  vapour  in  contact  with  it; 
and,  in  a  partially  frozen  soil,  the  free  energy  of  the  unfrozen  water  is 
equal  to  that  of  the  ice  formed  as  well  as  that  of  the  water  vapour  in  con¬ 
tact  with  it.  This  offers  a  convenient  way  of  determining  the  free  energy 
of  soil  water  by  measuring  the  pressure  of  water  vapour  in  equilibrium 
with  it. 


C.3.2.  Relation  between  vapour  pressure  and  free  energy: 
By  definition. 


f  s  e  +  P.v  -  T.s.  (3) 

In  differential  form, 

df  =  de  +  Pdv+v  dP  -  T  ds  -  sdT  (4) 

According  to  the  first  law  of  thermodynamics, 

de  r  dq  -  dw 

Separating  the  total  work  done,  dw,  into  two  component  parts  -  a  work 
of  expansion  against  a  constant  pressure,  P,  equal  to  P.dv,  and  a 
mechanical  work  done,  dwm,  Eq.  (1)  becomes 

de  B  dq  -  P.dv  -  dwm  (5) 

Further,  by  definition, 

dq 

ds  -  p  Or,  dq  =  T  ds  (Eq.  2) 

Substituting  these  into  equation  (4)  and  simplifying, 

dgr  vdP-sdT-dwm  (6) 

If  the  temperature  remains  constant,  and  if  no  mechanical  work  is  per¬ 
formed  in  the  process,  Eq.  (6)  reduces  to 


df  =  v  d  P 


(7) 
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Assuming  that  the  behaviour  of  water  vapour  can  be  approximated  by 
the  ideal  gas  equation, 


P.v  =  RT  (8) 

Eq.  (7)  may  be  rewritten  as 

df  =  RT  (9) 

P 


Integrating  between  two  phases  A  and  B, 

f,  -  f  -  RT  In  PA 

A  B  TIT 

Letting  the  subscript  B  refer  to  the  datum,  f 

equals  the  saturated  water  vapour  pressure, 

fA  =  RTlnPA 

P 


(10) 

g  equals  zero,  and 
designated  by  P  . 

(11) 


o 

where  in  f^  is  the  free  energy  of  the  soil  water  vapour, 
P^  is  the  pressure  of  soil  water  vapour,  and 
R  is  the  universal  gas  constant. 


P 

B 

Then, 


C.3.3  Effect  of  temperature  on  free  energy: 

From  Eq.  (6),  assuming  that  no  mechanical  work  is  performed, 
we  obtain,  for  conditions  of  constant  pressure, 

df  ;  -  s  d  T  (12) 


or, 


( 


3f 


) 


(13) 


ST  '  P 

Thus,  the  rate  of  change  of  free  energy  of  a  substance,  due  to  change 


in  temperature,  equals  the  negative  of  the  entropy  of  the  substance. 
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C.3.4.  Effect  of  hydrostatic  pressure  on  free  energy: 

This  effect  is  readily  seen  by  rewriting  Eq.  (7)  as 

3f 

— —  r  v  (14) 

ap 

Thus,  the  rate  of  change  of  free  energy  due  to  changes  in  pressure 
equals  the  specific  volume  of  the  substance.  In  particular,  if  we  have 
pure  water  and  ice  in  co-existence,  the  following  numerical  values 
apply: 

Bf 

-  z  1.000,  for  water,  and 

dP 

d  f 

-  -  1.090,  for  ice. 

dP 

An  increase  in  pressure  causes  a  greater  increase  in  the  free  energy  of 
ice  than  in  that  of  water  on  account  of  the  different  specific  volumes. 

As  the  free  energy  tends  to  be  a  minimum,  the  mass  of  ice  will  spontan¬ 
eously  melt  if  the  pressure  is  increased,  since,  by  doing  so  the  rate  of 
decrease  in  free  energy  is  kept  smaller.  On  the  other  hand,  the  rate 
of  decrease  of  free  energy  of  ice,  with  decreasing  pressure  is  greater 
than  that  of  water;  hence,  the  mass  of  water  will  freeze  spontaneously, 
if  the  pressure  is  reduced,  again  for  the  same  reason  as  above. 

C.3.5.  Effect  of  dissolved  salts  on  free  energy  of  a  solution: 

Rewriting  Eq.  (7)  in  terms  of  infinitesimal  quantities, 

Af=£B-fA  =jpB  vdP  (15> 
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Here  phase  A  refers  to  the  solution  while  phase  B  refers  to  the  pure 
solvent  and  the  above  equation  gives  the  free  energy  of  the  solution  with 
respect  to  that  of  the  pure  solvent.  In  our  particular  case,  the  'pure 
solvent1  is  the  stress-free,  pure  water  and  the  'solution'  is  the  soil-water 
containing  dissolved  salts.  If  the  external  pressure  acting  on  both  the 
pure  solvent  and  the  solution  is  the  same,  Pg  may  be  set  equal  to  zero, 
while  P^  equals  the  osmotic  pressure,  (  <^P0),  due  to  the  presence  of 
dissolved  salts.  Assuming  the  specific  volume  of  the  solvent  to  be 
constant  (the  compressibilities  of  normal  liquids  are  very  low),  equation 
(15)  reduces  to  r° 

Af  =  f  -  fA  -vjd? 

-  -v(AP0)  (16) 


Then, 


B(aP0)  ”  (17) 

or,  the  rate  of  change  of  free  energy  of  a  solvent  due  to  changes  in  the 
osmotic  pressure  of  the  solution,  equals  the  negative  of  the  specific 
volume  of  the  solvent.  Since  the  specific  volume  is  always  a  positive 
quantity,  it  is  seen  that  the  free  energy  always  decreases  with  increasing 
osmotic  pressure,  that  is,  with  increasing  quantity  of  dissolved  materials. 
In  contrast,  an  increase  in  hydrostatic  pressure  increases  the  free 
energy  of  the  substance. 


C.4.  Freezing  point  depression  of  soil  water: 

C.4.1.  Preliminary  considerations: 

It  has  been  found  (Bouyoucos  1917)  that  "some  of  the  water  in  the 


i.: 


hi.'.. 


: Ix 


L 


I 

\ 


i) 


\ 


..  :  > 


X- 


X  -r 


.■  > 


<,  - 


216 


soil  refused  to  freeze,  and  the  amount  varied  with  the  class  of  soil, 
degree  of  supercooling,  etc."  Subsequent  investigations  by  Bouyoncos 
and  others  have  revealed  that  the  freezing  point  of  water  in  soils  is  lower 
than  that  of  pure  water  under  the  same  conditions  of  pressure  and  that 
the  extent  of  depression  is  related  to  the  energy  state  of  the  soil  water. 
Nowadays,  the  freezing  point  depression  (FPD)  is  accepted  as  one  of  the 
common  means  of  estimating  the  free  energy  of  soil  moisture  (see 
Schofield  1935,  Bodman  and  Day  1938,  Schofield  and  da  Costa  1938).  The 
thermodynamic  interrelationships  between  the  free  energy  and  the  FPD 
have  been  lucidly  presented  by  Edlefson  and  Anderson  (1943).  In  the 
following  sections,  a  brief  outline  of  these  relationships  are  given;  for 
detailed  information,  the  reader  is  referred  to  the  paper  of  Edlefson  and 
Anderson. 

C.4.2.  Relation  between  FPD  and  hydrostatic  pressure: 

We  will  consider  the  case  of  a  pure  liquid  remaining  in  thermo¬ 
dynamic  equilibrium  with  its  solid  phase  under  atmospheric  pressure 
conditions;  the  temperature  of  the  system  will  then  be  freezing  point  of 
the  liquid  and  the  free  energies  of  the  two  phases  will  beequal  to  each 
other.  Due  to  an  increase  in  the  external  pressure  on  the  system,  the 
equilibrium  is  disturbed,  and  a  new  state  of  equilibrium  can  be  set  up, 
only  if  the  temperature  of  the  system  is  suitably  altered.  Then  the 
final  free  energies  of  the  two  phases  will  also  be  equal,  since  they  are 
again  in  equilibrium  with  each  other.  Thus,  the  starting  point  is  the 
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statement  that  the  initial  free  energies  and  the  changes  in  the  free 

energies  of  the  two  phases  are  separately  equal,  that  is, 

f  -  f  and  df  -  df 
1  "  2  1  2 


Substituting  for  df  from  equation  (6), 
mechanical  work  is  done  in  the  process. 


1 


dp  -  s^  dT  =  vo  dp  - 


Simplifying, 


dP  so  -  s 


1 _  I  As 


dT  v„  -  v 
2  1 


AV 


and  stipulating  that  no 

s2  dT  (18) 

(19) 


By  definition, 


T 


From  the  first  law  of  thermodynamics, 

de  =  dq  -  P.  dv  -  dw  (Eq.5) 

m 

and,  since  no  mechanical  work  is  performed,  this  may  be  rewritten  as 

dq  -  de  -V-P.dv  (20) 

-  d  (  e  +Pv)  -  dh  (21) 

In  equation  (21),  the  quantity,  h  t  e  +  P.v  ,  is  defined  as  the  'Heat 
content'  of  the  substance.  If  we  are  dealing  with  changes  of  state  of 
anyone  substance,  the  heat  content  will  become  the  commonly  known 
'latent  heat  of  fusion'  or  the  'latent  heat  of  condensation'  of  the  substance. 
From  Eq.  (21)  and  the  definition  of  entropy, 
s 2  -  -  as  -  Ah 

T"  (22) 
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Substituting  in  Eq.  (19), 


dP  - 


dT  av.  T 


(23) 


or. 


22.  P  -  Ah 


aT  av .  T 


(24) 


Thus,  the  relation  between  any  change  in  external  (hydrostatic)  pressure 
and  the  corresponding  change  in  the  equilibrium  temperature  of  a  system 
is  given  by  Eq.  (23);  in  this  equation, 


£>h  is  the  change  in  the  heat  content  of  the  substance, 
av  is  the  change  in  the  specific  volume  of  the  substance, 
T  is  the  absolute  equilibrium  temperature  (initial), 
a  P  is  the  change  in  the  external  pressure,  and 
aT  is  the  change  in  the  equilibrium  temperature. 


As  a  particular  case,  when  water  freezes  to  form  ice,  its  volume  in¬ 
creases,  and  so  av  is  positive;  but  an  amount  of  heat  equal  to  the  latent 
heat  of  fusion  of  ice  is  given  out  in  the  process,  and  so  Ah  is  negative. 
Thus,  Eq.  (23)  predicts  that  the  equilibrium  temperature  of  an  ice-water 
mixture  is  lowered  by  increasing  the  hydrostatic  pressure  on  the  system. 


Combining  equation  (23)  with  equation  (7),  we  obtain, 


v.  Ash 


f 


(24) 


which  gives  the  change  in  free  energy  of  a  substance  due  to  a  change  in 


the  hydrostatic  pressure,  as  measured  by  the  FPD. 
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C.4.3.  Relation  between  FPD  and  osmotic  pressure; 

Here  again,  we  start  with  the  same  initial  and  final  conditions, 
that  is,  f^  -  f^;  and  df^  -  df^, 

wherein,  is  the  initial  free  energy  of  a  pure  solvent  in  thermodynamic 
equilibrium  with  its  solid  phase,  whose  free  energy  is  f  .  The  addition 

Lu 

of  some  soluble  salts  to  the  solvent  decreases  its  free  energy,  as  given 
by  Eq.  (17).  Let  us  designate  the  increase  in  osmotic  pressure  of  the 
solvent  due  to  the  presence  of  dissolved  salts  by^Po. 

The  specific  free  energy  of  the  solid  phase  depends  only  on  the 
temperature,  since  the  external  pressure  conditions  are  assumed  to 
remain  unaltered;  consequently,  we  may  write, 

df  =  df9 

2  (  _ 2 _  )  dT  (25) 

ST  P 

and,  in  view  of  Eq.  (13), 

df  z  -  s  dT  (26) 


The  specific  free  energy  of  the  solution,  on  the  other  hand, 
depends  both  upon  the  temperature  and  the  osmotic  pressure;  hence, 

dfl  “  (  3£1  )  dT  3fx 


dT  Po 


•(- 


apo  t 


)  dPo  (27) 


Substituting  from  Eq.  (13)  and  (17), 


dfx  =  -  si  DT  -  vx  dPo 

df.  -  df  ,  we  also  have: 

-  s2  dT  =  -  s1  dT  -  v  dPo 


(28) 


Since 
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Or, 


Hence , 


or, 


Since 


we  have 


d  Po 
dT 


dh 
v.  T 


(29) 

(30) 


dPo  =  - 


dh 

vTT 


.  dT 


^Po  =  "  ^h,i,  .  aT  (31) 

v.  T 

f  -  -  v.  aPo  Eq.  (17), 


(32) 


which  expresses  the  change  in  free  energy  of  a  solution  with  respect  to 
that  of  the  pure  solvent,  in  terms  of  the  FPD  caused  by  the  presence  of 
dissolved  salts. 


C.5.  Methods  of  determining  free  energy  of  soil  water: 

Several  methods  have  been  evolved  for  the  determination  of  the 
so-called  capillary  potential  of  soil  moisture,  and  all  these  methods  may 
be  readily  used  for  the  evaluation  of  the  free  energy  of  the  soil  moisture. 
For  a  brief  description  of  the  principles  involved  in  the  various  methods, 
the  reader  is  referred  to  Chap.  II  of  this  dissertation;  a  comprehensive 
evaluation  of  the  merits  and  demerits  of  these  various  methods  has  been 
recently  published  (Croney  et.  al.  1952). 

C.6.  Movement  of  soil  moisture: 


In  the  preceding  sections,  some  of  the  characteristics  of  the 
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thermodynamic  function  "free  energy"  and  its  interrelationships  with 
other  physical  quantities  such  as  temperature,  hydrostatic  pressure  and 
osmotic  pressure  were  outlined.  In  this  closing  section,  the  manner  in 
which  these  relationships  may  be  utilized  to  predict  movement  of  water 
through  soils  will  be  outlined. 

Under  isothermal  conditions,  if  two  sections  of  moist  soil  in 
contact  with  each  other  do  not  have  the  same  free  energy,  moisture  will 
tend  to  migrate  spontaneously,  so  as  to  equalize  the  energies.  The 
direction  of  migration  is  specified  by  the  second  law  of  thermodynamics: 
it  will  be  in  the  direction  of  decreasing  free  energy.  An  incidental  point 
is  that  the  final  free  energy  of  the  system  as  a  whole  will  have  decreased 
by  a  certain  amount  since  the  change  is  spontaneous  and  not  perfectly 
reversible. 

In  analogy  with  the  "potential  concept"  of  driving  forces,  the 
negative  of  the  gradient  of  the  total  specific  free  energy  of  the  soil 
moisture  is  the  force  per  unit  mass  of  water  causing  movement.  Incor¬ 
porating  this  in  a  Darcy  type  equation, 

v  -  U  v(A-fs)  (33) 

where  v  is  the  steady  state  velocity  of  flow,  and 

k  is  the  proportionality  factor,  the  permeability. 

Combining  this  equation  with  the  continuity  requirement,  an  equation  rel¬ 
ating  the  moisture  content  of  the  soil  and  the  gradient  for  free  energy  may 
be  set  up;  the  solution  of  this  equation  for  particular  boundary  and  initial 
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conditions  of  the  problem  will  give  the  mode  of  distribution  of  water 
within  the  flow  region. 

There  is  one  difficulty  yet  to  be  overcome  in  this  analysis: 
the  final  equation  involves  moisture  density,  time,  the  rate  of  change  of 
free  energy,  and  the  permeability,  all  of  which  are  variables.  In  order 
to  solve  the  problem  completely,  four  independent  equations  are  required, 
or  else  it  is  necessary  to  eliminate  two  of  the  variables  by  some  other 
means.  Many  investigators  have  assumed  some  single -valued  relation¬ 
ships  between  moisture  density  on  the  one  hand,  and  the  permeability  and 
the  capillary  potential  (analogous  to  free  energy)  on  the  other  hand  (see 
Chap.  II).  On  account  of  the  lack  of  perfect  reversibility,  these  assum¬ 
ptions  are  not  strictly  justified.  Under  these  circumstances,  there  is  a 
great  need  for  expressing  thermodynamically  irreversible  functional 
relationships  in  terms  of  easily  measurable  physical  parameter;  once 
such  parameters  have  established,  we  would  have  taken  a  big  step  forward 


in  the  solution  of  this  complex  problem. 
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